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PREFACE 


The  research  report  contained  herein  was  supported  by  the 
Avionics  Laboratory,  Wright-Patterson  Air  Force  Base,  Ohio,  under 
Contract  No.  33615-77-C-1002,  “Electron  Device  Contact  Studies." 
Research  described  in  Sections  I -VI  was  performed  at  The  Ohio  State 
University,  Department  of  Electrical  Engineering.  Research  covered  in 
Section  VII  was  performed  at  The  Ohio  State  University,  Department  of 
Physics.  This  report  covers  a  period  from  November  1,  1979  to  June  30, 
1981. 


The  main  objective  of  this  research  program  is  to  investigate 
the  fabrication  of  low  resistance  contacts  to  gallium  arsenide  (GaAs) 
and  to  examine  their  behavior  with  respect  to  theoretical  models.  This 
report  covers  the  background  theory  used  in  the  modeling  of  the  con¬ 
tacts  and  gives  the  processing  steps  developed  for  the  fabrication  of 
Au-GaAs  contacts. 
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SECTION  I 
INTRODUCTION 

1.1  Purpose  and  motivation  for  this  research 

There  has  been  considerable  recent  interest  in  using  GaAs  as  the 
basis  for  high  speed  electronic  device  applications.  As  with  any 
practical  device,  making  reliable  contacts  to  the  operating  areas  of 
the  device  is  of  utmost  importance.  An  ideal  contact,  more 
specifically  an  ideal  ohmic  contact,  is  one  which  allows  current  to 
flow  in  either  direction  without  presenting  to  the  carriers  any 
resistance  at  the  interface.  If  the  interface  is  not  ideal,  the  con¬ 
tact  can  still  be  considered  ohmic  if  the  resistance  is  small  and  the 
potential  drop  developed  across  the  interface  is  essentially  linear 
with  respect  to  the  current. 

In  practice,  a  potential  barrier  is  nearly  always  present  at  the 
interface  and  presents  an  impediment  to  the  current  flow  preventing 
the  contact  from  being  ideally  ohmic.  Often  the  contact  is  highly 
non-ohmic  or  rectifying  as  a  result  of  a  non-linear  current-voltage 
(I-V)  characteristic.  The  object  of  contact  studies  is  to  determine 
where  contact  non-linearities  arise  and  how  they  affect  contact 
properties  and  to  use  that  knowledge  to  minimize  the  effects  of  the 
barrier  and  to  reduce  interface  resistance  to  as  small  a  value  as 
possible. 

As  will  be  discussed  in  detail  later  on,  contacts  which  present 
little  resistance  to  current  through  the  interface  are  referred  to  as 
“ohmic  contacts"  and  whatever  resistance  they  do  have  can  be  related  to 
the  contact  area  by  defining  a  quantity, "speci fic  contact  resistivity", 
Rq.  The  units  are  usually  in  J3-cm^  and,  for  a  contact  of  given  area, 
the  resistance  the  contact  presents  to  current  flow  can  be  easily 
determined.  Ideally  Rc  is  zero,  but  practically  R^  is  very  small, 

-10”®  n-cm^,  in  typical  contacts  being  made  today. 

The  difficulty  is  that  when  most  metals  are  applied  to  make  a 
contact  to  a  semiconductor  a  potential  barrier  forms  at  the  interface. 
Commonly  this  potential  barrier  causes  the  current  to  be  a  non-linear 
function  of  applied  bias  voltage  and  an  R^  can  be  defined  only  in  a 
piece-wise  sense.  The  usual  assumption  is  that  ohmic  contacts  in  the 
presence  of  such  a  barrier  are  made  possible  by  quantum  mechanical 
tunneling  through  the  barrier.  In  such  a  case  the  barrier  ceases  to  be 
significant  and  low  values  of  R^  cay  result. 
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Tunneling  is  initiated  in  a  metal -semi conductor  (MS)  contact  by 
doping  the  semiconductor  area  beneath  the  contact  to  a  high  concentra¬ 
tion.  This  high  doping  causes  the  formation  of  a  very  thin  space 
charge  region.  If  the  doping  is  high  enough  to  reduce  the  space  cnarge 
region  to  a  thickness  of  a  few  lattice  constants,  tunneling  can  occur. 

The  literature  on  ohmic  contacts  to  GaAs  is  quite  extensive  and 
the  amount  of  supporting  and  related  literature  on  contacts  to  other 
semiconductors,  especially  Si,  is  even  larger.  Several  particularly 
useful  reviews  have  been  given  by  Rideout  (1),  Yoder  (2),  and  Eckhardt 
(3)  with  an  introductory  textbook -monograph  having  been  written  by 
Rhoderick  (4). 

Research  on  the  theory  of  metal -semiconductor  (MS)  interfaces 
and  their  current  transport  characteristics  is  usually  considered  to 
have  begun  with  the  work  of  Mott  (5)  and  Schottky  (6).  The  diffusion 
theory  of  Schottky  and  the  thermionic  emission  theory  of  Bethe  (7)  were 
subsequently  combined  into  the  modern  thermionic  emission-diffusion 
theory  by  Crowell  and  Sze  (8).  A  good  explanation  of  the  essentials  of 
all  the  above  can  be  found  in  Sze  (9). 

The  above  research  was  concerned  with  MS  contacts  in  which 
emission  over  the  barrier  is  the  dominant  mechanism.  Ohmic  contacts, 
on  the  other  hands,  are  related  to  carriers  tunneling  through  interface 
potential  barriers,  A  classic  reference  for  this  type  of  tunneling  is 
given  by  Nordheim  (10)  .  The  most  extensive  theoretical  work  on  the 
I-V  characteristics  of  MS  tunneling  contacts  is  found  in  the  papers  of 
Padovani  and  Stratton  (11),  Crowell  and  Rideout  (12)  and  Chang  and  Sze 
(13).  Padovani  (14)  gives  a  good  summary,  with  additions  and 
corrections,  of  the  work  of  Padovani  and  Stratton.  Using  the  above,  Yu 
(15)  formulated  the  basic  theoretical  development  of  the  specific 
contact  resistivity  of  tunneling  MS  contacts.  Change,  et.  al.  (16) 
also  performed  a  numerical  analysis  for  based  on  Chang  and  Sze  (13). 

The  initial  experimental  work  on  developing  alloyed  type  ohmic 
contacts  to  GaAs  was  performed  by  Cox  and  Strack  (17).  Since  that 
time,  various  "recipes"  for  making  ohmic  contacts  have  been  examined. 
The  usual  technique  utilizes  the  alloying  of  a  Au-Ge  contact  mixture 
into  the  GaAs  surface  (18-22).  The  Ge  is  believed  to  move  into  the 
GaAs  surface  thereby  creating  a  highly  doped  region  which  forms  a  thin 
space  charge  layer  through  which  electrons  can  tunnel.  Usually  the 
alloying  is  performed  thermally  but  laser  annealing  and  electron  beam 
annealing  (3)  are  receiving  much  current  interest.  Ion  implantation 
has  also  been  tried  (23). 

Even  though  alloying  results  in  good  ohmic  contacts  with  low  R^;, 
many  problems  associated  with  alloying  are  present.  Some  arise  from 
the  non-uniform  way  in  which  the  alloying  occurs  (19).  However,  the 
most  serious  problems  arise  from  the  redistribution  of  the  Au,  Ge,  and 
Ga  during  the  alloying  (21,  24-26).  This  redistribution  can  occur  to 


quite  a  large  depth  into  the  semiconductor  and  has  severe  implications 
if  alloying  is  used  for  contacts  to  the  thin  electronic  structures 
employed  in  advanced  integrated  circuits.  This  problem  has  initiated 
much  recent  research  into  the  possibility  of  making  ohmic  contacts  via 
a  non-alloying  approach  (27-30).  The  research  described  in  this  report 
deals  in  part  with  the  possibility  of  fabricating  a  non-alloyed,  low 
R(;,  contact  to  GaAs  using  a  Sn  diffusion. 

One  large  problem  with  experimental  research  up  to  this  point  is 
that  beyond  giving  "recipes",  little  work  has  been  done  in  trying  to 
determine  whether  the  low  values  of  Rc  are  in  fact  due  to  tunneling 
through  the  barrier.  There  are  other  possibile  schemes  whereby  low 
R(.'s  might  be  obtained,  especially  since  the  alloying  step  can 
drastically  alter  the  morphology  of  the  interface.  Part  of  the  problem 
lies  in  the  fact  that  after  fabrication  of  a  low  R,-  contact  via 
alloying,  the  data  necessary  to  relate  contact  performance,  namely  Rc 
and  current-voltage  (I-V)  response,  to  theory  is  difficult  to  obtain. 
The  data  needed  are  primarily  the  potential  barrier  height,  ij>B.  at  the 
interface  and  the  semiconductor  doping,  Nn*  beneath  the  contact.  The 
value  of  No  especially  is  almost  impossible  to  determine  after  the 
alloyed  contact  has  been  fabricated. 

Models  have  been  developed  in  which  the  tunneling  I-V  response 
of  MS  contacts  can  be  theoretically  determined.  The  main  thrust  of 
this  research  was  to  attempt  to  relate  experimental  performance  of  some 
non-alloyed  Au-GaAs  contacts  to  those  models  to  see  whether  they  can 
accurately  predict  contact  behavior.  If  so,  then  contact  performance 
can  indeed  be  assumed  to  be  due  to  tunneling  phenomena  and,  more  im¬ 
portantly,  the  use  of  such  models  in  predicting  contact  performance 
under  various  conditions  is  justified.  The  application  of  an  accurate 
model  in  the  attempt  to  fabricate  contacts  that  are  even  more  ideally 
ohmic  would  allow  a  somewhat  more  scientific  approach  to  be  taken 
rather  than  the  "hit  and  miss  recipe"  approach  now  often  used. 


1.2  Fabrication  procedure 

In  this  research,  since  a  comparison  to  theoretical  models  was 
desired,  the  fabrication  process  needed  to  be  one  which  allowed  the 
important  parameters  of  (fg  and  Nq  to  be  determ  .led  accurately.  This 
meant  that  the  normal  approach  of  making  contacts  via  alloying  could 
not  be  used.  The  method  decided  upon  was  based  on  highly  doping  the 
GaAs  surface  so  that  carrier  tunneling  was  expected  but  doing  so  by 
diffusing  the  dopant  into  the  surface  first  and  then  applying  the  con¬ 
tact  metal.  This  method  of  fabrication  was  a  two-step  process  instead 
of  a  single  step  alloying  but,  of  utmost  importance,  it  allowed  the 
determination  of  the  Nq  parameter  before  the  contact  was  actually  made, 
and  therefore  allowed  for  accurate  comparisons  between  theoretical 
predictions  and  analysis.  It  was,  of  course,  necessary  that  the 
application  of  the  metal  should  not  affect  the  doping  at  depths 
critical  for  tunneling. 


The  diffusion  of  Sn  into  a  p-type  GaAs  wafer  was  used  to  create 
an  n-type  surface  layer  with  a  high  value  of  surface  dopant 
concentration,  C5,  and  since  Cs  is  large,  a  thin  space  charge  region  is 
formed  and  tunneling  can  occur.  The  theoretical  treatment  is  actually 
based  on  the  assumption  of  uniform  dopant  distribution,  Nq,  beneath  the 
contact.  However,  even  though  the  actual  dopant  distribution  has  a 
definite  profile,  with  C5  indicating  dopant  concentration  only  at  the 
surface,  the  value  of  C5  can  be  substituted  for  Nq  in  the  theoretical 
equations  with  little  error.  This  is  done  in  this  research,  but  both 
C5  and  N0  symbols  will  still  be  used  since  it  is  desirable  at  times  to 
keep  clear  the  designation  between  actual  dopant  profile  and  theoret¬ 
ical  assumption. 

The  choice  of  a  p-type  substrate  was  based  on  the  desire  to  form 
a  thin,  isolated  layer  for  conduction,  allowing  easy  determination  of 
the  data  needed  to  calculate  0$  and  Rc. 

1.3  Data  collection  and  interpretation 

To  determine  the  value  of  Cs  for  the  diffused  layer,  a  van  der 
Pauw  measurement  was  performed  and  a  value  of  diffused  layer  sheet 
resistivity  was  obtained.  Junction  lapping  and  staining  was  then  used 
to  determine  layer  thickness.  From  the  above  data  a  value  of  average 
resistivity  for  the  layer,  p,  could  be  determined.  An  "Irvin  curve",  a 
theoretical  curve  relating  p  of  the  diffused  n-type  layer  to  the  value 
of  Cs,  was  developed  and  used  with  the  measured  p  to  determine  an 
experimental  value  for  Cs* 

To  develop  the  value  of  418,  curve  fitting  of  experimental 
contact  I-V  data  to  theoretical  equations  for  the  contact  I-V  response 
was  performed.  The  curve  fitting  yeilds  both  4)8  and  Cs  parameters  and 
allowed  a  check  on  the  value  of  Cj  as  determined  from  the  p  measure¬ 
ments. 


Using  the  values  of  4)3  and  Cs  as  obtained  above,  theoretical 
values  of  Rc  were  calculated.  As  a  more  complete  check  on  the  applic¬ 
ability  of  the  model,  the  variation  of  Rc  with  temperature  was 
predicted  theoretically  and  then  compared  to  experimental  Rc  vs.  T 
data. 


The  above  model  and  procedures  were  also  applied  to  the  analysis 
of  annealed-al loyed  Au-Sn  diffused  layer  contacts  and  also  to  annealed- 
alloyed  Au-Ge  contacts.  Values  of  4>b  and  Cs  were  calculated  and 
compared  as  a  function  of  annealing  time  and  temperature.  This  part  of 
the  research  program  was  performed  to  see  whether  the  model  used  in  the 
analysis  of  non-alloyed  contacts  could  also  be  usefully  applied  to 
understanding  what  changes  alloyed  contacts  undergo  as  they  are 
fabricated  and  why  these  changes  result  in  the  contact  eventually 
exhibiting  improved  ohmic  behavior. 
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SECTION  II 


THEORETICAL  BACKGROUND 


2.1  Ohmic  contacts 

A  discussion  of  contact  performance  and  ohmic  contacts  must 
first  clarify  the  meaning  of  the  word  "ohmic".  To  define  the  term  it 
is  necessary  to  consider  the  structure  of  the  metal -semiconductor  (MS) 
i nterface. 

The  joining  of  a  metal  and  a  semiconductor  causes  the  formation 
of  a  potential  barrier  at  the  interface.  Figure  2.1  shows  the  band 
diagram  of  a  MS  contact  on  a  moderately  doped  semiconductor.  This  type 
of  barrier  is  commonly  referred  to  as  a  Schottky  barrier.  In  the 
figure,  E^  and  Ey  represent  semiconductor  conduction  and  valence  bands, 
respectively.  The  Fermi  level  energy  is  Ep.  The  potential  barrier 
height  is  given  by  41B  and  the  space  charge  region  width  is  W. 

Classically  ((>8  is  determined  from: 

i'B  =  -^M  -  Xs  ,  (2.1) 

where  (t>M  is  the  metal  work  function  and  Xs  is  the  electron  affinity  of 
the  semiconductor  (63).  If  >  <{15,  where  (p$  is  the  semiconductor  work 
function,  the  barrier  of  Fig.  2.1  results.  This  type  of  barrier 
presents  an  impediment  to  current  flow  and  even  for  relatively  low 
values  of  current,  rather  large  voltage  drops  appear  across  the 
interface.  If,  at  this  point,  a  piecewise  linear  approximation  of  the 
contact  I-V  response  is  made,  these  large  voltage  crops  result  in  a 
large  contact  resistance. 

I’f  't'M  <  'I’S*  ‘('8  become  zero  or  even  negative.  For  such  a 
case  the  carriers  crossing  the  junction  see  essentially  no  resistance 
since  no  barrier  exists.  Therefore,  the  potential  drop  across  the 
junction  is  very  small,  ideally  going  to  zero.  Such  a  contact  has  zero 
contact  resistance  and  is  the  classical  definition  of  an  "ohmic" 
contact.  (For  the  above  discussion,  a  MS  contact  to  an  n-type 
semiconductor  was  illustrated.  For  a  p-type  semiconductor  similar 
barriers  to  hole  flow  will  result  but  for  the  opposite  conditions  of 
and  (|)s.  Since  that  treatment  is  analogous,  and  since  n-type  samples 
were  used  in  this  research,  only  n-type  examples  will  be  discussed 
further. ) 
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metal 


N-TYPE  SEMICONDUCTOR 


For  device  fabrication,  the  most  desirable  contacts  are  ohmic 
contacts  as  the  word  is  defined  above.  However,  for  most  practical 
combinations  of  metals  and  semiconductors  the  (ti|v|,  Xs  relation  is  in  a 
direction  to  produce  a  barrier.  Moreover,  other  effects  such  as 
surface  states  act  to  pin  the  Fermi  level  at  a  fixed  position  resulting 
in  a  (j(B  independent  of  iti|«|  or  Xs  (31-35).  For  such  cases  the  formation 
of  a  Schottky  barrier  is  inevitable.  However,  ohmic  contacts  can  still 
result  from  other  conduction  mechanisms,  but  the  word  "ohmic"  takes  on 
a  slightly  different  meaning. 

In  the  most  general  case  of  contacts  to  devices,  the  goal  is  to 
have  the  voltage  drop  across  the  contact  Interface  much  smaller  than 
the  voltage  drop  across  the  bulk  region  of  the  device.  In  this  way  the 
contact  properties  or  contact  I-V  response  do  not  affect  overall  device 
performance.  In  addition,  the  I-V  response  of  the  contact  should  be 
essentially  resistive  or  linear.  However,  the  description  of  "ohmic" 
now  becomes  a  relative  one,  being  compared  to  the  bulk  resistance  of 
the  device  itself.  This  viewpoint  is  useful  and  for  many  practical 
contacts  is  perfectly  valid  in  characterizing  them  as  being  ohmic. 

2.2  Contact  resistivity 

A  more  quantitative  way  to  describe  contact  performance  as  it 
relates  to  ohmicity  is  in  terms  of  contact  resistivity,  R^,  which  is 
usually  given  in  units  of  Sl-cm^.  The  definition  of  Rq  implies  that  the 

effects  of  the  contact  interface  on  carrier  flow  can  be  lumped  into  a 

distributed  resistance.  If  the  I-V  response  of  the  interface  is  fairly 
linear,  the  value  of  Rc  is  applicable  over  a  wide  range  of  contact 
voltage  drops.  If  the  I-V  response  is  non-linear,  then  the  R^  value  is 
strictly  applicable  at  only  one  particular  value  of  I  and  v.  In  most 

cases,  the  value  of  R^  is  fairly  low,  making  the  voltages  drop  across 

the  contact  small,  and  typical  contact  I-V  response  over  this  region  of 
voltage  drops  is  fairly  linear. 

Lower  values  of  R^  are,  of  course,  desirable.  When  contact 
current  flow  is  perpendicular  to  the  interface,  the  actual  resistance 
of  a  contact  of  given  area.  A,  can  be  determined  from: 

Rc 

f^CONT  =  -  ♦  (2.2) 

A 

The  voltage  drop  across  the  contact  can  then  be  calculated  from: 

VCONT  “  IRcONT  •  (2.3) 

For  current  flow  in  a  planar  contact,  or  other  structure  where  current 
flow  is  not  all  perpendicular  to  the  interface,  equation  (2.3)  usually 
cannot  be  applied  directly.  In  such  cases  current  crowding  effects  are 
present  and  a  different  method  must  be  used  to  relate  the  value  of 
VconT  values  of  I  and  R^  (36-37).  However,  even  in  such  cases 

the  concept  of  Rc  is  valid  and  is  the  best  way  to  characterize  contact 
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performance  at  the  interface  level.  Therefore,  in  current 
terminology,  “ohmic"  contacts  refers  to  contacts  with  sufficiently 
small  values.  Values  of  in  the  range  ofl0“3n-cm2  to  are 

those  most  often  referred  to  in  the  literature  as  being  ohmic  (17-23, 
27-31). 

2.3  Practical  ohmic  contacts 

In  section  2.1  it  was  stated  that  for  most  cases  of  MS  contacts 
a  positive  will  be  formed.  For  Au-n-type  GaAs  contacts,  other 
workers  have  measured  a  <()b  i  n  the  range  of  0.8  to  1.0  eV,  with  0.9  eV 
being  the  value  most  often  reported  (31-35).  The  value  of  (()b  is  usual¬ 
ly  observed  to  be  constant  and  independent  of  GaAs  doping  level 
implying  that  the  Fermi  level  is  pinned  due  to  the  action  of  surface 
states. 

A  simple  model  explaining  Fermi  level  pinning  as  being  due  to 
surface  states  has  been  developed  by  Bardeen  (38).  There  is  still  some 
question  as  to  whether  the  surface  states  are  intrinsic  or  whether  they 
are  induced  by  the  deposition  of  the  contact  metal  but,  the  important 
point  is  that  when  an  Au  contact  is  made  to  GaAs  a  Schottky  barrier 
with  ijiB  “  expected.  Subsequently,  such  a  barrier  will 

result  in  a  poor  contact  unless  it  is  thin  enough  to  allow  carrier 
tunnel i ng. 

For  a  MS  contact  with  a  barrier  height  ipB  »  the  width  of  the 
insulating  space  charge  layer  beneath  the  contact  is  given  by: 

-  1/2 

w  =  —  (Vbi-  —  )  .  (2.4) 

q^D  q 


where  e  =  semiconductor  permittivity  =  CpCo 

q  =  electronic  charge 

N[)  =  semiconductor  donor  density 

k  =  Boltzmann's  constant 

T  =  temperature 

Vbi  =  built  in  voltage  (barrier  height  as  seen  by  electrons 
in  conduction  band  of  the  sani conductor) 

Vbi  =  4>b  +  np 

np  =  Fermi  potential  =  Ep  -  Ec 
kT 

or  if  HB  >>  np  ,  -r—  as  is  the  usual  case. 
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1/2 


(2.5) 


For  moderate  dopings,  W  is  too  large  to  allow  electrons  to 
tunnel  and  any  flow  must  be  over  the  barrier.  Such  a  contact  exhibits 
rectification  and  the  I-V  response  of  the  contact  is  given  by  the 
familiar  diode  equation: 


J 


♦  b 

A  exp  ( - ) 

kT 


qV 

exp  ( -  )  -1 

kT 


1/2 


(2.6) 


where 

A  =  Richardson's  constant  =  (2.7) 

h3 

h  =  Plank 's  constant 

m*  =  electron  effective  mass 


For  contacts  with  a  given  I-V  response,  the  value  of  specific 
contact  resistivity  Rc  is  defined,  in  the  limit  as  contact  voltage  goes 
to  zero,  as: 

(2.8) 

v+o 

or  equivalent,  since  the  I-V  response  is  single  valued  in  J  for  all  V, 


-1 

(2.9) 
v+o 

Applying  the  definition  of  R^  to  the  diode  equation,  applicable  for 
semiconductors  with  moderate  doping,  results  in  values  for  R^  on  the 
order  of  ki lo-ohms-cm^.  The  incremental  resistance  is,  of  course,  in 
parallel  with  a  depletion  layer  capacitance  and  the  result  is  clearly 
not  a  useful  contact  characteristic. 

However,  if  the  doping,  Nq  .  is  large,  the  band  diagram  as  seen 
in  Figure  2.2  is  present.  The  barrier  slope  is  quite  steep  and  W  is 
very  small.  For  values  of  W  on  the  order  of  a  few  lattice  spacings, 
electrons  can  quantum  mechanically  tunnel  through  the  barrier  under  an 
applied  voltage  and  give  rise  to  a  current.  The  tunneling  currents  can 
be  quite  large  even  for  a  small  contact  voltage  drop  and  therefore  the 
resistance  of  the  contact  can  be  quite  low. 
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Theoretical  equations  giving  the  I-V  response  of  tunneling 
contacts  have  been  developed  by  various  authors  such  as  Padovani  and 
Stratton  (11),  Crowell  and  Rideout  (12),  and  Change  and  Sze  (13).  Their 
results  are  all  similar,  but  the  work  of  Chang  and  Sze  is  based  on  a 
computer  numerical  solution,  using  Schroedinger' s  equation,  with  the 
least  number  of  simplifying  assumptions.  Subsequently,  Chang  and  Sze's 
results  are  not  given  in  a  closed  analytical  form  and  therefore,  a 
closed  expression  for  R^  cannot  readily  be  developed  from  their  graphs. 
It  is  easier  to  work  from  the  results  of  Padovani  and  Stratton  or 
Crowell  and  Pideout  where  I-V  equations  are  given.  The  most  compre¬ 
hensive  treatment  of  tunneling,  outlining  the  development  of  the 
tunneling  I-V  equations,  is  given  in  Padovani  (14).  Here  Padovani  has 
taken  the  earlier  work  of  himself  and  Stratton  (11),  as  well  as 
Stratton  (39),  and  presented  a  more  unified  development.  It  is  largely 
from  this  cource  that  the  equations  used  in  this  research  are  taken  and 
the  development  of  these  equations  is  reviewed  herein. 

2.4  General  aspects  of  tunneling 

The  basis  for  the  Development  of  the  I-V  response  of  a  tunneling 
contact  results  on  relating  current  flow  to  the  probability  of  electron 
transmission  through  a  thin  potential  barrier.  The  classic  work  of 
Nordheim  (10)  gives  an  expression; 


2(4iB-qV-nF) 

p  s  exp  -  -  ,  (2.10) 

3E« 

where  P  represents  the  transmission  probability  for  an  electron 
through  a  triangular  potential  barrier  of  total  height  (pB-qV-np. 

Where; 

V  =  applied  basis  voltage 
np  =  Fermi  level  referred  to 

The  term  E<»  represents  a  characteristic  tunneling  energy  and  is  given 
by: 


E.  ,  sL  _ 

2  em* 

where  =  Plank's  constant  /Zv 

Nq  =  Semiconductor  doping 
e  =  Semiconductor  permittivity 

m*  =  effective  mass  of  the  carrier  (in  this  case  an  electron). 
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The  tunneling  current  is  proportional ,  in  first  order,  to  P  and 
therefore  has  an  exponential  dependence  on  applied  forward  voltage.  The 
energy  E®  can  be  easily  shown  to  be  inversely  proportional  to  the  space 
charge  width,  W.  A  decrease  in  W,  meaning  a  thinner  barrier  and 
increased  probability  for  tunneling,  results  in  a  larger  value  for  E^, 
which  via  the  above  also  leads  to  a  larger  value  of  P.  An  increase  in 
N0  results  in  a  thinner  space  charge  region,  and  as  is  obvious  from  Eq. 
(2.11)  also  increases  E«.  Also,  since  the  electron  effective  mass,  m*, 
for  an  electron  in  GaAs  is  very  small, 

m*  =  0.068  mg 

where  (2.12) 

mo  =  9.1  X  lO"^''  kg  , 

the  value  for  E*  in  GaAs  is  much  larger  than  for  other  semiconductors 
where  m*  is  larger,  and  therefore,  for  similar  doping  levels,  electron 
tunneling  is  much  more  probable  for  barriers  made  to  GaAs.  This  is  why 
tunneling  is  believed  to  be  the  main  mechanism  giving  rise  to  low 
resistance  contacts  to  GaAs. 

Even  though  a  Schottky  barrier  can  be  approximated  to  first  order 
as  a  triangular  barrier,  the  actual  relationship  between  electron 
tunneling  and  current  flow  is  much  more  complex.  The  reason  is  that  the 
electrons  are  distributed  in  a  range  of  states  throughout  the  conduction 

band  and  the  contribution  of  all  the  electrons  must  be  considered.  The 

relationship  between  electron  tunneling  probability  and  current  flow  in 
the  X  direction  is  given  by: 

J  =  -T  o/  /  P(Ep  Py.  p2)dpy  dP2  dE^ 

(2.13) 

The  subscripts  1  and  2  refer  to  the  conductors  where  the  electrons 
originate  and  where  they  are  destined  respectively.  For  a  MS  contact 
with  an  n-type  semiconductor  under  forward  bias,  "1"  represents  the 
semiconductor  and  "2"  represents  the  metal.  The  Fi(Ei)  and  F2(E]^)  are 
Fermi  functions,  predicting  whether  or  not  an  electron  state  at  energy 
El  is  occupied.  The  transmission  probability  is  P(Ei,  py,  p^)  and  the 
integration  is  performed  over  all  energies  and  momenta  in  the  conduction 
band. 


One  large  problem  in  using  the  above  equation  lies  with  the 
formulation  of  P  for  a  given  barrier  shape.  For  any  specific  barrier 
shape,  an  analytical  solution  of  Eq.  (2.13)  using  an  extraction 
expression  for  P  is  extremely  difficult.  Instead,  P  is  usually  formu¬ 
lated  using  the  WKB  approximation  which  results  in: 

5  ^2  1/2 

P  =  exp  {-  y  /  [-(Px)^]  dx) 

n  xi 
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The  value  of  Px  is  the  x  component  of  electron  quasi -moment un  in  the 
space  charge  region.  For  a  total  momentum  of  p^  =  px^  +  py^  +  p^^,  the 
transmission  of  electrons  in  the  x  direction  is  dependent  on  Px.  Let¬ 
ting  pj^^  =  Py2  +  p^^  ,  Eq.  (2.14)  can  then  be  rewritten  as; 

?  *2  1/2 

F  =  exp  -  y  /  (Px^-P^)  dx  .  (2.15) 

n  x^ 

The  minus  sign  in  -(Px)^  in  Eq.  (2.14)  arises  from  the  WKB  approximation 
because  between  xi  and  X2,  the  classical  barrier  turning  points  where 
Px  =  0,  the  value  of  Px^  is  <0  since  the  electron  is  within  a  potential 
barrier  ())(x),  whose  energiy  is  larger  than  that  of  the  tunneling  elec¬ 
tron. 


One  large  analytical  difficulty  remains.  The  expression  for  InP, 
Eq.  (2.15),  is  still  difficult  to  use,  as  is,  in  Eq.  (2.13).  The 

simplification  used  to  overcome  this  problem  is  to  expand  InP  into  a 
Taylor  series.  This  procedure  has  been  done  by  Stratton  et.al.  (40), 
and  results  in  a  series  for  an  expansion  in  (pi)^  of 


(2.16) 


with  the  coefficients  b  and  P^^  defined  as: 


h  XI 


- 

/  p  dx 


(2.17) 


1  ^2  1 

Poi"  =  -  /  -^  dx  ,  (2.18) 

h  XI  p 

where  p  =  -p'^  (again  since  within  the  barrier  region  p^  <  0).  Since 
terms  above  px^  will  not  be  used,  K  need  not  be  given. 

The  further  simplification  of  Eq.  (2.16)  involves  expanding  b 
into  a  Taylor  series  around  the  energy  level  from  which  the  electrons 
are  tunneling.  This  will  result  in  two  different  regions  of  operation, 
field  emission  (FE)  and  thermionic  field  emission  (TFE)  which  will  be 
treated  separately  in  the  following  sections. 

2.5  Field  emission 

In  the  case  of  field  emission,  the  electrons  tunnel  from  states 
around  the  semiconductor  Fermi  level,  np.  as  shown  in  Figure  2.3.  The 
applied  forward  bias  is  V  and  (ti(x)  is  the  potential  barrier. 


Ill  liJ 


c 


Figure 


2.3  Electron  tunneling  under  field  emission  for  an 
applied  forward  bias. 
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Since  electrons  are  tunneling  around  np,  the  simplification  for 
Eq.  (2.16)  is  to  expand  b  in  terms  of  a  Taylor  series  around  the  Fermi 
level  (40).  The  series  is  given  by: 

b  =  bi  +  Ci(E-nF)  +  fl{E-nF)^  +  *•*  .  (2.19) 

where  the  coefficients  b,  ci,  and  f^  are  given  by: 


and 


bi  = 

f'F  ' 

0  ^2 

■p^  I  (P)nF  ‘^x 

Cl  = 

2 

h 

dp 

( — )nF  "^x  » 

d 

n  = 

1 

TT 

XI 

(^~^)nF  ‘^x 
dE2 

1 

1 

h 

X2  1 

/  {7)nF  dx  . 

P  2 

01 

X 

1  P 

where 

all 

integrals  are  evaluated 

coefficients  Into  Eq.  (2.19)  and  use  of  Eq. 
then  in  Eq.  (2.13)  for  J  gives: 


(2.20) 

(2.21) 

(2.22) 

(2.23) 

E  *  np*  Substitution  of  the 
(2.19)  in  Eq.  (2.16)  and 


j  =  2Tr£  (P  )2exp(-bi)  /*  CFi(E)-F2(E)]  exp  [-ci(E-nF)]  dE 

h3 


1  2ir  ~FM^  V 

X  [1-  ^  I  exp  ■( - )  dp]  (2.24) 

o  p  2 

^01 

which  is  valid  as  long  as  (42) 

-  Cl  <  (2fi)^/2  .  (2*25) 

The  next  step  in  the  analysis  is  to  choose  the  appropriate 
relationship  between  momentum,  p,  and  energy,  E,  so  the  coefficients  can 
be  evaluated.  The  assumption  used  here  is  that  of  a  parabolic  energy- 
momentum  relationship  in  all  directions  given  by: 

p2  =  2m*(E-nF)  .  (2.26) 
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For  electrons  tunneling  through  a  barrier,  (j)(x),  Eq.  (2.26)  becomes: 

=  2m*(E-nF-4>)  (2.27) 

and  within  the  barrier  (space  charge  layer)  E  >  np  but  E<(i)  ,  there¬ 
fore  using  p"^  =  -p2  as  the  quasi  momentum  in  the  barrier: 

p^  =  2m*[(|)-(  E-np)]  (2.28) 

It  is  from  Eq.  (2.28)  that  the  coefficients  pQi,  bj^,  cj,  and  f^  are 
determined.  The  ({)  term  also  includes  a  dependence  on  applied  bias  V, 

Padovani  (14)  states  that  with  the  assumption  of  a  parabolic 
energy -momentum  relationship  and  with  the  assumption  that  m*  can  be 
used  to  describe  the  electron  while  it  is  tunneling  through  the 
barrier,  for  most  values  of  applied  voltage,  energy  levels,  and  barrier 
heights: 

-  =  C^E  .  (2.29) 


This  relationship,  if  applied  in  Eq.  (2.24),  results  in: 

A  ircikT 

j  - - exp(-bi){ -  [1-exp  (-ciqV)]  -  ciqVexp(cinF)} 

(cikT)2  sin(TfCikT) 

(2.30) 

valid  when : 

(np-qV) 

exp  -  >  >  1  .  (2.31) 

kT 

The  previous  Eq.  (2.7)  gives  the  value  of  the  Richardson  constant. 

As  seen  from  Eq.  (2.31),  the  above  expression  for  J  is  valid 
only  for  the  case  of  small  applied  biases,  and  is  somewhat  less  valid 
at  room  temperature  (higher  kT)  than  for  lower  temperatures  where 
tunneling  conduction  is  usual ly- considered  to  be  dominant. 

If 

(np-qV) 

exp  -  <  <  1  (2.32) 

kT 
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as  kvould  be  the  case  for  larger  forward  biases,  then  Eq.  (2.24) 
becoines: 

A  TTClkT 

J  = - _exp(-bi){ -  -  (l+cinF)exp(-cinF)}  (2.33) 

(cikT)^  sin(itcikT) 

The  coefficient  bj,  cj,  and  fj  used  in  Eq.  (2.25)  are  the  same 
for  both  Eq.  (2.3l))  and  Eq.  (2.33),  and  must  be  determined  after  the 
shape  of  the  potential  barrier,  <^{x)  is  defined.  If  a  standard 
parabolic  Schottky  barrier  of  height  (jig  formed  on  material  with  a 
doping  of  Np  is  used  as  4i(x),  then  the  coefficients  becomes: 


Hr  (4iB-dV  + 

bi  =  {  C - ] 

Eoo  n  p 


•In  {  [- 


np 


1/2 


]  } 


Cl  =  —  In  {  [- 
E« 


hF 


Tpr 


] } 


and 


4E«nF 


'{’B 


-np) 


1/2 


(2.34) 


(2.35) 


(2.36) 


The  term  E«,  is  the  characteristic  tunneling  energy  and  is  given  by  Eq. 
(2.11).  The  applied  voltage,  V,  is  in  volts,  and  all  other  energies 
E®  .  np  are  given  in  joules. 

In  most  practical  cases  np  <<  (|>B"q'^  therefore  the 
coefficients  will  simplify  to: 


^B  -  qV 

bi  •  - 

Eoo 


(2.37) 


Cl 


4(,B-qv) 


2E,o 


(2.38) 
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and 


1 

4E»nF 


(2.39) 


f 


1 


The  above  coefficients  are  substituted  in  Eq.  (2.30)  or  Eq.  (2.33)  to 
get  the  form  of  the  I-V  response  for  a  given  range  of  bias  voltage. 

For  the  case  of  large  applied  bias,  (Eq.  2.33)  )  it  is  seen  that 
the  dominant  term  for  the  voltage  response  is  the  exp  (qV/E*)  term.  It 
is  then  possible  to  rewrite  Eq.  (2.33)  as  : 


J  =  Js  exp  (  ) 


(2.40) 


where 


Js 


A  <{iB 

- -  exp  ( - ) 

(cik)2  Eoo 


uqkT 

- } 

sin(nqkT) 


-  (1  +  cinF)exp(-ciriF)} 

(2.41) 


with  Js  representing  a  saturation  current.  Note  that  in  this  case  J 
does  not  go  to  zero  for  V  =  0.  This  is  because,  for  low  bias,  Eq. 
(2.30)  must  be  used  instead,  which  does  give  J  =  0  for  V  =  0.  It  is 
interesting  to  note  the  similarity  between  Eq.  (2.30)  and  Eq.  (2.10) 
from  Nordheim.  The  same  dependence  on  exp(qV/E  )  is  seen,  as  well  as 
the  exp  (-<{>b/E*)  dependence  which  was  incorporated  into  0$  in  the 
above. 


The  theoretical  definition 
Eq.  (2.9)  as: 


Since  Rc  is  given  in  the  limit  of 
use  to  compute  is  Eq.  (2.30). 
resulting  equation  giving  Rc  for 

Airq  <((B 

Rc  =  -  exp( - )  - 

kTsin(TrqkT)  E* 

The  coefficients  bj,  cj,  and  fj  a 


of  contact  resistivity  was  given  by 


(2.9) 


V+0  ,  the  appropriate  J-V  equation  to 
This  has  been  done  by  Yu  (15)  and  the 
FE  type  contact  is; 

Agq  (()B 

- eqp( - cihf)]  (2.42) 

(clkT)2  E« 

e  also  evaluated  at  V  =  0  and  are: 
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(2.43) 


$ 


and 


Cl  - 


In 


2£« 


1 

fl  -  - 


(2.44) 


(2.45) 


With  values  of  Nq  ^nd  ^ig  chosen,  bi,  ci,  and  fi  are  calculated  and  Eq. 
(2.25)  is  used  to  determine  whether  the  contact  is  conducting  in  the  FE 
range.  If  Eq,  (2.25)  is  not  satisfied  then  Eq.  (2.42)  cannot  strictly 
be  used. 


Chang  and  Sze  (13)  have  calculated  theoretical  I-V  curves  for  a 
MS  contact  based  on  an  equation  similar  to  Eq.  (2.13).  However,  they 
consider  a  total  response  of  both  tunneling  and  thermionic  emission 
components  together  instead  of  treating  separate  areas  of  conduction  as 
outlined  with  the  conditions  given  by  Eq.  (2.25).  They  also  solve  for 
the  transmission  coefficient,  P(E),  by  a  numerical  integration  of 
Schroedinger's  equation  within  the  barrier  region.  Their  results  are 
therefore  not  in  closed  form,  as  is  Eq.  (2.30)  and  Eq.  (2.33),  and 
cannot  be  used  to  directly  derive  an  expression  for  Rq.  However,  they 
do  make  a  comparison  between  their  numerical  solution  for  P(E)  and  a 
WKB  approximation  for  P(E)  and  show  that  if  electron  tunneling  occurs 
at  points  other  than  near  the  top  of  the  barrier,  the  value  of  P(E)  as 
predicted  by  the  WKB  approximation  is  very  close  to  the  exact  value. 
This  means  that  the  WKB  approximation  is  valid  over  most  of  the  range 
of  doping,  Nq,  and  applied  bias,  V,  used  in  this  research. 

Change,  Fang,  and  Sze  (16)  have  made  a  numerical  calculation  of 
Rc  based  on  the  computer  analysis  of  Chang  and  Sze  and  have  presented 
curves  giving  Rc  for  various  values  of  Nq  and  <(>b  as  well  as  curves  of 
Rc  vs.  T  for  MS  contacts  ot  GaAs.  Their  results,  as  expected,  are  very 
similar  to  the  values  of  Rc  as  predicted  by  Eq.  (2.42).  The  disadvan¬ 
tage  of  the  computer  calculation  is  the  necessity  of  performing  the 
tedious  numerical  integration  for  each  change  in  parameter  Ng,  ifig,  or 
T.  Therefore,  the  simplification  in  calculation  and  interpretation 
available  with  the  use  of  Eq.  (2.42)  is  worth  whatever  small  loss  of 
accuracy  results  from  using  the  WKB  approximation. 
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2.6  Thermonic  field  emission 

If,  for  a  particular  choice  of  Nq,  and  T,  Eq.  (2.25)  is  not 

satisfied,  then  the  contact  may  be  conducting  in  the  thermionic  field 
emission  (TFE)  range.  This  range  can  also  be  called  thermally  assisted 
tunneling.  Figure  2.4  illustrated  TFE  conduction.  The  electrons  are 
not  being  emitted  from  states  around  Ep,  as  they  are  in  FE,  but  rather 
from  states  higher  up  in  the  band  around  a  maximum,  E^. 

The  theoretical  calculation  of  J  for  TFE  involves  the  same 
procedure  used  in  the  base  of  FE,  except  that  instead  of  expanding  the 
Taylor  series  for  P  around  Ep,  it  is  expanded  around  E|v|.  Padpovani  and 
Stratton  (11),  working  from  Stratton  (41),  performed  such  a  calculation 
on  a  Schottky  barrier  with  the  results. 

^  np  Em  It  1/2  1/2 

J  =  -  (exp  [(_)-bM.(— )]}  (— )  [l-erf(Emfm  ]  (2.46) 

2irkT  kT  kT  fm 


where  E^  is  the  peak  of  the  energy  distribution  of  the  emitted 
electrons  and  is  given  by: 


CmkT  =  1  (2.47) 

[Note;  Padovani  and  Stratton  include  an  extra  —  factor  in  Eq.  (2.46) 
which  is  not  in  Stratton's  original  article.]  Equation  (2.46)  is 
applicable  to  forward  biases  greater  than  a  few  kT/q. 


The  calculation  of  the  coefficients  bn,,  Cn,,  and  f^  follow 
analogously  from  the  previous  equations  given  for  the  FE  case  after  an 
assumption  as  to  the  shape  of  the  emitted  electron  distribution  has 
been  made.  The  assumption  of  a  Gaussian  distribution  with  half  width: 


(ln2)^^^  ( 


1/2 


(2.48) 


was  used  by  Padovani  and  Stratton  (11)  and  for  a  forward  biased 
Schottky  barrier  the  coefficients  were  calculated  to  be  : 


1 


1/2 


E  Em 


bm  =  —  [  (4b-<1V  +  nF)*^‘'(«t)B-qV  +  np  +  En,)^^*"-  ]  , 

E  kT 


1  C  (4B-qV  +  +  np  +  Ep)^*^^  ] 


Cm  =  —  In  { 


Em 


1/2 


(2.49) 


(2.50) 
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EMISSION 


Figure  2.4  Electron  tunneling  under  thermionic-fields 
emission  for  an  applied  forward  bias. 


and 


wi  th 

(fB-qv  +  np) 

Em  ^ - 

cosh^  {—) 
kT 

E«  being  as  defined  before. 


(2.51) 


(2.52) 


The  use  of  the  above  is  also  limited  to  cases  where  is 
between  Ep  and  the  top  of  the  barrier.  This  gives  conditions: 

cikT  >  1  (2.53) 

P(Em)  <  exp(-l)  (2.54) 

where  cj  is  given  by  Eq.  (2.38)  and  P(En,)  is  the  transparency  of 
the  barrier  at  E^.  The  condition  of  Eq.  (2.54)  is  from  Murphy  and  Good 
(42).  Substituting  the  coefficients  in  Eq.  (2.46)  and  neglecting  the 
erf  term  gives: 


qV 

J  *  vJs  «P  (— )  (2.55) 

Eo 


where 

Ec. 

Eo  =  E«  coth  ( — )  (2.56) 

kT 

and 


Js 


A[TrE«(i>B-qV  +  np)]^'^^ 
- exp 

cosh2(— ) 
kT 


'IF 

kT 


4*8  ’^F 

Eo 


(2.57) 


However,  this  equation  is  still  applicable  only  for  values  of 
forward  bias,  V,  greater  than  a  few  kT/q.  It  is  clear  that  Eq.  (2.55) 
does  not  result  in  J  =  0  for  V  =  0. 


A  more  recent  development  of  the  I-V  response  of  an  MS  contact 
in  the  TFE  region  has  been  given  by  Crowell  and  Rideout  (12).  Using 
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assumptions  similar  to  the  above  for  the  shape  of  the  emitted  electron 
distribution  (Gaussian)  and  the  region  where  TFE  dominates,  they  give 
an  equation  which  is  written  as: 


J 


qV 

Js  exp  ( — ) 


qV 

-  exp  (— ) 


kT 


(2.58) 


and 


Js 


A  [nE,(4>8-qV  +  np)] 


kT 

cosh  (~) 
kT 


E» 

tanh( — ) 
kT 


1/2 


np  (t>B  +  np 

exp  (  —  -  -  )  (2.59) 

kT  £o 


This  equation  differs  from  Eq.  (2.55)  in  the  extra  [l-exp(-  |y)]  term 
in  Eq.  (2.58).  Also,  Eq.  (2.59)  has  an  extra  Ctanh(|y)]^^^  term. 
Otherwise  the  equations  are  the  same. 

However,  the  recent  presentation  by  Padovani  (14)  of  Padovani 
and  Stratton's  (11)  results  show  that  whereas  Eq,  (2.46)  is  applicable 
only  for  V  greater  than  a  few  kT/q  it  can  be  made  to  apply  at  all 

forward  biases  with  the  inclusion  of  a  [l-exp(-  ■^)]  term.  Also 

given  in  Padovani  (14)  is  a  new  equation  for  f^  ,  Eq,  (2,51),  which 
becomes: 

fn,  =  [4E«E^  (1- — )]’^  .  (2.60) 

<>8 

The  new  term,  (1  -  |^),  is  reported  to  be  a  "correction"  term  to  the 
original  Eq.  (2.51). 

If  the  new  f^,  coefficient,  Eq.  (2.60),  and  the  original  values 
of  bm.  Cm,  and  Em,  Eqs.  (2.49),  (2.50),  and  (2.52),  are  used  in  Eq. 

(2.46),  and  if  Eq.  (2.46)  is  multiplied  by  the  [(1  -  |y)]  term,  the 
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result  is  the  same  as  Eq.  (2.59).  That  is,  the  results  of  the  I-V 
response  in  the  TFE  case  as  derived  by  Crowell  and  Rideout  will  match 
those  of  Padovani  and  Stratton  if  the  above  mentioned  corrections  are 
carried  through.  The  same  assumption  used  in  the  FE  case,  that 
()>8  >>  V  +  np,  is  used  here. 

It  should  be  noted,  however,  that  even  though  Padovani  (14) 
outlines  the  above  correction  to  Padovani  and  Stratton,  he  still  uses 
Padovani  and  Stratton's  original  equation  (Eqs.  (2.55),  (2.56),  and 
(2.57)  )  in  describing  the  TFE  response.  Obviously,  the  more  correct 
form  is  that  of  Eqs.  (2.58)  and  (2.59)  and  these  equations  will  be  used 
in  this  research. 

The  equations  for  the  TFE  I-V  response  also  have  a  range  of  T, 
((>8,  and  Nq  over  which  they  are  adequate  approximations.  This  range  is 
given  by  Eq.  (2.53)  and  Eq.  (2.54).  Using  Eq.  (2.38)  for  ci,  and  Eq. 
(2.53) : 


kT  >  2E  {  In 
00 


4(<t>B-qv) 

’IF 


-1 

} 


(2.61) 


must  be  satisfied  in  order  to  have  TFE  conduction  occuring. 

The  other  limitation  to  TFE  is  Eq.  (2.54).  Padovani  shows  that 
this  limit  results  in  ; 

cosh2( — )  2(4ib  +  np  -  V 

- -  <  -  .  (2.62) 

sinh^( — ) 
kT 

If  Eq.  (2.62)  is  not  satisfied,  the  contact  I-V  properties  are  governed 
by  thermionic  emission  (TE)  rather  than  thermionic  field  emission.  In 
the  TE  range  the  I-V  response  is  that  of  a  MS  diode  with  rectifying 
properties.  This  type  of  contact  does  not  exhibit  ohmic  properties, 
and  therefore  operation  in  the  TE  range  will  not  be  considered  in  this 
research.  For  the  most  part,  with  the  values  of  (t)B»  Np  and  T  used  in 
this  research,  Eq.  (2.62)  is  always  satisfied. 

Some  additional  clari fication  must  be  made  concerning  the 
conditions  expressed  in  Eq.  (2.25)  and  Eq.  (2.53).  If  Eq.  (2.25)  is 
rewritten  slightly  it  becomes: 

cikT  <  1  -  kT(2fi)^''^  •  (^-^S) 
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Therefore,  if  c^kT  >  1  -  kT(2fi)  '  ,  the  contact  is  not  operating 
strictly  in  the  FE  range.  For  certain  values  of  T,  f]^,  and  cj^, 

1  -  kT(2fi)^/^  <  cikT  <  1  (2.64) 

is  satisfied.  The  left  hand  inequality  means  that  the  contact  is  not 
in  the  FE  range  but  at  the  same  time  the  right  hand  inequality  means 
that  the  contact  is  not  yet  strictly  into  the  TFE  range. 

The  problem  here  is  not  that  neither  FE  or  TFE  conduction  is 
occuring,  but  that  in  this  transition  range  the  Taylor  expansion  used 
in  determining  P(E)  is  not  strictly  accurate.  The  actual  type  of 
emission  consists  of  both  FE  and  TFE  occuring  simultaneously.  In  such 
a  case  the  strict  use  of  either  FE  or  TFE  equation  for  the  I-V  response 
will  lead  to  slight  errors. 

For  conditions  that  depart  only  slightly  from  Eq.  (2.25)  the 
equation  for  the  FE  I-V  response  can  be  used  as  an  approximation. 
Otherwise,  the  TFE  I-V  response  can  be  used  and  will  result  in  slight 
errors  until  the  conditions  are  such  that  Eq.  (2.53)  becomes  satisfied. 


The  determination  of  an  expression  for  in  the  TFE  range,  as 
given  in  Yu  (15),  is  again  based  on  the  use  of  Eq.  (2.4).  Using  this 
with  the  I-V  response  given  by  Eq.  (2.58)  and  ignoring  the  minor 
voltage  dependence  in  the  term  of  Eq.  (2.59),  the  value  of  R^  for 
the  TFE  range  is  given  by: 


kT  kT 

Rc  =  —  - yT? 

qA  CifE„(<l>B+nF)]^' 
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cosh(— )  exp 
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E„  kT 


(2.65) 


with  Eo  being  given  by  Eq.  (2.56).  This  equation  is  applicable  if  Eq. 
(2.61)  is  satisfied  and  can  be  used  to  give  a  fairly  close  approxima¬ 
tion  until  the  condition  given  by  Eq.  (2.25)  is  approached.  A  transi¬ 
tion  to  Eq.  (2.42)  can  then  be  made. 


2.7  Theoretical  calculation  of  Rf  based  on  FE  and  TFE  equations 

Using  Eq.  (2.42)  and  Eq.  (2.65)  with  T  =  zyz  and  as  a  para¬ 
meter,  Figure  2.5  results.  The  range  of  Nq  is  from  10^ ^/cm^  to  10^° 
/cm^.  The  range  in  ())b  of  0.7  eV  to  1.1  eV  covers  the  values  usually 
reported  in  the  literature  for  Au  contacts  on  GaAs.  The  most  often 
reported  value  of  ijig  is  •  0.9  eV.  The  curves  in  the  upper  left  portion 
of  the  figure  are  for  the  TFE  range  of  operation  and  those  in  the  lower 
right  are  for  the  FE  range.  The  dotted  line  in  the  center  represents 
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the  transition  region  where  the  TFE  and  FE  equations  for  are  not 
strictly  valid.  The  dotted  line  is  merely  a  smooth  fit  curve  joining 
the  solid  lines  together. 

It  is  evident  from  this  figure  that  for  415  =  0.9  eV,  an  Np  on 
the  order  of  10^^/cm^  can  initiate  FE  type  tunneling  but  in  order  for 
Rc  to  become  very  small,  on  the  order  of  lO'^si-cm^  ,  a  value  of  Nq 
near  10'^°/cm^  is  necessary.  Also  evident  from  the  figure  is  that  the 
value  of  Rc  is  extremely  sensitive  to  Nq.  In  the  near  10^ ^/cm^  range, 
a  factor  of  two  change  in  Np  gives  rise  to  nearly  an  order  of  magnitude 
change  in  R^.. 

Based  on  the  same  equations.  Figure  2.6  gives  the  temperature 
variation  of  R^  in  the  FE  range  with  410  =  0.9  eV  and  Nq  as  a  parameter. 
For  Rc  to  remain  in  the  FE  range  at  room  temperature,  with  4)8  =  0.9  eV, 
NG  needs  to  be  above  «•  2  x  10^ '^/cm^.  Figure  2.7  gives  the  temperature 
response  or  R^-  in  the  FE  range  for  Nq  =  6  x  10^^/cm^  with  4)3  as  a 
parameter.  It  is  evident  from  both  the  figures  that,  as  the  contact 
properties  are  pushed  more  into  the  FE  range  by  larger  Nq  of  lower  43, 
the  value  of  R^  becomes  more  independent  of  temperature.  This  is  the 
expected  characteristic  of  a  contact  where  carrier  transport  is  by 
tunneling  instead  of  by  a  thermally  assisted  emission  process. 

Figure  2.8  gives  the  temperature  response  of  Rc  in  the  TFE 
range  of  410  =  0.9  eV  and  Nq  as  a  parameter.  It  should  be  noted  that 

toward  the  upper  end  of  Nq,  (Nq  »  8  x  10^®/cin^)  the  TFE  equation  for  R^ 
is  not  strictly  valid  and  the  values  shown  in  this  region  are  the 
approximation  made  by  using  the  TFE  equation.  Figure  2.9  gives  the 
temperature  response  for  R<;  in  the  TFE  range  for  Nq  =  6  x  10*/cm^,  the 
upper  end  of  the  TFE  region,  with  410  as  a  parameter.  In  both  cases  it 
is  seen  that  since  TFE  is  a  thermally  assisted  tunneling  process,  there 
exists  a  larger  temperature  variation  in  R^  . 

The  theoretical  equations  for  R,-  under  TFE  and  FE,  and  their 
resulting  plots  of  Rc  vs.  Nq  and  Rq  vs.  T  ,  are  the  predictions  with 
which  the  experimental  data  in  this  research  will  be  compared.  If  good 
correlation  between  theoretical  prediction  and  experimental  data  is 
observed,  then  the  Schottky  model  would  appear  to  be  adequate  for 
describing  contact  properties  and  could  therefore  be  used  with  more 
confidence  in  future  contact  studies. 
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Figure  2,6  Variation  of  contact  resistivity  with  temperature 
and  doping  level.  FE  range  with  Pg  =  0.9  eV. 
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Figure  2.7 


Variation  of  contact  resistivity  with  temperature 
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and  barrier  height.  FE  range  with  Nq  =  6  x  10  /cm  . 
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Figure  2.8  Variation  of  contact  resistivity  with  temperature 
and  doping  level.  TFE  range  with  0r,  =  0.9  eV. 
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SECTION  III 


EXPERIMENTAL  METHOD 


3.1  Fabrication  and  characterization  of  Sn  doped 
surf^ace  layers 

The  contacts  investigated  in  this  research  were  fabricated  on 
n-type  layers  diffused  into  the  surface  of  a  p-type  GaAs  wafer. 

Diffused  layers  were  chosen  over  bulk  n-type  GaAs  material  because  it 
was  necessary  to  be  able  to  vary  the  concentration  of  dopant  in  the 
GaAs  to  determine  how  this  variation  affected  the  value  of  contact 
resistance.  For  tunneling  contacts,  the  n-type  layer  surface  dopant 
concentration,  Cs,  is  the  parameter  used  in  place  of  Ng  in  the 
theoretical  equations.  A  p-type  substrate  was  used  because  the 
junction  formed  between  the  n-type  region  and  substrate  would  provide 
isolation  for  the  actual  current  carrying  layer  and  would  allow  for  a 
straightforward  determination  of  Cg. 

The  fabrication  of  the  diffused  layers  utilized  a  Sn  bearing 
oxide  layer  spun  on  the  p-type  substrate  surface.  The  formation  of 
n-type  diffused  layers  in  GaAs  using  Sn  in  an  oxide  layer  source  has 
been  demonstrated  by  other  workers.  To  overcome  this  problem, 
diffusions  are  usually  performed  by  enclosing  the  GaAs  wafer  in  an 
evacuated  and  sealed  or  a  partially  sealed  quartz  box  or  ampoule 
containing  some  elemental  As  or  extra  pieces  of  GaAs.  In  such  a  closed 
environment  the  As  escaping  from  the  additional  pieces  of  As  and  GaAs 
creates  an  overpressure  of  As  which  prevents  the  escape  of  As  from  the 
wafer  surface. 

In  this  research,  open  tube  diffusions  were  tried  initially  and 
the  problems  of  As  loss  and  undesirable  surface  disruption  were 
observed.  It  was  therefore  decided  that  a  closed  chamber  type  of 
environment  was  necessary.  An  evacuated  and  sealed  ampoule  structure 
was  not  used  because  of  the  relatively  difficult  fabrication.  Instead 
a  "semi -closed  chamber"  (SCC)  type  of  vessel  was  used. 

Figure  3.1  shows  a  picture  of  the  SCC,  made  from  a  quartz 
inner-outer  joint,  the  ends  of  which  were  closed  as  in  a  test  tube. 

For  the  diffusion,  the  dopant  coated  GaAs  wafers  were  placed  inside  the 
SCC  along  with  some  pieces  of  elemental  As  and  extra  pieces  of  GaAs. 

The  two  halves  of  the  SCC  were  pushed  together  and  the  SCC  was  placed, 
into  an  open  quartz  boat  which  allowed  for  easy  insertion  and  removal 


from  the  diffusion  furnace.  During  the  diffusion  an  Ar  flow  was  main¬ 
tained  through  the  furnace. 


During  the  initial  period  of  the  diffusion,  the  As  in  the  SCC 
vaporized  and  the  increased  pressure  forced  the  SCC  apart.  This  As 
vapor  escaped  from  the  furnace  tube  and  was  exhausted  through  a 
"scavenger  box"  at  the  furnace  mouth.  After  this  occured,  the  SCC  was 
obviously  not  totally  enclosed  but,  with  only  a  slight  separation 
between  the  halves  of  the  SCC,  the  internal  atmosphere  was  As  saturated 
and  stagnant  and  this  was  enough  to  improve  the  surface  quality  very 
much  over  the  open  tube  diffusion  case.  In  a  few  cases  some  slight 
pitting  was  still  seen  but  the  surfaces  generally  remained  specular  in 
appearance. 

After  the  diffusion,  pieces  of  the  wafers  were  scribed  out  for 
lapping  and  staining.  The  test  pieces  were  hand  lapped  on  a  jig  with 
an  angle  of  -  2°.  The  junction  was  then  stained  with  a  commercial 
staining  solution  (48)  which  delineated  the  n-type  surface  layer  from 
the  p~type  bulk.  The  depth  of  the  n  layer  (the  junction  depth),  xj, 
was  measured  by  examining  the  lapped  and  stained  sample  under  a 
microscope  that  was  equipped  with  a  Watson  interferometer.  The  value 
of  Xj  was  calculated  from; 

X 

Xj  =  n(-)  (3.1) 

2 

where  X  is  the  wavelength  of  the  illumination  (in  this  case 
X  »  0,546  u),  and  n  was  the  number  of  whole  fringe  shifts  from  the 
original  surface,  over  the  lapped  edge,  to  the  delineated  junction. 

With  a  little  practice,  shifts  of  half  a  fringe  were  easily  apparent 
with  shifts  of  1/4  of  a  fringe  being  harder  to  estimate.  An  accuracy 
of  ^  1/2  a  fringe  translates  to  •>  jf  0. 14  u  accuracy  in  the  measurement 
of  Xj.  The  diffusion  process  parameters  of  time  and  diffusion 
temperature  were  chosen  to  attempt  to  realize  xj  values  of  “  lu  and 
therefore  the  accuracy  to  which  xj  could  be  determined  with  the  above 
measurement  is  on  the  order  of  15%. 

The  remaining  main  portion  of  the  diffused  wafer,  on  which  the 
actual  contact  structures  were  made,  was  then  tested  to  determine 
values  for  sheet  resistivity,  Rs,  average  resistivity,  p  ,  average  Hall 
mobility,  uhi  ^nd  average  carrier  concentration  N^.  The  measurement 
was  made  by  the  van  der  Pauw  (49)  technique.  Contact  to  the  n-type 
layer  was  made  with  small  dots  of  In-Ga  eutectic  on  the  wafer  surface 
along  the  periphery  of  the  wafer.  The  relative  ohmicity  of  these 
contacts  was  checked  by  examining  the  I-V  curve  of  the  diffused  layer 
between  all  contacts  on  a  transistor  curve  tracer.  A  straight  line 
represented  adequate  contact  ohmicity.  (Actually,  due  to  the 
potentiometric  nature  of  the  measurement  used  in  the  van  der  Pauw 
technique,  relatively  small  contact  non-ohmicity  is  unimportant.) 


The  resistance  values  obtained  from  the  van  der  Pauw 
measurements  were  used  to  determine  a  value  of  p  for  the  layer.  The 
determination  of  p  involved  the  use  of  the  value  of  Xj  determined 
earlier.  Van  der  Pauw  type  Hall  effect  measurements  were  made  to 
verify  that  the  layer  was  n-type,  and  to  determine  values  of  ph  and 
N[).  The  average  mobility  and  electron  concentration  values  are  not 
directly  useful  in  determining  Cj  but  do  give  an  order  of  magnitude 
check  on  similar  values  observed  by  other  researchers  for  GaAs.  It 
should  also  be  noted  that  gives  the  average  ionized  electron  density 
and  not  the  average  donor  dopant  density,  unless  full  ionization  can  be 
assumed. 

The  measurement  of  p  is  very  important  because  p  data  is  used 
to  determine  the  value  of  dopant  surface  concentration  Cg,  and  Cg  is 
the  most  sensitive  parameter  controlling  R^.  This  is  readily  seen  from 
Figure  2.5.  Therefore,  equally  important  is  the  way  in  which  Cj  is 
determined  from  p  data.  The  method  used  in  this  research  was  to  relate 
p  to  Cs  through  an  “Irvin"  curve  (50). 

An  “Irvin"  curve  is  a  plot  of  Cg  vs.  p  and  is  developed  by 
considering  how  p  is  related  to  Cj,  C0  and  the  shape  of  the  diffused 
dopant  profile.  Because  of  the  importance  of  the  Cg  term  in 
accurately  comparing  theoretical  predictions  of  Rc  to  experimental 
data,  the  rest  of  this  section  details  the  development  of  an  "Irvin" 
(p-Cg)  curve  for  Sn  diffusions  in  GaAs  wherein  several  important  facts 
unique  to  this  research,  are  considered. 

The  development  of  the  p-Cg  "Irvin"  curve  initially  involves 
considering  that  for  an  n-type  sample  the  resistivity  is  given  by: 

1 

-  =  qun"  (3.2) 

p 


where  p 


resistivity 


q  =  electronic  charge 
Up  =  electron  mobility 
n  =  electron  concentration 


In  a  diffused  layer  n  is  a  function  of  distance  into  the  wafer 
and  therefore  p  is  also  a  function  of  distance.  The  effect  of  the 
profile  is  to  create  an  n-type  layer  which  has  an  average  p  value. 
Indicated  by  p,  and  given  by: 


1 


p 


q 


/  ^  Wn  (N)N(x)dx 


(3.3) 
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Eq.  (3.3)  represents  the  effect  of  all  individual  p(x)  values  taken  in 
parallel  and  averaged  over  the  entire  surface  layer  from  x  =  0  to  x  = 

^j* 


The  integral  for  p  is  not  easily  calculated  because  the 
mobility,  up,  is  a  function  of  doping  and  is  therefore  also  a  function 
of  X.  The  usual  way  to  calculate  Eq.  (3.3)  is  to  use  a  computer  and 
approximate  the  integral  with  a  numerical  integration.  However,  the 
integral  can  first  be  simplified  by  removing  the  dependence  on  xj. 

The  electron  concentration,  N(x),  is  determined  from  the  dopant 
distribution. 


X 

C(x)  =  Cs  erfc  — -  (3.4) 

2  /  Dt 

where  0  =  diffusion  coefficient  for  Sn  in  GaAs 
at  the  diffusion  temperature  used 

t  =  time  of  the  diffusion 

Eq.  (3.4)  is  the  standard  profile  expected  for  a  diffusion  from  an 
unlimited  source  as  was  the  case  for  this  research.  A  simple 
derivation  of  Eq.  (3.4)  can  be  found  in  Grove  (51). 

If  a  substitution  of 

y  =  — ^ -  (3.5) 

2  ✓~Dt 

is  made  then 


dx  =  2  /  Dt  dy 

yj  =  —37 

2  /  Dt 

and  Eq.  (3.3)  becomes 

1  ‘I  ,  -yj 

-  “  ,  _ _  Un(N)N(y)  2  /  Dt  dy 

yj(2/Tt} 

or  simpl i fied: 


(3.6) 

(3.7) 


(3.8) 
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J  un(N)N(y)  dy 


(3.9) 


q 

n  ° 

where  yj  is  given  by  the  condition  that  at  the  junction  the  free 
electron  concentration  equals  the  background  concentration  of  holes  in 
the  bulk  p-type  substrate.  That  is: 

N(yj)  =  Cb  (3.10) 

The  simplification  in  Eq.  (3.9)  over  Eq.  (3.3)  is  that  for  the 
calculation,  specific  information  on  D  or  t  is  no  longer  needed  since 
the  averaging  process  removes  the  dependence  on  these  variables.  For 
any  particular  fabrication  run  in  which  0  or  t  may  vary,  Eq.  (3.10) 
must  still  be  satisfied  and  therefore  the  given  values  of  Cb  and  C5 
alone  determine  p  .  Most  importantly,  the  curve  developed  from  Eq. 

(3.9)  is  valid  for  any  particular  combination  of  D  or  t,  as  long  as  the 
assumptions  concerning  the  shape  of  the  profile  are  the  same  and  Eq. 

(3.10)  is  satisfied. 

In  the  development  of  the  p-Cj  curve  from  Eq.  (3.9)  the  value  of 
Uf,(N)  does  not  have  to  be  given  with  respect  to  x  or  y  because  in  the 
computer  calculation,  at  each  y  increment  points,  N(y)  is  calculated 
first  and  then  un(N)  is  determined.  Therefore,  a  ''S.  N  variation  is 
sufficient.  The  relationship  for  ^'’<1  N  used  in  this  calculation  was 
obtained  by  fitting  an  equation  for  pf,  in  terms  of  N  to  the 
experimental  data  for  p^  vs.  N  given  in  Sze  (9). 

(Some  important  facts  should  be  noted  concerning  the  pp  vs.  N 
graph  for  electrons  in  GaAs  as  given  in  Sze.  The  graph  is  labeled  as 
Pn  vs.  "impurity  concentration"  instead  of  "electron  or  carrier 
concentration".  However,  in  the  original  article  of  Sze  and  Irvin  (5), 
from  which  the  graph  was  taken,  the  statement  is  made  that  the  plot  is 
actually  of  pp  vs.  carrier  (electron)  concentration  for  various  bulk 
n-type  GaAs  samples.  The  authors  used  the  assumption  that  in  GaAs  at 
room  temperature  all  donor  impurities  were  ionized  and  therefore 
carrier  concentration  was  equivalent  to  impurity  concentration.) 

From  this  point,  further  development  of  the  p-Cj  curve  depends 
on  the  formulation  of  the  N(x)  term  in  Eq.  (3.9).  It  is  therefore 
necessary  to  consider  the  relationship  between  the  Sn  dopant  concentra¬ 
tion  profile,  C(x),  and  active  electron  concentration  profile,  N(x). 
Specifically  ,  an  assumption  must  be  made  as  to  the  degree  of  dopant 
ionization.  Tuck  and  Badawi  (44)  confirmed  that  the  Sn  concentration 
in  GaAs  can  be  accurately  characterized  by  an  erfc  profile  as  given  in 
Eq.  (3.4).  However,  they  also  made  some  measurements  which  implied 
that  not  all  Sn  donors  were  ionized,  meaning  that  H(x)  was  not  the  same 
as  C(x),  as  it  would  be  for  full  ionization. 
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The  problem  of  detennining  whether  or  not  all  Sn  donors  are 
ionized  can  be  analyzed  by  considering  the  physics  of  electron 
conduction.  For  this  research,  diffusions  with  rather  large  dopant 
concentrations,  in  the  10^®  to  10^°  per  cm  ®  range,  were  attempted. 

In  this  range  one  would  expect  that  the  dopant  is  not  fully  ionized. 
Therefore,  a  calculat'on  must  be  made  to  theoretically  determine  the 
Fermi  level  and  degree  of  ionization.  This  straightforward  calculation 
is  based  on  the  charge  balance  equation  (53): 

n  =  No+  «  p  -  Na-  (3.11) 

where  for  an  n-type  layer,  p  •  0  and  Na”  =  crucial 

assumption  used  at  this  point  is  that  even  at  these  very  high  levels  of 
dopping  the  Sn  donors  can  be  assumed  to  lie  at  a  discrete  energy  level, 
Eq  •  5  meV  below  £<-  (54).  The  number  of  free  electrons,  n,  can  then  be 
calculated  using  the  expression: 

2  ^  1 

"  =  Net - )  r - ^ - ST—  >^1=  (3-12) 

^  °  1  rf) 

where 

2nn*kT  212 

Nc  *  2  ( - )  (3.13) 

h2 

is  the  effective  density  of  states  in  the  conduction  band 
e  *  a  energy  referenced  to  Eq  «  0 

np  =  Fermi  level  referenced  to  E^  =  0 


The  above  equation  assumes  a  parabolic  density  of  states  vs.  energy 
distribution  for  electrons  in  the  conduction  band. 

The  ionized  donor  density  is  given  by: 

1 

Nd+  =  No  -  (3.14) 

1  +  2exp  (^) 

where  Nq  =  concentration  of  donor  atoms  . 

With  the  above,  Eq.  (3.11)  becomes: 


41 


(3.15) 

Eq.  (3.15)  is  solved  to  yield  a  value  of  Fermi  level,  np,  and  subse¬ 
quently  gives  the  value  of  electron  concentration,  n,  for  the  chosen 
conditions  of  Nq,  Cg,  Eq,  and  T  within  the  given  assumptions. 

When  the  p-Cj  (“Irvin")  curve  was  being  initially  calculated  to 
use  in  this  research  it  was  believed  that,  as  opposed  to  a  calculation 
performed  by  Galiga  (55)  in  which  the  full  ionization  assumption  of 
N(x)  =  C(x)  was  made,  partial  ionization  as  included  in  Eq.  (3.15) 
needed  to  be  considered.  Therefore,  a  computer  numerical  integration 
of  Eq.  (3.9)  was  performed  and  involved  the  following  steps; 

1.  For  a  given  Cj  and  T  the  dopant  concentration 
Nq  in  Eq.  (3.15)  was  given  by 

ND(y)  =  C(y)  =  Cs  erfc  y 

2.  At  each  integration  point,  y-j,  Eq.  (3.15)  was 
used  to  determine  np  and  subsequently  N(yi). 

The  value  of  the  Fermi -Dirac  integral  for  a 
particular  value  of  np  was  evaluated  using  the 
series  approximation  of  Battocletti  (56). 

3.  Using  the  un  N  graph  in  Sze  (9),  from  which 
a  fitted  equation  of 

un  -  .1500  [N(y)]  +  30,000 (3.16) 


was  obtained,  a  value  of  un(yi}  calculated. 

4.  The  product  of  un(yi)'N(yi)  gave  the  integrand 
value  at  yi. 

5.  The  integral  was  calculated  using  a  Simpson's 
approximation  (70)  involving  the  summation  of 
all  the  un(yi)*N(yi)  products  up  to  a  point 
where  N(yi)  »  Cg.  namely  the  junction  at  y  =  yi. 


6.  p  was  then  calculated  using  Eq.  (3.9). 
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This  sequence  of  calculations  was  performed  for  a  range  of  Cj 
and  Cg  values  and  the  graph  of  Figure  3.2  resulted.  This  graph  and  the 
preceding  calculations  were  used  to  attempt  to  determine  Cj  for  some  of 
the  experimental  wafers  but  severe  problems  became  apparent.  As  can  be 
seen  from  the  experimental  data  section  in  this  report,  experimental 
values  of  p  in  the  low  lO"^  range  were  measured.  If  the  calculations 
given  in  Figure  3.2  were  correct,  C5  values  very  much  greater  than  10^^ 
are  indicated.  The  curve  of  Cc  vs.  p  is  increasing  so  rapidly  in  the 
low  p  region  that  values  of  lCK°/cm^  or  lO^Vcm^  would  be  likely. 

Given  that  the  manufacturers  specification  for  the  dopant 
concentration  of  the  undiluted  spin  on  solution  was  only  10^°/cm^,  the 
accuracy  of  above  results  and  Figure  3.2  is  highly  doubtful.  Also,  if 
Cs  values  of  10^°/cm^  or  more  were  present,  the  GaAs  would  be  extremely 
highly  doped  and  drastic  changes  in  band  structure  could  be  expected, 
leading  to  inaccuracies  in  the  assumptions  used  to  derive  the  equations 
for  Cs  vs.  p. 

For  these  reasons  it  was  concluded  that  Figure  3.2  was  not 
correct.  The  main  di screpaancies  lay  with  the  assumption  of  partial 
ionization  of  the  Sn  dopant  at  a  discrete  level  of  Eq. 

An  explanation  of  electronic  conduction  in  heavily  doped  semi¬ 
conductors,  as  given  in  various  sources  such  as  Fistul  (57),  Madelung 
(58),  or  Blakemore  (59),  point  to  the  formation  of  an  impurity  band 
with  associated  impurity  band  conduction,  especially  for  GaAs  at  the 
dopant  levels  used  in  this  research.  The  extremely  small  value  of  Eq 
for  Sn  could  quickly  lead  to  band  tailing  such  that  a  continuous  exten¬ 
sion  of  Ec  down  through  Eq  and  into  the  band  gap  is  possible. 

As  Sze  (9)  details,  there  have  been  some  attempts  at 
quantitatively  determining  the  changes  in  the  band  structure  resulting 
from  high  doping.  These  calculations  are  extremely  difficult  in  them¬ 
selves  and  incorporating  such  calculations  into  the  calculation  of  p 
and  then  0$  would  clearly  be  intractable.  However,  as  Madelung  (58) 
states,  when  such  band  tailing  is  present,  the  effect  is  to  reduce  Eq 
to  zero  and  a  reasonably  accurate  assumption  is  to  consider  that  all 
the  donor  states  are  ionized. 

With  this  in  mind,  the  calculation  of  Cs  vs.  p  was  changed  to 
incorporate  the  assumption  of  full  ionization  of  all  Sn  donors.  This 
results  in  the  use  of  N(y)  =  C(y)  =  Cs  erfc(y)  in  Eq.  (3.9).  The  value 
of  un  still  given  by  Eq.  (3.16).  The  use  of  the  Sze  graph  for  up 
vs.  N  seemed  valid  since  it  is  an  experimental  dataplot  of  actual 
mobilities.  Regardless  of  the  mechanism  responsible  for  determining 
the  actual  ionization  level  of  donors  and  their  subsequent  mobilities, 
in  impurity  bands  or  otherwise,  Sze's  graph  gives  the  actual  value  of 
Up  measured  for  a  measured  value  of  carrier  concentration  N.  There¬ 
fore,  the  limiting  assumption  in  calculating  the  Cs  vs.  p  plot  was  not 
the  use  of  Sze's  graph  in  light  of  band  tailing  and  the  resulting 


Figure  3.2  Variation  of  p  and  for  an  erfc  profile,  n-type 
diffusion  in  GaAs.  P-type  substrate,  background 
concentration,  Cg.  Partial  ionization  of  Sn 
donor  atoms. 


change  in  conduction  specifics,  but  rather  the  assumption  that  doping  a 
sample  with  C(y)  Sn  atoms  resulted  in  N(y)  =  C(y)  free  electrons. 

When  the  full  ionization  assumption  was  used  and  the  same 
integral  calculation  procedure  was  performed  on  Eq.  (3.9),  Figure  3.3 
resulted.  Since  this  assumption  and  general  calculation  procedure  was 
the  same  as  used  by  Baliga  (55)  the  results  should  be  the  same.  A 
comparison  of  his  Cc  vs.  p  graphs  and  Figure  3.3  shows  that  the  values 
are  nearly  identical. 

As  to  the  validity  of  the  assumption  of  full  ionization, 
experimental  data  given  later  on  will  show  that  the  value  of  Cs  deter¬ 
mined  from  a  measurement  of  p  and  the  use  of  Figure  3.3,  agrees  well 
with  the  values  of  0$  determined  from  I-V  measurements  on  the  test 
samples.  The  agreement  of  these  two  Cj  values  indicates  that  the 
assumption  of  full  ionization  is  valid.  This  point  will  be  discussed 
again  in  later  sections. 


3.2  Fabrication  of  contact  structures 

The  contact  structure  used  in  this  research  is  shown  in  Figure 
3.4.  The  squares  at  both  ends  are  the  contacts  to  be  tested.  The 
interior  strips  are  voltage  pick-off  strips  which  allow  the  measurement 
of  the  IR  drop  along  the  n-type  layer.  This  structure  allowed  the 
determination  of  Rc  by  the  transfer  length  method  as  will  be  detailed 
in  the  next  section. 

The  fabrication  of  the  contact  structures  utilized  a  photolitho¬ 
graphic  lift-off  technique.  After  the  n-type  layers  had  been  diffused 
and  the  dopant  film  removed,  a  layer  of  photoresist  was  spun  on  the 
surface.  The  photomask  used  to  expose  the  photoresist  consisted  of  a 
photographic  emulsion  plate  negative  having  five  identical  contact 
structures.  This  allowed  the  fabrication  of  five  test  structures  at 
one  time  on  the  same  diffused  layer.  After  exposure  and  development  of 
the  photoresist  the  resulting  pattern  consisted  of  holes  in  the 
photoresist  layer,  exposing  the  GaAs  surface  where  contacts  were  to  be 
made. 


The  next  step  was  to  evaporate  a  layer  of  Au  onto  the  photo¬ 
resist  covered  wafer.  An  Au-GaAs  MS  contact  was  formed  in  those  areas 
where  there  were  holes  in  the  photoresist  layer. 

To  remove  the  unwanted  Au  and  to  form  the  final  contact  pattern 
of  Figure  3.4,  the  wafer  having  the  photoresist  coating  and  Au  film  was 
placed  in  a  beaker  filled  with  acetone.  The  acetone  dissolved  the 
photoresist  and  lifted  off  the  Au  not  in  contact  with  the  GaAs,  leaving 
behind  the  desired  contact  structure  pattern. 

After  the  contact  structures  were  formed,  each  separate 
structure  was  scribed  apart  to  make  an  individual  test  element.  Each 


of  these  structures  was  composed  of  the  Au  pattern  on  top  of  an  n-type 
diffused  layer,  separated  by  p-n  junction  isolation  from  the  p-type 
substrate  beneath.  The  scribed  structure  was  chosen  over  a  mesa-etched 
type  of  separation  because  of  the  ease  with  which  pieces  of  GaAs  can  be 
scribed  apart  and  the  relatively  smooth  cleave  edge  which  results. 
Etching  apart  the  individual  structures  would  have  also  have  involved 
another  masking  step.  Overall,  etching  would  have  resulted  in  more 
complicated  processing  without  any  overriding  advantage. 

The  formation  of  individual  n-type  test  islands  in  the  p-type 
substrate  was  also  attempted  but  the  process  was  not  used  because  of 
difficulties  associated  with  the  initial  oxide  layer  needed  to  mask  the 
Sn  diffusion  into  the  non-island  areas  of  the  substrate. 


3.3  Transfer  length  measurement 

The  transfer  length  measurement  (60),  also  known  as  the  Schottky 
method,  is  actually  an  extension  of  the  general  transmission  line 
method  (TLM)  (61)  of  measuring  Rj-.  In  the  transfer  length  method,  the 
actual  measured  quantity  is  the  transfer  length,  Lj  ,  which  gives 
through  the  expression 

Rc  =  RsLt^  (3.17) 

The  derivation  of  Eq.  (3.17)  and  an  understanding  of  the  physical 
meaning  of  transfer  length  can  be  based  on  the  model  of  an  MS  contact 
with  planar  current  flow  as  shown  in  Figure  3.5(a).  Figure  3.5(b) 
shows  the  variation  of  the  voltage  in  the  conducting  layer  versus 
distance  from  the  contact  in  Figure  3.5(a).  The  conducting  layer  (in 
this  case  n-type)  is  characterized  by  a  sheet  resistance,  Rj.  The 
non-ideal ly  conducting  properties  of  the  interface  are  lumped  into  the 
term  Rc  which,  in  the  model,  represents  a  distributed  resistance  of  the 
contact.  Theoretically  Rc  is  given  by  Eq.  (2.9). 

In  Figure  3.5(a),  for  purposes  of  this  initial  calculation,  the 
contact  length  is  assumed  to  be  infinite.  Practically,  this  assumption 
is  valid  as  long  as  the  Ly  value  is  much  less  than  the  contact  length. 
The  contact  edge  defines  the  x  =  0  point  with  x  <  0  representing  the 
non-contacted  diffused  n-type  layer.  The  contact  metal  is  assumed  to 
be  perfectly  conducting  and  therefore  at  a  uniform  potential,  in  this 
case  ground. 

Current  is  flowing  into  the  layer  under  the  contact  from  the 
left.  For  X  <  0  the  voltage  drop,  V(x),  with  respect  to  x  is  dependent 
only  on  the  IR  drop  along  the  Rs  of  the  layer.  Or: 

Zs  dV(x) 

I  =  -  —  -  (3.18) 

Rs  dx 
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Figure  3.5  Model  of  a  MS  contact  used  to  determine 
transfer  length,  (a)  Lumped  element 
model,  (b)  Voltage  distribution. 
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where  Zs  is  the  width  of  the  contact  structure  through  which  the 
current  flows. 


For  X  <  0,  a  plot  of  V(x)  vs.  x  yields  a  straight  line  of  slope 
=  -IRs/Zg.  Such  an  experimental  plot  of  V(x),  as  measured  from  the 
voltage  pick  off  strips,  versus  the  separation  of  the  strips,  is  called 
the  transfer  length  plot  and  the  value  of  Rj  can  be  calculated  from  the 
slope  since  I  and  Zs  are  known. 

Eq.  (3.18)  is  valid  up  to  the  contact  edge  at  x  =  0.  At  some 
point  xj  >  0,  some  current  is  lost  to  the  contact  through  R^.  In  the 
interval  dx  around  X2  this  amount  of  current  is: 

Z 

dll  - - dxV(x)  (3.19) 

Rc  x  =  xi 

and  the  current  remaining  in  the  conducting  layer  at  x  =  xi  is 

Z  dV(x) 

II  =  I(x  =  xi)  - - -  (3.20) 

“3  X  .  XI 

Taking  the  derivative  of  Eq.  (3.20) 


dl 

Z 

d2V(x) 

dx 

Rs 

dx^ 

X  =  X  2 

or  rewritten: 

dl  1 

(  X  =  XI 

=  - 

Z  d2v(x) 

—  -  dx 

Rs  dx^ 

which  must  be  equal 

to  Eq. 

(3.19).  Therefore 

(3.21) 


(3.22) 


Z  d2v{x)  Z 

-  —  -  dx  =  —  dxV(x)  (3.23) 

Rs  dx2  Rc 


1  d2v(x)  1 

—  -  -  —  V(x)  =  0  (3.24) 

Rs  dx2  Rc 

which  has  the  solution: 
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(3.25) 


V(x)  =  v(0)exp(--^} 

where  L7  =  [Rc/^s]^^^  which  is  Eq.  (3.17). 

The  slope  of  Eq.  (3.25)  at  x  =  0  is: 


(3.26) 


and  must  be  the  same  as  dV/dx  for  x  <  0. 

Therefore,  if  a  plot  of  V(x)  vs.  x  is  made  and  the  line  for 
X  <  0  is  extended  with  the  same  slope  past  x  =  0,  the  intersection  of 
this  extended  line  with  the  V(x)  =  0  axis  will  occur  at  x  =  Lj.  A 
typical  transfer  length  plot  showing  the  value  of  Lj  is  given  in  Figure 
3.6.  In  this  figure  the  right  hand  contact  is  ground  and  V(x)  is 
measured  with  respect  to  this  point.  The  voltage  at  each  pick  off 
strip  is  assumed  to  be  the  voltage  at  a  point  along  the  surface  which 
is  the  center  of  the  strip.  This  does  not  give  rise  to  large  errors 
since  the  width  of  the  strip  is  small  compared  to  their  separation. 

A  value  of  Lj  for  the  left  hand  contact,  where  the  current 
enters,  can  also  be  found.  In  this  case  Lj  is  the  difference  between 
the  contact  edge  and  the  point  on  the  extended  line  where  V(x)  is  equal 
to  the  applied  voltage,  V^pp*  Any  asymmetry  in  contact  fabrication  due 
to  processing  will  show  up  as  different  values  of  Lj. 

In  the  above  derivation  any  voltage  drop  in  the  vertical 
direction  due  to  vertical  current  flow  in  the  conducting  layer  has  been 
implicitly  neglected.  Schuldt  (36)  analyzed  the  exact  case  and  showed 
that  for  thin  layers,  as  in  the  case  here,  the  above  approximation 
leads  to  very  little  error. 

An  interesting  point  should  be  noted  concerning  the  actual 
resistance  of  the  contact,  RcONT*  A  non-zero  value  for  RcONT  gives 
rise  to  a  finite,  measurable  voltage  drop  across  the  contact,  which  is 
the  V(0)  given  by  Eq.  (3.25).  V(0)  can  be  calculated  from  Eq.  (3.26) 

and  from  the  slope  dV/dx  for  x  <  0  which  is 


(3.27) 


Substituting  gives: 


51 


V,  in  mV 


V(0) 


iRsLt 

z 


or; 

V(0)  RsLt 

>^CONT  =  -  =  - 

I  Z 


and  substituting  Eq.  (3.17)  for  Rg  gives: 


Rcont 


zlt 


(3.28) 


(3.29) 


(3.30) 


Normally,  when  current  flow  is  perpendicular  to  the  interface,  the 
contact  resistance  RcONT  determined  from  R^  as: 

Rc 

f^CONT  =  -  (3.31) 

AREA 

However,  for  the  above  case,  with  lateral  current  flow  into  a 
planar  contact,  the  apparent  contact  area  is  ZLy  and  is  independent  of 
actual  contact  length.  Therefore,  in  a  planar  contact  with  a  value  of 

Rc  giving  rise  to  a  small  Lj  .  the  current  will  crowd  into  the  edge  of 

the  contact  within  an  effective  distance  of  Lt  •  In  order  to  lower 
RCONT  such  contacts,  the  width,  Z,  can  be  made  wider,  but 

increasing  the  contact  length,  and  with  it  the  apparent  area,  will  not 

reduce  RcONT  through  Eq.  (3.31)  as  is  commonly  thought.  Fang,  et.  al. 
(37)  discuss  this  idea  further. 

According  to  the  previous  development,  if  Rc  =  0  then  Ly  =  0. 
However,  if  Rc  =  0  then  the  vertical  IR  drop  through  the  conducting 
layer  is  no  longer  negligible  in  comparison  to  the  IRc  drop  and  the 
previous  treatment  is  in  error.  Hower  (60)  and  Berger  (61)  give  a 
correction  factor  but  in  most  cases  of  non-zero  Rc  it  can  be  neglected 
with  little  resulting  error. 

For  the  previous  development  an  important  assumption  was  that 
the  contact  length  was  infinite.  Practically,  if  a  transfer  length 
plot  is  made  and  yields  a  value  of  Ly  which  is  very  much  shorter  than 
the  physical  length  of  the  contact,  then  Eq.  (3.17)  is  valid  to  use  to 
determine  Rc. 

However,  for  cases  in  which  Rc  is  not  sufficiently  small,  as  was 
the  case  for  many  contacts  in  this  research,  a  transfer  length  plot  can 


53 


yield  an  Ij  value  which  Is  greater  than  the  physical  contact  length, 
xc.  Clearly  In  such  a  case  the  above  assumption  of  an  Infinitely  long 
contact  Is  not  valid  and  a  more  exact  analysis  must  be  performed. 

The  more  exact  analysis  Is  based  on  the  model  given  in  Figure 
3.7.  This  model  is  very  similar  to  that  of  Figure  3.5  except  for  a 

finite  contact  length  x^  and  a  change  In  coordinates.  In  this  case, 

since  it  will  be  easier  to  perform  an  Integration  of  current  from  right 
to  left,  the  origin  of  the  coordinate  system  is  at  the  right  end  of  the 
contact. 

The  basic  assumptions  used  previously  still  hold.  The 
conducting  layer  is  described  by  a  sheet  resistivity  of  Rg.  The  layer 
is  thin  and  voltage  drops  due  to  the  vertical  flow  of  current  are 
negligible.  The  contact  itself  is  perfectly  conducting  and  is  at 
ground  potential.  The  effect  of  the  MS  interface  is  lumped  into  the 
distributed  specific  contact  resistivity  term  R^.  Current  enters  the 
contact  from  the  left  of  x  =  x^  and  begins  to  be  taken  to  ground 

through  Rj.  The  contact  width  is  equal  to  Z. 

At  any  point,  X,  the  current  lost  to  the  contact,  dl(x),  is 
given  by: 

V(X) 

dI(X)  =  -  Z  dX  (3.32) 

Rc 

Since  all  the  current  entering  from  the  left,  It0T»  i^^^t  eventually  go 
into  the  contact  between  x  =  Xc  and  x  =  0: 

Xc  V(X) 

ItOT  =  /  dl(x)  =  /  -  Z  dX  (3.33) 

''c 

and  this  gives  a  boundary  condition  on  V(x). 

If  Vg  represents  V(X  =  0) ,  a  constant  which  will  be  calculated 
later,  the  voltage  V(a)  at  any  point  0  <  a  <  x^  ,  is  the  sum  of  Vg  and 
the  total  Rjdl  voltage  drops  up  to  that  point  giving: 

a  Rj 

V{a)  =  Vg  +  /  dI{X)--(a-X)  (3.34) 

“  2 

The  [Rs/Z](o-X)  term  is  the  resistance  along  the  conducting  layers  up 
to  the  point  x  =  a  ,  through  which  I(x)  flows  creating  a  voltage  drop. 

If  Eq.  (3.34)  is  rewritten  using  Eq.  (3.32)  the  result  is: 
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(3.35) 


V(a) 


Vo 


r  «  R5(a.X)  dA 


and  this  is  the  equation  which  must  be  solved  to  determine  V(A). 

If  the  derivative  with  respect  to  a  is  taken  on  both  sides  of 
Eq.  (3.35)  the  result  is; 


d^V(a) 

da^ 


—  V(a) 


(3.36) 


which  is  the  same  as  Eq.  (3.24).  This  is  to  be  expected  since  the  only 
difference  between  the  previous  model  and  this  more  exact  model  lies  in 
the  boundary  conditions.  Therefore,  the  same  relationship  of 


Lt 


r  ^  il/2 

L  it  J 


(3.37) 


can  be  made  to  define  a  transfer  length  relating  R(-  to  Rj. 

The  solution  to  Eq,  (3.36)  with  the  condition  of  a  finite  Vq  at 
A  =  0  is: 


V(a)  =  Vq  cosh  (-^j 


(3.38) 


If  Eq.  (3.38)  is  applied  in  the  boundary  condition  of  Eq.  (3.33) 
for  the  total  current,  the  result  is: 


Vo 


I  TOT  Rc 


ZLTSinh(-r^] 


and 


V(a)  = 


Lt^ 

I  TOT  Rc 


ZLTSinh( 


(3.39) 


cosh  (^} 


(3.40) 


If,  for  comparison,  tne  coordinate  system  is  changed  back  to  the  same 
as  that  in  the  previous  development,  then  : 


and 


a  =  Xc  -  X 


V(x) 


I  TOT  Rc 


ZLTSinh(-^j 


cosh  (- 


Xr-X 


l-T 


(3.41) 

(3.42) 
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Eq.  (3.42)  is  the  more  exact  expression  for  Eq.  (3.25).  Given  that: 


and 


2 


(3.43) 


cosh  {^) 


exp  (•^)  exp  (■^)  +  exp  (^)  exp  (■^) 


2 

(3.44) 


for  the  case  of  an  infinitely  long  or  semi -infinite  contact  in  which 
Lt  «  xc  : 


»  1 


then 


Xc 

sin  h  — 
t-T 


2 


and 


cosh  (^^~^) 


exp  (•^) 

- exp(-^) 

2 


(3.45) 


(3.46) 


(3.47) 


and  if  Eqs.  (3.46)  and  (3.47)  are  substituted  into  Eq.  (3.42)  the 
result  is: 


V(x) 


or 


ItotRc 


(3.48) 


ItotRc  ,  X 

V(x)  =  -  exp  (-7--)  (3.49) 

ZLt 
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Si  nee  from  Eq.  (3.37) : 


The  transfer  length  plot  yields  a  value  of  transfer  length  when 
the  V(x)  line  is  extended  past  the  point  x  =  0,  with  a  slope  given  by 
Eq.  (3.52),  and  intersected  with  the  V{x)  =  0  axis.  The  equation  of 
this  extended  line  is  given  by: 


V(x)  =  - 

ItOT^s  , 

-  (x) 

Z 

+  V(0) 

(3.57) 

where  V(0)  depends  on  whether  Eq. 
the  case  of  an  infinite  contact. 

(3.25)  or  Eq.  (3.42)  is  used, 
from  Eq.  (3.25): 

For 

V(0)  =  - 

ItOTRsLt 

z 

(3.58) 

therefore: 

V(x)  =  - 

IT0T«s 

-  (X) 

z 

^TOT^sLt 
+  - 

z 

(3.59) 

and  when  V(x)  =  0,  x 

=  Lt. 

For  the  more  exact  case,  using  Eq.  (3.42): 


V(0) 


ItOT^c 

ZLt 


cosh  (^) 
sinh  (■^) 


and  therefore  (using  Rj  =  Rc/Lj^in  the  -IjOT^^s/^  term): 


(3.60) 


Itot^c  ItotRsLt 

V(x)  =  -  -  (x)  +  - 

ZLt2  ZLt 


and  when  V(x)  =  0, 

X  »  xlt 


cosh  (•^) 
sinh  (^) 


cosh  [^) 
sinh  (■^) 


(3.61) 


(3.62) 


where  xlj  is  referred  to  in  this  research  as  the  apparent  transfer 
length. 
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It  is  clear  from  Eq.  (3.62)  that  for  a  given  combination  of 
and  Rj,  a  transfer  length  plot  can  yield  a  value  of  xlj  greater  than 
X(;.  In  such  cases  a  calculation  of  : 

«c  =  Rs(xlt)^  (3.63) 


would  be  in  error. 


The  solution,  for  cases  in  which  the  assumption  is  not  valid,  is 
to  take  the  value  of  xlj  as  determined  from  a  transfer  length  plot  and 
using  the  measured  value  of  x^  in  Eq.  (3.62),  calculate  the  proper 
value  of  Lj  to  use  in  Eq.  (3.50). 


Figure  3.8  is  a  graph  of  Eq.  (3.62)  giving  the  relationship 
between  xjj  and  Lj  for  an  observed  value  of  xjj  with  respect  to  x^. 
The  figure  was  determined  by  first  rewriting  Eq.  (3.62)  as: 


I  cosh  ("r^) 

xlt  Lj 

xc  xc  sinh 


cosh 

l^)  sinh  [^) 


(3.64) 


Then,  for  various  values  of  Xc/Lj  .  the  right  hand  side  of  Eq.  (3.64) 
was  calculated  to  give  xlj/Xc*  The  product  of  the  chosen  Xc/Lj 
value  and  the  resultant  xlt/xc  value  gives  the  value  of  xjt/^c  value 
gives  the  value  of  xlt/Lt  Trom: 


xc  xlt  xlt 

“  •  ■  ■'  ' 

Lt  Xc  Lj 


(3.65) 


The  figure  then  represents  a  plot  of  xlj/xc  vs.  xlt/Lj. 

To  use  Figure  3,8,  first  a  measurement  of  xlj  is  made  and  a 
value  of  xlj/Xc  is  calculated.  From  the  figure,  the  corresponding 
value  of  xlj/Lj  is  found  and  Lj  is  calculated  and  used  to  find  Rq/ 

The  other  approximation  extreme  is  that  Xc/Lj  >>  1  or  even  x^/lj 
<  1.  This  would  occur  if  Rc  was  very  much  larger  than  Rs  leading  to  a 
large  Lj.  In  such  a  case,  in. the  limit  when  xq/Lj  <<  1: 


cosh  - 

1 

(3.66) 

sinh  - 

xi: 

Lt 

(3.67) 
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These  substituted  into  Eq.  (3.60)  give: 


V(0) 


ZLr 


simplifying  to: 


V(0) 

or 

Rc 


^TOTRc 

2xc 


V(0) 

-  (Zxc) 

I  TOT 


(3.68) 

(3.69) 


Eq.  (3.69)  is  a  statement  that  for  large  in  comparison  to  Rg,  the  IR 
drop  in  the  contact  is  almost  entirely  across  Rq.  In  such  a  case 
current  crowding  at  the  contact  leading  edge  is  minimal  and  most  of  the 
IjOT  current  flow  is  perpendicular  to  the  contact.  With  V(0)  being  the 
voltage  drop  across  this  contact  due  to  Ijqt  ^xc  being  the  total 
contact  area,  Eq.  (3.69)  is  the  same  as  Eq.  (3.31)  giving  RcoNT  ^ 
nonplanar  contact. 


In  Figure  3.8,  the  straight  line  represents  a  plot  of  the 
equation: 


XLT 

Xc  Lj 


(3.70) 


which  is  for  the  case  when  Xc/Lj  <<  1.  The  derivation  of  Eq.  (3.70) 
from  Eq.  (3.64)  is  made  using  the  above  assumption  and  Eqs.  (3.66)  and 
(3.67). 

From  Figure  3.8  it  is  seen  that  if  the  value  of  xlj  <  xq  , 
then  the  infinite  Xc  approximation  of  xl7=  L7  is  valid  to  about  1%. 

Also  if  X17  >  10  Xc  ,  Eq.  (3.69)  can  be  used  in  a  simple  calculation 
using  V(0)  ,  IjqT  ^nd  contact  area  to  find  R^  accurate  to  about  IX. 

Berger  (61),  using  a  transmission  line  model  similar  to  Figure 
3.7,  developed  equations  for  V(x)  and  I(x)  which  could  be  used  to  cal¬ 
culated  Eq.  (3.42)  and  the  other  subsequent  equations.  However,  the 
development  given  in  this  section  closely  parallels  the  development 
presented  on  infinite  length  contacts,  and  therefore  comparisons  to 
to  that  case  are  more  straightforward  than  if  the  method  of  Berger  was 
used. 
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For  the  case  of  Xq/Lj  <<  1,  it  is  very  easy  to  determine  from 
Eq.  (3.69)  and  from  a  nueasurement  of  the  voltage,  Vi  of  the  pick-off 
contact  strip  just  before  the  end  contact.  The  voltage  drop  across  the 
contact,  V(0),  is  then  given  by: 

ItotRs^L 

V(0)  =  Vi - ,  (3.71) 

Z 

where  AL  is  the  separation  between  the  voltage  pick-off  strip  and  the 
contact  edge. 

The  model  as  shown  in  Figure  3.5  and  Figure  3.7  strictly  depicts 
current  flow  for  a  reverse  biased  MS  contact  when  an  n-type  conducting 
layer  is  used.  Actual  contact  performance  in  terms  of  Rc  and  the  I-V 
response  were  measured  on  the  opposite,  forward  biased  contact,  of  the 
contact  structure.  However,  all  the  equations  developed  herein  can 
also  be  applied  to  the  forward  biased  contacts  since  the  models  used 
are  still  the  same.  The  only  changes  are  in  the  sign  of  the  current 
and  the  reference  for  the  contact  voltage,  being  V^pp  instead  of 
ground. 


If  an  arbitrary  reference  of  V/\pp  =  0  is  used,  then  V(x)  for  the 
forward  bias  contact  will  be  negative  and  will  represent  the  voltage 
difference  between  Vapp  and  the  voltage  in  the  conducting  layer.  The 
value  of  transfer  length  Lj  or  xlt  will  have  the  same  meaning  and  can 
be  used  in  all  the  equations  developed  previously. 

3.4  Current-voltage  (I-V)  measurement 

Another  important  experimental  measurement  used  in  this  research 
was  the  variation  of  current  with  voltage.  It  is  easily  seen  from  the 
I-V  equations  derived  in  Chapter  II  that  the  FE  or  TFE  tunneling  con¬ 
tact  exhibits  an  exponential  I-V  response.  For  a  forward  bias  greater 
than  a  few  kT/q  the  slope  of  the  I-V  curve  gives  the  value  of  E  or  Eo 
and  from  these  the  doping  level  under  the  contact  can  be  deternffned. 

The  extension  of  the  I-V  curve  to  the  V=0  axis  gives  Js>  from  which  a 
value  of  <»B  can  be  calculated.  Therefore,  the  major  parameters  used  in 
detailing  theoretical  contact  performance,  $b  and  0$,  can  be  readily 
determined  with  I-V  measurements. 

The  actual  procedure  used  to  determine  4b  ^s  method 

was  first  to  measure  I-V  data  on  the  various  contact  structures  and 
then  to  curve  fit  the  theoretical  equations  to  the  experimental  data. 

A  suitable  fit  meant  that  the  general  form  of  the  equations,  and 
therefore  the  general  aspects  of  the  Schottky  model,  were  valid  to  use 
in  explaining  contact  performance. 

A  stronger  argument  for  the  use  of  the  tunneling  model  comes 
from  relating  the  values  of  4b  ^s  »  determined  above,  to  values 
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of  (t>B  and  Cg  taken  from  other  measurements.  With  the  contact 
structures  used  in  this  research,  the  other  commonly  used  methods  of 
measuring  <{>0  were  not  applicable  and  a  further  discussion  of  this  fact 
will  be  given  later.  However,  an  independent  measurement  of  Cs,  using 
the  van  der  Pauw  technique,  did  yield  a  value  which  could  be  compared 
to  that  determined  from  an  I-V  measurements.  If  the  values  of  Cg  as 
determined  by  the  two  different  techniques  are  the  same,  then  the 
contacts  are  presumed  to  operate  as  tunneling  contacts.  A  further 
verification  of  tunneling  operation  was  made  when  the  41B  and  Cg  values 
from  the  I-V  fit  were  used  in  the  theoretical  equations  for  and 
compared  to  the  actual  variation  of  vs.  temperature. 

One  problem  in  relating  the  theoretical  I-V  equations  to 
experimental  data  is  that  the  current  flow  in  the  actual  contact 
structures  is  planar,  giving  rise  to  current  squeezing  at  the  edge  of 
the  contact.  The  effect  is  analogous  to  the  base  current  crowding  in  a 
bipolar  transistor  and  can  be  analyzed  using  the  same  equations  (62). 

However,  for  this  research,  a  development  of  the  effective  I-V 
response  for  the  contact  was  performed  using  a  model  similar  to  those 
given  in  the  development  of  R^.  The  only  difference  is  that  in  the 
model  instead  of  a  distributed  Rc,  distributed  diodes  with  a  response: 

J  =  Jg  exp  (3.72) 

were  used.  The  basic  form  of  Eq.  (3.72)  is  valid  for  both  FE  and  TFE 
operation,  as  can  be  seen  by  comparison  to  Eq.  (2.40)  and  Eq.  (2.58) 
(for  biases  greater  than  a  few  kT/q).  For  the  FE  case: 

E'  =  E  (3.73) 

00 

and  for  the  TFE  case: 

E'  =  Eo  (3.74) 

There  is  a  slight  dependence  in  Jg  on  V  but  since  the  effect  of 
V  through  the  exponential  term  is  the  dominant  factor,  the  dependence 
in  the  Og  term  can  be  neglected  with  only  sight  error. 

Figure  3.9  shows  the  model  of  the  contact  as  changed  for  use  in 
the  development  of  the  I-V  equations.  The  figure  shows  the  case  of  a 
forward  bias  contact  made  on  an  n-type  conducting  layer  with  sheet 
resistivity  Rg.  As  before,  voltage  drops  due  to  vertical  current  flow 
in  the  conducting  layer  are  neglected.  For  convenience,  the  voltage  on 
the  metallic  contact  of  length  Kq  is  taken  as  ground.  Therefore,  the 
values  derived  for  voltage  in  the  conducting  layer,  V(x),  will  be  nega¬ 
tive. 

The  total  current  flowing  through  the  contact  is  IjoT  enters 
the  conducting  layer  a  point  x  =  x^.  Since  the  voltage,  V(x)  ,  for 
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I 


I 


Figure 


I  1 

XsQ  XSX, 

V(0) 


3.9  Model  used  to  derive  I-V  equations  for  the 
planar  contact  structure. 
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0  <  X  <  Xc  cannot  actually  be  measured,  experimental  data  consists  of 
a  measurement  of  IjoT  vs.  V(x  =  x^)/  Actually,  V(xc)  also  cannot  be 
directly  measured  but  can  be  easily  determined  from  the  measured 
voltage,  V^,  of  the  first  voltage  pick-off  strip  to  the  right  of  the 
contact.  Namely: 

IjotML 

V(xc)  =  Vi  +  -  (3.75) 

Z 


Eq.  (3.75)  is  exactly  the  same  as  Eq.  (3.71)  taking  into  account 
that  the  direction  of  current  flow  is  reversed.  The  aL  is  the  differ¬ 
ence  in  distance  between  the  x  =  x^  point  and  the  location  of  the 
voltage  pick-off  strip. 


The  incremental  current  flow  through  the  contact  at  any  point  x 
is  given  by: 


I(x)  =  J(x)Zdx  (3.76) 

where  J(x)  is  given  by  Eq.  (3.72).  The  voltage  at  any  point  0  <  x  <  Xc 
is  the  difference  between  V(x  =  0)  and  the  total  IRj  drops  up  to  point 
X.  Therefore,  analogously  to  Eq.  (3.34): 


V(x)  =  V(0)  -  /  J(y)Z  ^(x-y)dy 


(3.77) 


Since  V(0)  is  a  constant,  taking  the  derivative  of  both  sides  gives: 


dWxl 

ox 


-  /  J(y)Rs^y 

o 


(3.78) 


Taking  the  derivative  again  gives: 

=  -d(x)Rs  (3.79) 

and  using  Eq.  (3.72)  gives: 

=  .Rjjjexp  (-  ^^)  (3.80) 

and  this  is  the  equation  which  must  be  solved  for  V(X).  (If  the 
assumption  is  made,  as  used  in  the  calculation  of  Lj,  that  J(x)  is 
given  by  a  linear  V(x)/Rc  term,  then  Eq.  (3.36)  results.  ) 

The  boundary  conditions  are  determined  from  the  fact  that  all 
incremental  dl(x)'s  must  add  up  to  lyoT  x  =  xq.  Therefore: 
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/  J(y)Zdy  *  I70T 


(3.81) 


and  using  Eq.  (3.78)  gives  the  conditions: 

=  0 


dV(x) 

dx 


X  =  0 


and 


dx 


X  =  X, 


>^slTnT 

Z 


The  solution  for  V(x)  is  easily  calculated  assuming  a  form: 

V(x)  =  A  In[u(x)] 

with  A  a  constant  and  u(x)  of  the  form: 
u(x)  =  K  cos  y  X 

With  the  above  assumptions,  the  solution  is  given  by  : 


(3.82) 


(3.83) 


(3.84) 

(3.85) 


V(x) 


2E'  qRsJs  cos  y  x 

- in  ( - )  - 

q  2E'  Y  X 


(3.86) 


The  parameter,  y  ,  depends  on  the  boundary  conditions  as  given  in  Eq. 
(3.83)  and  is  determined  from: 


dlTOT^s 

Y  tan  (yxc)  =  - 

2E'Z 


(3.87) 


where  Z  =  contact  width 
Xc  =  contact  length 
q  =  electronic  charge 

As  a  check  on  Eqs.  (3.86)  and  (3.87)  the  limiting  case  of  Rs  fl 
can  be  examined.  For  R-  0  ,  the  IRj  voltage  drops  go  to  zero  and  the 
current  crowding  at  the  edge  is  reduced.  In  this  case  the  dominant 
voltage  drop  is  V(x)  across  the  contact  itself  and  lateral  current  flow 
does  not  have  a  large  effect.  The  contact  itself  and  lateral  current 
flow  does  not  have  a  large  effect.  The  contact  I-V  should  tend  to  that 
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for  a  contact  with  perpendicular  current  flow. 


From  Eq.  (3.87),  when  Rg  +  0  then  y  *  0  and  therefore  yX(-  -*•  0. 

Therefore; 


sin(YXc)  -  YXc 
and 


cos(yxc)  +  1 
and 


tan(YXc)  YXc 

substituting  in  Eq.  (3.87)  gives: 
,  qlJOTRs 

Y^Xc  =  - 

2E'Z 


or 


Y 


qlTOTf^s 

2E'XcZ 


(3.88) 

(3.89) 

(3.90) 

(3.91) 

(3.92) 


The  voltage  called  the  "contact  voltage"  in  an  experimental  measurement 
is  V(xc).  Therefore,  substituting  Eqs.  (3.92)  and  (3.89)  in  Eq.  (3.86) 
for  V(xc)  gives: 


V(xc) 


2E'  qRgJs  2E'Zxc 

—  In  (— ]  ( - ) 

q  2E'  glTOT^^c 


(3.93) 


or  simplified  and  rearranged: 


I  TOT  =  Js(Zxc)  exp 


qV(xc) 

E' 


(3.94) 


since  Zx^  is  the  contact  area  then,  as  expected,  Eq.  (3.94)  is  the  I-V 
equation  for  a  MS  diode  with  perpendicular  current  flow  for  biases 
greater  than  a  few  kT/q. 
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In  using  Eq.  (3.86)  and  (3.87)  in  curve  fitting,  the  terms  of 
Rj,  Z,  and  Xj.  were  given  data.  Terms  Z  and  x,;  were  determined  from 
dimensional  measurements  and  Rj  was  determined  from  the  van  der  Pauw 
measurement. 

For  an  assumed  value  of  E'  and  Jj,  y  was  calculated  using  Eq. 
(3.87)  for  a  specified  IjoT-  Then,  using  y  and  Eq.  (3.86),  the 
corresponding  contact  voltage,  V(xc)  ,  for  that  value  of  IjoT  cal¬ 
culated.  These  steps  were  performed  for  a  range  of  Ijqt  values  and  the 
resulting  Ijqt  I'ssponse  was  compared  to  the  experimental  I-V 

data.  The  initial  assumptions  of  E*  and  0$  were  then  changed  and  the 
procedure  repeated  until  a  good  fit  between  theoretical  and  experi¬ 
mental  data  was  achieved. 

In  using  the  values  of  E'  and  to  find  ((>b  and  Cj,  the  type  of 
emission,  FE  or  TFE,  needed  to  be  determined.  The  value  of  E'  was 
first  set  equal  to  Eo  and  Eq.  (2.56)  was  used  to  find  £»  .  This  meant 
that  the  emission  was  initially  assumed  to  be  TFE.  Using  the  resulting 
value  of  E«,  in  Eq.  (2.11)  the  doping  level,  Cj  ,  (actually  Np)  was 
calcul ated. 

Using  the  calculated  value  of  Cj,  a  value  for  the  Fermi  level, 
np.  was  calculated.  Using  these  values  of  E»  and  np  and  using  Eq. 
(2.45),  a  value  of  fj  was  calculated.  Since  Eq.  (2.25)  : 


—  -  Cl  >  (2fi) 

kT 


(2.25) 


must  be  satisfied  in  order  for  the  emission  to  be  FE,  the  calculated 
value  of  fi,  was  used  to  find  Cimax  Tor  the  given  T  from: 


ClMAX 


1 

kT 


-  (2f:) 


(3.95) 


This  value  of  Ci(y|/\x  was  used  in  Eq.  (2.44): 


In(M] 

^np  ■' 

Cl  =  -  (2.44) 

2E« 

to  determine  a  value  of  4bmAX*  Since  the  value  of  Js  is  inversely  pro¬ 
portional  to  (t»B,  use  of  4BMAX  Cq.  (2.41)  to  determine  under  FE, 
results  in  the  minimum  value  of  Js,  Jsmtn*  would  be  possible  and 

still  have  the  emission  be  FE. 
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Using  the  value  of  E'  obtained  from  the  best  fit,  a  value  of  E* 
was  calculated  and  determined.  If  the  value  of  from  the  curve 

fit  is  less  than  JsMIN*  a  $5  >  (fBMAX  ’s  needed.  However, 

vB  ^  'fBMAX  results  in  c>  ci,,..  which  violates  Eq.  (2.25)  for  the  given 
Eco  .  Therefore,  if  J5  from  tne  curve  fit  is  less  than  Jsmii\|  deter¬ 
mined  above,  the  emission  cannot  be  FE. 

If  J5  from  the  curve  fit  is  greater  than  ^  4’B  *^0^5 

exist  in  the  range  that  allows  Eq.  (2.25)  to  be  satisfied  and  still  fit 
J5  and  E„  to  the  FE  equation.  In  such  a  case,  the  value  of  Js  is  used 
in  Eq.  (2.41)  to  determine  a  value  of  (tis  •  This  4)b  and  the  value  of  C5 
calculated  before  then  represent  the  Schottky  barrier  conditions 
necessary  to  result  in  the  experimental  I-V  data  for  that  particular 
sampl e. 


The  initial  assumption  of  E  ‘  =  Eo  instead  of  E  '  =  £00  causes  no 
problem  in  the  above  procedure,  even  though  for  the  case  where 
'-’s  ^  '^smIN  ’  emission  is  FE,  the  E'  actually  represents  Eoo  and 

not  Eo.  This  is  because,  in  the  range  where  Eq.  (2.25)  is  satisfied 
and  FE  is  indicated,  the  appropriate  E®,  values,  when  put  into  Eq. 

(2.56)  result  in  Eo  =>  E«.  Therefore,  an  initial  assumption  of  E‘  =  Eo, 
instead  of  E '  =  ,  causes  no  problem  even  if  FE  is  indicated  in 

further  calculation. 

From  some  actual  experimental  data  at  room  temperature,  if  the 
above  calculation  indicates  that  the  emission  is  FE,  then  the  value  of 
Eo  from  Eq.  (2.56)  is  equal  to  £«  well  within  1%.  Therefore,  the 
initial  assumption  of  E '  =  Eo  instead  of  E '  =  E^  is  only  off  by  a  very 
small  factor. 

For  the  case  where  the  above  comparison  of  and  the  calculated 
JsMi|\|  results  in  J5  <  Js|v|rw  .  the  emission  is  not  FE  and  Eq.  (2.41) 
was  not  used.  This  meant^that  the  emission  may  be  TFE.  Using  the  same 
values  of  E  '  =  Eo,  and  calculated  values  of  £«,  ,  np  .  and  Np,  Eq. 
(2.59)  with  V  =  U  was  used  to  calculate  a  to  match  the  measured 
to  the  Js  expected  under  TFE. 

To  check  whether  the  calculation  for  TFE  was  strictly  valid,  the 
calculated  values  of  E„,  np  .  and  4)5  were  used  in  Eq.  (2.61).  If  Eq. 
(2.61)  was  satisfied  then  the  use  of  Eq.  (2.59)  was  strictly  valid.  If 
Eq.  (2.61)  was  not  satisfied  then  there  is  a  combination  of  FE  and  TFE 
occuring  and  the  value  of  Js  as  calculated  using  Eq.  (2.59)  is  slightly 
in  error. 

3.5  Other  measurements 


In  addition  to  an  I-V  measurement,  there  are  two  other  types  of 
measurements  which  are  commonly  made  to  find  4>b,  a  C-V  measurement, 
and  a  photoelectric  measurement  (9). 
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The  C-V  technique  measures  the  variation  in  space  charge  region 
capacitance  of  a  Schottky  barrier  with  reverse  bias  voltage.  A  plot 
of  (l/c2)  vs.  Vr  gives  a  line  whose  slope  is  related  to  doping,  Cj,  and 
whose  extrapolated  intercept  on  the  (1/c^)  =  0  axis  gives  ((10. 

For  this  research  the  C-V  technique  was  attempted  but  was 
unsuccessful.  The  source  of  difficulty  was  the  large  tunneling 
current.  In  a  normal  diode  structure  the  only  current  flow  under 
reverse  bias  is  due  to  leakage  and  saturation  currents  and  is  very 
small.  In  this  research,  because  C5  is  large  in  order  to  realize 
tunneling  and  give  a  low  R^,  a  large  current  flows  under  reverse  bias 
and  the  presence  of  this  current  appears  as  a  large  leakage  resistance. 
Most  instruments  cannot  measure  the  variation  in  C  when  a  large 
effective  leakage  resistance  is  present  and  such  was  the  case  here. 

The  photoelectric  technique  is  used  to  measure  the  variation  in 
photo-generated  current  versus  photon  energy  for  a  Schottky  barrier 
illuminated  by  a  monochromatic  light  source.  The  mechanism  involved  is 
the  photoelectric  excitation  of  an  electron,  at  the  Fermi  level  in  the 
metal,  up  over  the  Schottky  barrier  thereby  contributing  to  current 
f  1  ow. 


The  problem  encountered  when  such  measurements  were  tried  on  the 
contact  structures  used  in  this  research  was  that  no  detectable  photo¬ 
current  was  observed.  The  reason  is  that  the  relatively  thick  metal 
layer  does  not  allow  the  penetration  of  incident  photons  down  to  the 
region  next  to  the  MS  interface  where  they  could  excite  electrons. 
Electrons  excited  away  from  MS  interface  recombine  before  they  can 
reach  the  junction  and  therefore  do  not  give  rise  to  any ’measurable 
photo  current. 

In  addition  to  the  problem  of  the  thickness  of  the  Au,  the 
relationship  between  photon  energy  and  4>b  in  the  presence  of  electron 
tunneling  through  a  thin  barrier  is  not  well  known.  In  the  theory 
usually  used  with  the  photoelectric  technique  it  is  assumed  that  the 
photocurrent  is  due  only  to  electrons  emitted  over  the  barrier.  How¬ 
ever,  for  this  research,  the  contacts  exhibited  appreciable  current 
flow  due  to  tunneling.  Electrons  excited  thermally,  as  in  TFE,  or 
photoelectrical ly  would  be  expected  to  be  able  to  tunnel  in  a  similar 
manner  and  therefore  photoelectric  current  data  for  contacts  in  which 
tunneling  occurs  may  yield  erroneous  values  of  41 5  if  the  simple 
emission  theory  is  used. 
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SECTION  IV 


EXPERIMENTAL  DATA 


4.1  Diffusions 


The  GaAs  starting  substrates  were  obtained  from  Crystal 
Specialities  (64)  and  the  manufacturer's  specifications  were: 


p-type 

Carrier  concentration 
Mobi 1 1 ty 
Orientation 
Resistivi ty 
Wafer  thickness 
One  side  polished 


Zn  doped 
5  X  lO^S/cm^ 
378.6  cm  /  V-sec 
(100) 

0.37  n  cm 
18  mils 


Small  sections  of  each  wafer  were  scribed  out  and  cleaned  with  a 
trichloroethylene,  acetone,  methanol,  and  H20  wash  sequence.  Ohmic 
contacts  were  made  to  the  periphery  of  the  polished  side  by  alloying 
small  dots  of  In--  2%  Zn  alloy  in  an  annealing  furnance  at  225°  for  10 
min.  The  I-V  curve  of  the  sample  and  contacts  was  observed  on  a  curve 
tracer  to  check  for  ohmicity. 

Each  main  wafer  was  then  retested  to  determine  average 
resistivity,  p  ,  carrier  concentration,  Cg,  and  mobility,  pg,  with  the 
van  der  Pauw,  Hall  effect  technique.  Table  4.1  lists  the  results  of 
this  measurement.  The  entries  in  Table  4.1  differ  from  the  manufactur¬ 
er's  data  but  these  measured  values  were  taken  to  be  the  more  valid 
characterization  of  the  starting  wafers. 

Before  each  diffusion,  the  wafers  to  be  diffused  were  given  a 
standard  wash  and  etch  preparation  as  listed  in  Table  4.2.  The  MB  etch 
was  reported  to  result  in  the  least  number  of  surface  defects  (21). 

After  cleaning,  the  Sn  dopant  solution  was  spun  on  the  wafer 
surface.  As  explained  in  Section  II,  the  spin-on  dopant  solution  was  a 
mixture  of  Sn  bearing  Tinsilicafilm  and  undoped  Silicafilm.  The  pro¬ 
portions  of  the  two  solutions  were  varied  to  give  different  values  of 
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dopant  surface  concentration.  After  spin  on,  the  film  was  densified 
into  a  Sn  bearing  Si2  layer  by  hearing  the  wafer  for  20  min  to  1  hour 
at  200°C  in  a  furnace.  A  normal  air  atmosphere  was  present. 

The  dopant  coated  wafers  were  placed  in  a  closed  but  not  tightly 
sealed  capsule  along  with  extra  As  and  GaAs  pieces.  The  diffusion  was 
at  950°C  for  6  hours.  An  Ar  flow  was  maintained  in  the  diffusion 
furnace.  The  arsenic  released  during  the  diffusion  was  kept  out  of  the 
room  by  means  of  a  scavenger  box  located  at  the  mouth  of  the  furnace. 
Figure  4.1  is  a  photograph  of  the  diffusion  furnace  set  up. 
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Table  4.1 


GaAs  Starting  Substrate  Specifications 


Sample  # 

p  n-cm  X  10"^ 

PH  cm^/V-sec  X  10^ 

Cft/cm^  X  10^ 

3A-1 

1.65 

1.937 

1.95 

4A-1 

2.40 

2.367 

1.10 

5A-1 

1.46 

2.135 

2.00 

6A-1 

1.44 

2.220 

1.96 

7A-1 

1,46 

2.155 

1.98 

8A-1 

2.52 

2.417 

1.03 
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Table  4.2 


Wafer  Preparation  Process 

1.  1  min.  wash  in  trichloroethylene  in  ultrasonic  cleaner. 

2.  1  min.  wash  in  acetone  in  ultrasonic  cleaner. 

3.  1  min.  wash  in  methanol  in  ultrasonic  cleaner. 

4.  Rinse  and  soak  in  high  purity  (18  Mfi-cm)  H^O. 

5.  MB  etch  for  2  min. 

1  HF  :  1  HCL  :  4  H^O  +  1  drop  H^O^/  10  ml  of  solution 

6.  Rinse  in  18  M  fl-cm  H^O. 

7.  Blow  dry  with  prepurified  N^. 

8.  Store  in  petri  dish. 
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Table  4.3 


Diffused  Layer  Data  Determined  from  van  der  Pauw  Measurements 


Di f fused 
Sample  # 

microns 

p  n-cm  X  10“^ 

WH  cm^/V-sec  x  10^ 

Cr/cih^  X  10^^ 

4A-1 

1.64 

1.35 

1.561 

2.96 

5A-1 

1.37 

1.31 

1.431 

3.33 

6A-1 

1.47 

1.54 

1.350 

3.00 

7A-1 

1.02 

2.46 

1.365 

1.87 

8A-1 

1.64 

3.99 

1.412 

1.11 

5-2 

1.58 

1.62 

1.611 

2.40 

5-3 

1.83 

2.06 

1.609 

1.89 

5-4 

1.73 

2.13 

1.675 

1.75 

6-1 

1.91 

1.16 

1.678 

3.21 

6-2 

1.68 

1.77 

1.631 

2.16 

6-3 

1.67 

2.33 

1.709 

1.57 

6-4 

1.44 

1.78 

1.482 

2.37 

6-5 

1.09 

2.09 

1.937 

1.54 
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4.1  Pnotograph  of  the  diffusion  furnace  equipment. 


•ft 


m 


Pieces  of  the  diffused  wafers  were  hand  lapped  with  0.3  p 
alumina  powder  on  a  glass  plate  and  the  junction  was  delineated  using  a 
commercial  staining  solution  (48)  and  exposure  to  light  from  a  micro¬ 
scope  lamp.  Junction  depths,  xj,  were  then  measured  with  a  Watson 
i nterferometer. 

The  n-type  layers  were  then  characterized  with  van  der  Pauw  Hall 
effect  measurements,  and  values  of  layer  average  resistivity,  p  , 
average  Hall  mobility,  uh  ,  and  average  carrier  concentration  Rp  ,  were 
determined.  The  results  of  the  junction  depth  measurements  and  van  der 
Pauw  measurements  are  given  in  Table  4.3. 

Using  the  curve  of  p  vs.  0$  for  full  ionization.  Figure  3.3,  and 
the  values  of  p  from  the  van  der  Pauw  measurements,  the  value  of  Cj  for 
each  sample  can  be  determined.  These  values  are  listed  in  Table  4.4. 

Cs  can  also  be  determined  using  the  equation: 

X 

C(x)  =  Cjerfc  -  (4.1) 

/Dt 

if  the  diffusion  coefficient,  D,  for  Sn  in  GaAs  at  the  diffusion 
temperature  is  known.  Assuming  full  ionization,  when  x  =  xj,  then  C(x) 
=  Cg.  Therefore: 


Cg  *  Cj  erfc  - 

2/W 


(4.2) 


and  Cs  cam  be  determined. 

The  value  of  D  for  Sn  in  GaAs  has  been  investigated  by  various 
researchers  and  is  usually  given  in  the  form: 

Eo 

D  =  Do  exp [ - )  (4.3) 

kT 

and  Fane  and  Goss  (66)  gave  data  from  which  values  of 

Do  *  4.957  X  10-2  (4.4) 

Eo  =  2.725  eV  (4.5) 

can  be  calculated.  These  are  also  reported  in  Kendall  (67)  and  Casey 
(68).  Kendall,  gives  a  plot  of  Eq.  (4.3)  based  on  Fane  and  Goss's 
data,  but  Casey,  stating  that  his  listed  values  are  also  from  Fane  and 
Goss's  work,  nevertheless  gives  different  values  of  Do  and  Eo.  Casey 
gives  values  of: 
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Do  = 
Eo  = 

Goldstein  and 
Do  = 
Eo  = 


3.8  X  10-2 
2.7  eV 

Keller  (69),  (also  reported  In  Madelung  (58)  ), 
6  X  10-**  cm2 /sec 
2.5  eV 


(4.6) 

(4.7) 
give: 

(4.8) 

(4.9) 


valid  over  their  experimental  temperature  range  of  1069“  C  to  1215®  C, 
but  possibly  Eq.  (4.3)  can  be  extrapolated  with  their  values  down  to 
the  temperature  range  less  than  1069®  C. 


Tuck  and  Badawi  (44)  report  that  the  diffusion  coefficient 
depends  on  the  nature  of  the  substrate,  giving  slightly  different 
values  if  the  substrate  is  undoped  or  previously  doped  n-type.  They  do 
not  report  data  for  a  p-type  substrate  but  their  data  for  the  undoped 
substrate  case  can  be  used  as  a  comparison.  Their  values  are: 


Do 

= 

3.5  cm2 /sec 

(4.10) 

Eo 

= 

3.3  eV 

(4.11) 

Using  the 
950®  C  can 

above  data,  the  value  of  D  for 
be  calculated,  resulting  in: 

a  diffusion  temperature  of 

0l 

s 

2.99  X  10-^3  cm2 /sec 

with  Fane  and  Goss 

(4.12) 

O2 

s 

3.056  X  10-1**  cm2/sec 

with  Goldstein 
and  Keller 

(4.13) 

03 

= 

9.066  X  10-1**  cm2 /sec 

with  Tuck  and  Badawi 

(4.14) 

The  research  of  Fane  and  Goss  actually  showed  non-erfc  profiles  and 
their  data  was  determined  by  a  curve  fit  to  a  portion  of  their  observed 
profiles.  Goldstein  and  Keller  observed  Gaussian  type  profiles  because 
of  their  limited  source  of  Sn.  Tuck  and  Badawi  observed  good  erfc 
profile  fits  over  a  large  portion  of  the  Sn  diffusion  profile. 

Using  the  values  of  Di,  Dg,  D3  in  Eq.  (4.3),  along  with  t  =  6 
hours,  a  value  of  Cs  can  be  determined  for  a  given  Cg  and  Xj.  The 
results  of  such  calculations  are  presented  in  Table  4.4.  The  results 
will  be  discussed  in  Section  V. 

4.2  Contact  fabrication 


After  the  diffusion  had  been  characterized  via  the  above,  the 
wafers  were  prepared  for  Au  contact  fabrication.  The  wafers  were  given 
the  same  wash  sequence  as  outlined  in  Table  4.2,  except  for  the  etch 
step.  Each  wafer  was  then  coated  with  photoresist  and  exposed  to  the 
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Values  of  Surface  Concentration  of  the  Diffused  Layers 
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Figure  4*2  Contact  mask. 


82 


Figure  4.. 3  Photograph  of  the  test  apparatus  used  to 

determine  R^,  from  a  measurement  of  transfer 
length. 
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Figure 


a,  STRIP  SPACING 


4.4  Example  transfer  length  plot.  Sample  #  6- IB 
T  =  18.8®C.  I  *  50  pA.  (a  =  3.66  x  10"^  cm 
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contact  mask.  The  structure  of  the  mask  is  shown  in  Figure  4.2.  After 
the  pattern  was  exposed  and  developed  the  photoresist  layer  had 
openings  where  the  Au  contacts  would  be  made. 

The  patterned  wafers  were  given  a  coating  of  approximately  1500A 
of  evaporated  Au.  The  evaporation  used  6"  of  20  mil  diameter  Au  wire 
placed  in  an  evaporation  boat  12"  away  from  the  substrates.  The 
electrical  power  applied  to  the  evaporation  boat  was  gradually 
increased  until  evaporation  began  and  was  waintained  at  a  level  which 
resulted  in  all  the  gold  in  the  boat  being  evaporated  in  two  or  three 
mi nutes. 

After  the  evaporation,  the  Au  coated  wafers  were  soaked  in 
acetone  for  up  to  one  hour  or  until  lift  off  of  the  unwanted  Au 
occured.  If  a  clean  lift  off  was  not  obtained  after  the  one  hour,  the 
beaker  was  given  a  short  dip  on  an  ultrasonic  cleaner  to  force  the 
unwanted  Au  to  lift  off.  In  most  cases  this  technique  gave  acceptable 
resul ts. 

After  the  Au  patterning,  each  wafer  was  given  a  low  temperature, 
“  250°  C,  anneal  in  N2  for  5  minutes  to  promote  better  contact  adhesion 
and  to  try  and  eliminate  any  possible  interface  nonuniformities.  The 
individual  contact  structures  were  then  scribed  apart,  ready  for 
testing  with  the  transfer  length  measurement. 


4.3  Measurement  of  contact  resistivity 

For  the  transfer  length  measurements,  a  constant  DC  current  was 
supplied  through  the  sample  using  a  Keithly  225  current  source.  The 
current  level  was  chosen  so  that  the  voltage  drop  across  the  contact 
would  be  on  the  order  of  a  few  millivolts.  This  would  give  a  good 
approximation  to  the  V  =  0  condition  used  in  defining  R^.  The  voltages 
were  measured  with  respect  to  the  contact  pad  where  the  current  exited 
using  a  Keithly  191  OVM.  Contact  to  the  pads  and  pick-off  strips  was 
made  with  tungsten  tipped  micromanipulator  probes.  Figure  4.3  is  a 
photograph  of  the  test  set  up.  The  probe  stage  rested  on  a  thermo¬ 
electric  cold  plate  which  was  used  for  subsequent  measurements  of  Rc 
vs.  T. 


As  an  example.  Figure  4.4  is  the  transfer  length  plot  for  sample 
#6-lB.  The  plot  is  for  a  temperature  of  18.8°  C,  and  the  contact 
structure  had  a  width,  Zs  =  0.0937  cm. 

The  pick  off  strip  spacing,  a,  is  equal  to  3.66  x  lO"^  cm  and 
the  contact  pad  width  and  length  is  «•  7  x  10"^  cm.  Only  voltages  from 
every  other  pick  off  strip  were  taken.  The  current  was  I  =  50  uA. 

The  straight  line  drawn  on  Figure  4.4  is  a  least  squares  fit  to 
the  data.  The  equation  for  the  line  is: 
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V{a)  =  -0.154  a  +  1.834  mV 

given  in  units  of  a.  The  total  applied  voltage  was: 

Vapp  =  2.132  mV  (4.16) 

Therefore,  the  apparent  transfer  length,  xlj  ,  for  the  forward  bias 
contact  is  given  as  ia(  when  V(a)  »  Vapp*  Solving  for  a  gives: 


xlT  =  1“1 


2.132  mV  -  1.834  mV 
-0.154  mV/a 


XLT  “  1.94  a  (4. 17) 

and  for  the  reversed  bias  contact  (voltage  reference  contact),  xlj  is 
found  from  a  when  V(a)  =  0.  Or: 

1.834  mV 

Xlj  =  -  -  10  a 

0.154  mV/a 

XLJ  =  1.92  a  (4. 18) 

It  is  seen  that  the  values  of  xlj  are  nearly  the  same  in  both 
cases.  Actually,  Rc  for  both  forward  and  reverse  biased  contacts 
should  become  equal  in  the  limit  as  V  0,  since  the  I-V  curve  is 
continuous  at  V  =  0.  If  xlt  values  for  the  two  contacts  are  very 
different,  then  contact  non -uniformity  due  to  processing  variations  is 
probably  the  cause. 

The  value  of  Rx  can  be  calculated  from: 
dV  Zs 

Rg  =  . -  (3.18) 

dx  I 

For  sample  #6-lB  from  Eq.  (4.15): 
dV 

—  s  -0.1538  mV/a  =  -4.20  x  lO’^  V/cm 
dx 

since  a  =  3.66  x  lO"^  cm.  Therefore,  with  Zs  =  0.0937  cm  and 
I  =  50  uA: 
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In  most  cases,  the  transfer  length  plot  yielded  a  very  linear 
V(x)  variation.  Therefore,  a  much  simpler  method  to  find  xlj  was  used. 
For  any  particular  contact  structure,  the  width,  Z5,  was  measured  and 
Rc  was  calculated  from  the  van  der  Pauw  value  of  p  and  x,-.  Then  Eq. 
(3.18)  :  ^ 


dV  IRg 

dx  Zs 


(4.25) 


When  Eq.  (4.25)  is  given  in  units  of  a  ,  then  4V  =  -dV/do  ‘lOa 
where  aV  is  the  IR  drop  along  the  entire  conducting  layer.  The  voltage 
drop  across  the  contacts  due  to  was  assumed  to  be  equal  and 
therefore: 


Vapp-AV 

VCONT  =  - 

2 


The  value  of  xlj  is  then  calculated  from: 


(4.26) 


VCONT 

xLT  =  -  (4.27) 

dV 

■ 


In  most  of  the  cases  in  this  research,  the  level  of  Cg  was  such 
that  XL7  »  Xq  ,  implying  large  values  of  Rf  Since  the  asymtotic 
approximation  from  Figure  3.8  is  valid,  Rc  can  be  calculated  from: 

VCONT 

Rc  =  - XcZ  (4.28) 

I 

where  x^  is  the  contact  length  and  Z  is  the  contact  width. 

Experimental  values  of  Rc  measured  at  T  =  23°  C  for  various 
contact  samples  are  shown  in  Figure  4.5.  They  are  plotted  against  the 
value  of  Cg  for  that  particular  sample  as  determined  from  the  van  der 
Pauw  measurement  of  p  (Table  4.4). 

Immediately  apparent  is  the  fact  that  not  all  the  points  lie  on 
a  single  (>b  line  as  would  be  expected  if  4»g  were  governed  by  Fermi 
level  pinning.  Since  this  result  suggested  that  ())b  varied  with 
processing,  I-V  curves  for  various  contacts  were  measured  to  determine 
independently  whether  or  not  415  was  varying. 
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4.4  Current -voltage  measurements 

The  measurement  of  the  I-V  response  of  the  contacts  was  perform¬ 
ed  by  the  method  explained  in  Section  III.  The  DC  current  was  supplied 
with  a  Keithly  225  DC  current  source  and  the  voltage  drops  were 
measured  with  a  Keithly  191  OVM,  The  voltage  drop  across  the  contact 
was  taken  to  be  the  value  of  V(X(.)  and  was  determined  using  Eq.  (3.75). 
This  voltage  was  measured  for  a  range  of  Ijqt  resulting  data 

gave  the  I-V  response  of  the  contact. 

Figure  4.6  and  Figure  4.7  give  the  I-V  response  of  various 
contact  structures.  As  expected,  the  I-V  curves  have  an  exponential 
response  although  the  slope  decreases  at  higher  current  because  of  the 
effect  of  Rs  of  the  conducting  layer.  The  slopes  are  not  all  the 
same,  indicating  differing  values  of  Cs*  Different  values  of  also 
appear  to  be  present  since  the  effective  values  of  Jg  are  not  the  same. 

Figures  4.8  through  4.13  show  the  I-V  response  of  some  of  the 
samples  in  Figures  4.6  and  4.7  along  with  curves  fitted  using  the 
theoretical  equations  given  in  Chapter  III.  The  central  position  of 
the  I-V  curves  was  taken  to  give  the  best  fit.  At  the  higher  current 
levels  the  fit  is  poorer  but  at  these  points  the  experimental  data  are 
probably  less  exact  because  of  sample  heating. 

Table  4.5  gives  the  value  of  Jj  and  £’  determined  from  the  curve 
fit.  These  values,  along  with  data  for  Rs  ,  Xq  ,  and  Z  for  any 
particular  contact,  were  used  in  Eqs.  (3.86)  and  (3.87)  to  give  the 
theoretical  I-V  response.  Also  given  are  the  values  of  Cs  and  ({>g  as 
calculated  using  E*  and  Js. 

The  label  EST  in  the  emission  column  means  that  for  this 
particular  combination  of  Cs  and  (j>g  ,  the  emission  was  out  of  the 
strict  TFE  range  and  moving  into  the  FE  range.  These  values  of  Cs  and 
(^g  are  calculated  using  the  FE  equation  as  an  approximation  and  are 
therefore  estimates.  However,  the  values  for  Cs  and  (|)g  are  calculated 
using  the  FE  equation  as  an  approximation  and  are  therefore  estimates. 
However,  the  values  for  Cs  should  be  very  close  to  actual  values  with 
the  values  of  (jig  being  slightly  more  in  error. 

The  variation  in  4g  ,  the  reason  for  the  observed  ())g  being 
breater  than  0.9  eV,  and  comparison  of  the  Cs  values  in  Table  4.5  with 
those  of  Table  4.4  will  be  discussed  in  Section  V. 

4.5  Temperature  variation  of  Rr 

Using  the  probe  station  mounted  on  the  thermoelectric  cooler,  as 
shown  in  Figure  4.3,  transfer  length  measurements  of  R^  versus 
temperature  were  made.  The  experimental  temperature  range  was  from 
-  -10“  C  to  -  +25“  C.  The  temperature  was  measured  using  a 
copper-constantan  thermocouple  counted  in  the  prove  support  plate 
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Figure  4.11  I-V  response  of  Au-Sn  diffused  sample  #  7A=1B. 

T  =  23°C.  Theoretical  TFE  fit  with  0g  =  1.324  eV 
Cc  =  9.64  X  10^®/cm^. 


^CONTACT  VOLTS 


Figure  4.12  I-V  response  of  Au-Sn  diffused  sample  #  8A-1A. 

T  =  Z3°C.  Theoretical  FE  fit  with  0r  =  1.417  eV. 

19  3  ^ 

C5  =  1.18  X  lO'Vcm'’. 
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directly  beneath  the  sample.  The  vs.  T  variation  for  selected 
samples  is  given  in  Figure  4.14  to  4.19.  The  solid  lines  represent  the 
theoretical  prediction  of  the  Rr  vs.  T  response  using  the  values  of  Cj 
and  i>p  as  given  in  Table  4.5.  The  fit  is  fair  and  the  trends  are 
clearly  indicated.  More  discussion  will  be  given  in  Chapter  V. 


I-V  Response  Oata  and  Contact  Parameters 


M 


CJ 


Figure  4.15  Temperature  variation  of  R^..  Sample  it  6-2D. 
Theoretical  fit  parameters  same  as  for 
Figure  4.9. 
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Figure  4.16  Temperature  variation  of  R  .  Sample  #  6-3B. 
Theoretical  fit  parameters  same  as  for 
Figure  4.10. 


103 


SECTION  V 


DISCUSSION 

5.1  Contact  resistivity  and  Sn  dopant 
surface  concentration 


When  this  research  was  initiated,  it  was  felt  that  experimental 
measurements  of  Rc  on  samples  with  varying  Cj  would  yield  a  curve 
corresponding  to  a  single  41B  line  in  Figure  4.5.  This  was  because 
previous  researchers  had  found  strong  evidence  of  Fermi  level  pinning 
in  Au-GaAs  Schottky  barriers  resulting  in  a  relatively  invariant  (j>B  “ 
0.9  eV. 

The  experi^mental  data  for  R^,  plotted  against  the  values  of  Cg 
determined  from  p  measurements  using  Figure  3.3,  do  not  line  along  a 
single  ((ib  line.  The  scatter  in  data  points  could  be  due  to  two 
factors.  One  being  that  the  effective  (jio  is  not  pinned  and  is  in  fact 
varying.  The  other  is  that  the  values  of  Cg  as  determined  from  p 
data  are  not  accurate. 

To  resolve  this  uncertainty,  contact  I-V  response  was  analyzed 
and  yielded  the  data  of  Table  4.5.  Immediately  apparent  from  this 
table  is  the  fact  that  the  values  of  ((ig  appear  to  be  varying.  They  are 
also  somewhat  larger  than  the  0.9  eV  usually  associated  with  the 
Au-GaAs  contact.  Also  evident  is  the  fact  that  the  values  of  Cs  are 
greater  than  the  values  of  Cj  determined  by  p  measurements. 

For  this  research  the  value  of  Cj  is  one  of  the  important 
parameters  in  contact  characterization  using  the  tuneling  model. 
Therefore,  various  methods,  as  explained  in  Section  IV,  were  used  to 
try  to  determine  an  accurate  value  and  to  provide  cross  checks  of  that 
value.  However,  as  seen  by  the  Cs  data  in  Table  4.4  and  Table  4.5,  the 
Cs  values  calculated  using  various  methods  are  quite  different  for  the 
same  sample. 

To  determine  which  Cs  value  is  the  most  accurate,  it  is 
necessary  to  consider  the  manner  in  which  Cs  was  calculated  in  each 
case.  The  use  of  Eq.  (4.2)  with  D  given  by  either  Di  =  2.99  x  10”^ ^ 

D  =  3.056  X  10"^**  ,  or  03  =  9.066  x  10’^“*  cm^/sec  resulted  in  the 
widest  variation  in  Cs-  This  is  to  be  expected,  since  the  value  of  Cs 
is  quite  sensitive  to  the  choice  of  diffusion  coefficient,  D. 
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Since  the  values  of  C5  as  detennined  the  the  p  measurements  were 
in  the  right  range,  relative  to  the  I  -V  data,  the  values  of  C5  as 
calculated  using  D  =  D^,  or  D  =  D2  appear  to  be  inaccurate.  Values  of 
Cs  in  the  10^^/cm^  range,  as  obtained  using  Di,  are  clearly  too  low  to 
allow  appreciable  tunneling  to  take  place.  Values  for  C5  of  10^^  to 
10^^/cm^,  as  obtained'using  D2»  seem  too  high  since  such  values  are 
above  the  atomic  concentration  of  GaAs.  Also  the  Sn  in  the  dopant  film 
was  only  on  the  order  of  10^°/cm^.  For  these  reasons,  the  use  of 
and  D2  in  determination  of  C5  is  not  appropriate. 

The  C5  values  obtained  using  0=03  and  the  values  of  C5 
obtained  from  p  measurements  are  in  slightly  better  agreement.  Also,  a 
comparison  of  the  Cs  values  determined  from  I-V  data  to  those  obtained 
using  0=03  show  fair  agreement.  The  values  of  Cs  as  determined  from 
the  above  three  methods  are  plotted  in  Figure  5.1  against  the  value  of 
p  as  determined  from  van  der  Pauw  measurements.  This  value  of  p  is 
valid  since  it  was  verified  during  the  transfer  length  measurements. 

The  data  are  also  listed  in  Table  5.1. 

It  is  evident  from  Figure  5.1  that  the  data  for  Cs  from  D3  and 
that  from  the  I-V  measurements  appear  quite  scattered.  However,  in 
general  the  Cs  values  from  the  I-V  measurements  lie  above  those  from 
the  p  measurements  and  the  Cs  values  from  the  0=03  calculation  are 
the  highest. 

Using  only  Figure  5.1  it  would  be  impossible  to  determine  which 
Cs  is  the  most  accurate  value.  Therefore,  a  measurement  of  Rq  versus 
temperature  wa«  made  in  order  to  determine  which  value  of  Cs  was  the 
most  accurate.  As  seen  from  the  Rc  vs.  T  plot.  Figures  4.14  to  4.19, 
the  use  of  the  Cs  and  <(ib  values  as  obtained  from  the  I-V  response  give 
acceptable  fits.  The  fit  is  actually  better  than  it  first  appears 
because  the  change  in  R^-  for  a  given  choice  of  ijiB  and  Cs  is  very 
sensitive  and  therefore  by  only  a  slight  change  in  the  values  of  Cs  or 
90  used  in  the  calculation,  the  experimental  and  calculated  points 
could  be  made  to  fit  even  more  closely.  The  point  is  that  the  Cs 
values  as  determined  from  the  I-V  curve  fit  do  accurately  predict  the 
Rc  vs.  T  response,  and  therefore  these  Cs  values  should  be  considered 
to  be  the  most  accurate. 

With  this  as  a  basis,  the  values  of  Cs  as  calculated  from  Eq. 
(4.2)  with  0=03  are  too  large,  meaning  that  the  actual  D  is  not  quite 
equal  to  D3.  Also  the  values  of  Cs  as  determined  by  a  measurement  of  p 
are  too  small.  However,  in  order  for  these  conclusions  to  be  valid  the 
question  of  the  scatter  in  the  data,  as  seen  in  Figure  5.1,  must  also 
be  addressed. 

The  scatter  in  the  Cs  vs.  p  data  may  be  partly  due  to  inaccuracy 
in  the  value  of  p  against  which  the  Cs  data  was  plotted.  The  p  values 
were  determined  from  a  van  der  Pauw  measurement  as  explained  in  Section 
III.  However,  the  van  der  Pauw  measurement  actually  yields  a  value  of 
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Table  5.1 


Values  of  Cs  Determined  Using  Various  Methods  of  Calculation 


Sample  # 

p  D-cm 

X  10-3 

Cg  rom  p 

/cm^  X  10^^ 

Cj  from  D3 
/cm}  X  10^^ 

Cg  from  I-V 
/cm3  X  10^® 

6-lB 

1.16 

7.7 

8.60 

16.8 

6 -2D 

1.77 

4.0 

2.70 

5.45 

6-3B 

2.33 

2.7 

2.57 

4.68 

7A-1B 

2.46 

2.5 

1.81 

9.64 

8A-1A 

3.99 

1.4 

1.17 

11.8 

4A-1-28 

1.35 

6.5 

1.25 

16.5 

5-3B 

2.06 

3.2 

5.78 

6.19 

5-4B 

2.13 

3.0 

3.50 

8.60 

6-4C 

1.78 

4.0 

9.18 

6.23 
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p  -  CM 


Figure  5.1  Variation  of  vs.  'p.  determined  using  various 
calculations.  ^  determined  from  van  der  Pauw 
measurements. 


resistivity,  Rj  ,  from  which  p  is  calculated  if  the  thickness  of  the 
conducting  layer  is  known.  Therefore,  p  was  not  determined  directly, 
but  rather  through: 

p  ■  ^s*j  (5*1) 

where  xj  was  the  layer  thickness  (junction  depth).  The  measurement  of 
Rj  is  quite  accurate,  as  was  the  measurement  of  xj  using  the  interfer¬ 
ometer.  However,  in  order  to  delineate  the  junction  a  lapping  and 
staining  technique  was  used  and  some  degree  of  variability  is  inherent 
with  this  step. 

Therefore,  if  the  value  of  xj  was  somewhat  in  error  due  to 
inaccurate  junction  delineation,  use  of  a  p  value  in  determining  C5 
through  Figure  3.3  would  result  in  errors  in  Cg.  From  the  van  der  Pauw 
measurement,  an  error  of  10%  in  Xj  would  give  an  error  of  10%  n  p  . 
Using  Figure  3.3,  this  error  of  10%  in  p  would  then  give  an  error  of 
about  10%  in  Cs- 

An  error  in  the  value  of  xj,  or  even  03,  greatly  affects  the 
value  of  Cs  as  calculated  using  Eq.  (4.2)  because  of  the  sensitivity  of 

the  erfc  term.  An  error  of  10%  in  either  xj  or  /  D3  results  in  errors 
in  Cs  on  the  order  of  80%.  Therefore,  the  scatter  in  these  data  points 
i s  even  greater. 

If  the  data  points  for  samples  7A-1B  and  8A-1B  are  neglected  in 
Figure  5.1,  due  to  probable  errors  in  Xj,  the  other  values  show  a 
general  trend  of  decreasing  Cs  for  increasing  p  .  If  the  curve  from 
Figure  3.3  is  shifted  up  by  approximately  a  factor  of  two  in  Cs  for  all 
points,  then  a  crude  fit  is  obtained.  Given  the  previously  mentioned 
difficulty  in  accurate  determination  of  xj,  the  scatter  around  this 
shifted  line  appears  to  be  a  more  reasonable  experimental  data  point 
scatter. 

The  set  of  Cs  data  points  determined  from  the  I-V  measurements, 
taken  as  being  the  most  accurate  and  lying  consistently  above  the  p-Cj 
curve  and  following  the  same  trend,  implies  that  the  basis  concepts 
used  in  the  calculation  of  the  p-Ce  curve  are  valid.  Therefore,  if 
slight  changes  were  made  in  the  calculation,  the  p-Cg  curve  could  be 
shifted  upward  and  be  made  to  correspond  well  to  the  experimental  data. 


One  change  could_lie  with  the  assumption  of  full  ionization. 

The  calculation  of  the  p-Cs  curve  was  based  on  the  assumption  that  the 
dopant  had  an  erfc  profile_and  that  all  donors  were  ionized.  An 
electrical  measurement  of  p  detects  only  the  free  carrier  concentration 
but,  under  this  assumption,  the  free  carriers  have  the  same  profile  as 
the  dopant  atoms  and  therefore  a  measurement  of  p  along  with  the  use  of 
Figure  3.3  will  give  an  accurate  value  of  the  dopant  concentration,  Cg. 
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For  trte  case  of  partial  ionization,  though  not  as  incomplete  as 
that  implied  by  Figure  3.2,  the  number  of  free  electrons  is  less  than 
the  donor  dopant  concentration.  Therefore,  for  a  given  erfc  donor  atom 
distribution,  the  value  of  p  with  partial  ionization  would  be  higher. 

Or  alternatively,  for  a  given  p  value,  the  value  of  C5  must  be  larger 
in  the  case  of  partial  ionization,  since  not  all  donors  contribute  a 
free  electron. 

This  is  exactly  what  is  seen  in  Figure  5.1.  For  a  given  value 
of  p  ,  the  value  of  Cj  as  determined  from  I-V  data,  lies  above  the 
value  of  C5  predicted  by  the  P-C5  curve  wherein  full  ionization  was 
assumed.  It  should  be  pointed  out  that  the  I-V  measurements  yield 
a  response  based  on  the  characteri sites  of  the  MS  contact  space  charge 
region.  In  this  region,  all  the  donors  are  ionized  due  to  the  band 
bending  caused  by  ({ig  and  therefore  the  value  of  Cg  calculated  from  IV 
data  refers  to  donor  atom  concentration  directly  and  not  just  bulk 
ionized  donor  concentration. 

The  fitting  of  the  I-V-C5  data  in  Figure  5.1  to  a  p-Cg  curve  for 
incomplete  ionization  would  follow  a  similar  calculation  to  that  used 
to  obtain  Figure  3.2.  The  degree  of  partial  ionization  would  not  be 
the  same  but  would_be  closer  to  full  ionization  since  the  difference  in 
I-V-Cg  values  and  p-Cg  values  is  only  about  a  factor  of  two. 

One  problem  with  the  above  argument  is  that  the  reasoning 
applied  to  reject  the  use  of  Figure  3.2  in  favor  of  full  ionization 
applies  here  also.  For  high  levels  of  doping  impurity  band  formation 
and  band  tailing  are  highly  likely  and  therefore  full  ionization  is  the 
case.  However,  an  upward  shift  of  the  p-Cg  curve  to  fit  experimental 
data  can  also  be  realized  by  a  change  in  the  assumption  concerning 
carrier  mobility. 

The  calculation  of  Figure  3.2  was  based  on  Eq.  (3.9)  which  used 
a  value  of  electron  mobility,  un  (9)-  These  values  of  ^0^ 

a  given  carrier  concentration  Nq  were  for  bulk  GaAs  samples.  In  this 
research  the  conduction  was  in  an  n-type  diffused  layer  instead  of  bulk 
material.  Therefore,  it  is  possible  that  Sze's  values  are  somewhat 
higher  than  is  actually  the  case.  For  diffused  areas  in  which  a  large 
amount  of  dopant  has  been  incorporated,  as  is  the  case  here,  the 
mobility  of  the  electrons  is  very  likely  less  than  that  for  bulk 
material  due  to  stress  to  the  lattice  incurred  during  the  diffusion. 

If  the  actual  mobility  is  less  than  that  represented  by  Eq.  (3.16), 
then  to  obtain  a  given  value  of  p  ,  more  electrons,  meaning  a  larger 
Cg,  are  needed  than  for  the  case  of  Figure_3.3.  The  result  is  to  shift 
the  p  -Cg  curve  upward  for  all  values  of  p. 

However,  since  the  values  of  I-V-Cg  in  Figure  5.1  are  about  a 
factor  of  two  larger  than  the  values  from  the  p-Cg  curve,  a  mobility 
decrease  by  a  factor  of  two  is  needed  to  give  twice  the  Cg  for  a  given 
p.  A  very  rough  check  can  be  obtained  by  comparing  the  values  of 


112 


average  mobility,  uh,  and  average  carrier  concentration,  ,  as  listed 
in  Table  4.3,  to  those  from  Sze.  From  Sze's  graph,  in  the  range  Nq  » 

2  X  10^^/cm^  ,  the  mobility  is  about  2500  cm"^/V-sec.  While  this  is  not 
exactly  twice  the  values  of  in  Table  4.3,  it  is  significantly 
larger.  Therefore,  an  upward  shift  in  the  p-Cj  curve  could  be  realized 
by  assuming  a  lower  value  of  mobility  for  the  electrons  in  the  diffused 
layer.  This  would  account  for  much  of  the  disagreement  between  the  p- 
Cs  and  I-V-C5  values  in  Figure  5.1  and  would  still  be  within  the 
assumption  of  full  ionization. 

If  impurity  band  conduction  is  occuring  due  to  the  high  doping 
levels  then  it  is  probable  that  the  electron  impurity  band  mobility  is 
less  than  the  normal  conduction  band  mobility.  However,  Sze's  graph 
actually  takes  this  into  account  indirectly  since  his  plot  is  with 
respect  to  actual  measured  carrier  concentrations.  The  plotted 
mobilities  are  experimentally  observed  values  at  that  particular  doping 
level.  The  only  difference  in  the  case  of  this  research  is  that 
diffused  layers  instead  of  bulk  material  are  used. 

A  final  check  on  the  validity  of  the  I-V-Cg  data  and  the 
assumption  of  full  ionization  is  in  relating  the  I-V-Cg  data  and  xj  to 
Eq.  (4.2).  From  Figure  5.1  it  is  evident  that  the  actual  diffusion 
coefficient  is  not  equal  to  D3,  given  by  Eq.  (4.14),  since  the  0$ 
values  predicted  using  0  =  O3  are  too  large.  If  the  I-V  data  values  of 

Cs  from  Table  4.5  are  used  with  values  of  xj  from  Table  4,3  and  Cg  from 
Table  4.4,  a  value  of  appropriate  diffusion  coefficient  can  be 
calculated  using  Eq.  (4.2).  The  results  of  such  a  calculation  are 
presented  in  Table  5.2.  The  diffusion  time,  t,  was  taken  as  6  hours. 
The  diffusion  temperature  was  950°  C. 

From  Table  5.2  it  is  evident  that  an  assumption  of  full 
ionization  and  an  erfc  dopant  profile  does  give  consistent  results. 

The  value  of  diffusion  coefficient,  0  «  1.5  x  10"^ ^om^/sec,  is  very 
similar  for  all  samples.  This  value  of  D  is  slightly  larger  than  D3  = 
9.066  X  10”^‘*  but  is  in  the  same  range. 

The  value  of  D  for  sample  7A-1B  is  unusually  low,  indicating  a 
possibility  of  error  in  the  measurement  of  xi.  The  values  of  D  for 
sample  SA-IB  and  4A1-2B  are  close  to  each  other  and  to  D3  but  are 
lower  than  for  the  other  samples.  This  could  be  due  to  some  dependence 
of  0  on  substrate  doping.  The  substrate  doping  for  the  above  two 
samples  was  less  than  for  the  others  by  about  a  factor  of  two. 

In  summary,  the  discussion  presented  in  this  section  shows  that 
a  Sn  diffusion  in  GaAs  can  be  described  by  an  erfc  profile  with  all 
donor  atoms_assumed  to  be  fully  ionized.  Also,  a  theoretically 
determined  p-Cg  curve,  relating  layer  resistivity  to  surface 
concentration,  is  applicable  if  a  correction  for  surface  effects  is 
made  for  the  values  of  electron  mobility  used  in  the  calculation. 
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Table  5.2 


Diffusion  Coefficient  for  Sn  in  GaAs  Calculated  Using  Cg  Values 
Determined  from  I-V  Measurements 


Sample  # 

*3 

Ml  crons 

Cb  /cm3 

X  1017 

Cc  /  cm^ 

X  1018 

D  cm^/sec 
X  10-13 

6-18 

1.91 

1.97 

16.8 

1.33 

6-2D 

1.68 

1.97 

5.45 

1.48 

6-3B 

1.67 

1.97 

4.68 

1.56 

6-4C 

1.44 

1.97 

6.23 

1.04 

5-3B 

1.83 

2.0 

6.19 

1.69 

5-4B 

1.73 

2.0 

8.60 

1.34 

7A-1B 

1.02 

1.98 

9.64 

0.45 

8A-1B 

1.64 

1,03 

11.8 

0.90 

4A-1-2B 

1.64 

1.1 

16.6 

0.84 

1-1 


5.2  Variation  of  ((ir 

Previous  research  on  MS  contacts  and  barrier  heights  has  led  to 
the  general  conclusion  that  an  Au-GaAs  contact  barrier  is  subject  to 
Fermi  level  pinning  with  a  resultant  relativity  constant  in  the 
range  0.9  to  1.0  eV.  The  exact  mechanism  of  pinning  are  not  clear  but 
the  effect  is  too  make  <(13  invariant  with  respect  to  doping  level  in 
semi  conductor. 

As  regards  the  research  presented  here,  the  above  implies  that 
the  variation  of  should  be  due  only  to  a  variation  in  with  a 
fixed  1(15  ■  0.9  eV.  The  I-V  data  as  given  in  Table  4.5  indicated  that 
(})B  was  also  varying,  although  most  values  were  near  1.0  eV.  However, 
use  of  these  unexpectedly  larger  41B  values  in  the  theoretical  equations 
resulted  in  good  fits  to  the  I-V  curves  and  quite  close  predictions  to 
the  Rc  vs.  T  response.  This  means  that  the  Schottky  barrier  model  is 
generally  valid  in  describing  contact  tunneling  properties  but  requires 
the  use  of  a  somewhat  larger  value  of  ^b  than  usually  quoted  in  the 
1 iterature. 

Since  there  is  very  much  experimental  evidence  supporting  Fermi 
level  pinning  in  Au-GaAs  barriers,  another  explanation  for  the  large  i^b 
was  sought.  The  most  probable  cause  is  due  to  the  effect  of  an 
additional  interfacial  oxide  layer  between  the  Au  and  the  GaAs  surface 
(71).  Some  researchers  have  suggested  that  this  oxide  may  be  formed 
during  the  contact  evaporation  step  (72).  The  effect  of  this 
additional  layer  is  to  reduce  the  value  of  the  saturation  tunneling 
current,  Jj  ,  and  therefore  to  present  an  effectively  higher  barrier. 

The  model  for  such  a  structure  is  shown  in  Figure  5.2.  This 
model  is  the  one  proposed  in  the  Bardeen  model  to  explain  Fermi  level 
pinning  (38).  The  actual  barrier  height  in  the  semiconductor,  ({>3  » 
can  still  be  subject  to  pinning  as  long  as  the  oxide  layer  is  thin 
enough  to  allow  the  metal  Fermi  level  to  affect  the  position  of  the 
semiconductor  Fermi  level  and  band  structure  through  the  difference  in 
work  functions.  For  a  barrier  thin  enough  to  allow  electron  tunneling 
to  occur,  the  above  assumption  is  quite  valid. 

This  means  that  the  effect  of  any  surface  states  in  pinning  (>3 
are  still  felt  and  that  a  413  ■  0.9  eV  would  be  formed.  The  band 
bending  in  the  semiconductor  and  the  shape  of  the  Schottky  barrier  are 
exactly  the  same  as  for  a  pure  MS  contact  with  the  same  4)3.  Therefore, 
electron  tunneling  through  the  space  charge  region,  the  Schottky 
barrier  portion  of  the  diagrams,  is  accurately  given  by  the  theoretical 
equations  presented  in  previous  chapters. 

The  effect  of  the  oxide  layer  is  to  present  an  additional  width 
of  insulator,  acting  as  an  effectively  wider  space  charge  region,  in 
the  path  of  the  tunneling  electrons.  This  effective  widening  of  the 
insulating  region  acts  to  reduce  the  value  of  contact  current  at  all 
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Figure  5.2  Model  of  a  MS  contact  with  an  interfacial 
oxide  layer. 


6 


applied  biases.  For  a  contact  I-V  characteristic  which  is  given  by  an 

equation  such  as  Eq.  (2.40)  or  Eq.  (2.55),  the  result  is  a  reduction  in 

the  Js  terms.  The  slope  of  the  I-V  response  in  the  exponential  region 
is  not  affected  as  long  as  the  barrier  is  thin  enough  so  that  the 
applied  contact  bias  is  taken  up  mostly  across  the  space  charge  region 
and  not  across  the  oxide  layer.  If  a  small  portion  of  the  applied  bias 
is  taken  up  across  the  oxide,  the  effect  is  seen  as  an  increase  in  the 
value  of  Eo  or  E*  in  Eqs.  (2.40)  and  (2.55). 

The  reason  that  an  effective  decrease  in  Jg  leads  to  an 
apparent  increase  in  (|ig  may  be  understood  by  considering  the  process  by 

which  ((ig  is  obtained  from  the  curve  fitting  process.  In  the  curve  fit, 

a  Schottky  barrier  model  without  an  interface  oxide  is  assumed.  As 
explained  in  Chapter  II,  this  gives  rise  to  an  I-V  response  as  given  by 
Eq.  (2.40)  or  Eq.  (2.55).  The  parameters  of  Jg  and  Eo  and  E*  are  then 
chosen  to  give  a  good  fit.  The  value  of  Cg  is  determined  from  the  E* 
term.  The  Jg  term  is  proportional  to  various  physical  constants  and 
also  to  the  value  of  E*  or  E®  through  an  ex(-i|)g/Eo)  term.  Therefore, 
once  Eo  and  Jg  have  been  determined,  a  value  of  4ig  is  chosen  to  fit 
the  theoretical  value  of  Jg,  through  Eq.  (2.41)  or  Eq.  (2.57),  to  the 
experimental  value. 

If  an  interfacial  oxide  layer  exists,  the  physical  effect  is  for 
the  contact  to  exhibit  a  value  of  Jg  lower  than  would  be  observed  with 
no  oxide  layer.  Since  in  the  Schottky  barrier  model,  a  lowering  of  Jg 
evidenced  by  an  increase  in  (jig  through  the  exp (-<)»g/Eo )  term,  an  experi¬ 
mentally  fit  value  of  Jg  lower  than  expected  would  result  in  a  value  of 
(fig  higher  than  expected,  as  observed  in  Table  4.5. 

An  estimate  of  the  thickness  of  the  interfacial  oxide  layer  can 
be  obtained  by  considering  the  reduction  in  Jg  from  that  expected  for 
barrier  with  <(ig  =  0.9  eV.  If  Fermi  level  pinning  results  in  (i)g  =  0.9 
eV  in  the  semiconductor  portion  of  the  contact  then  without  any  oxide 
layer  the  tunneling  electrons  would  experience  only  the  <t)g  =  0.9  eV 
Schottky  barrier  and  the  value  of  Jg  would  be  determined  by  Eqs.  (2.41) 
or  (2.57).  Since,  in  the  model  of  Figure  5.2,  it  is  assumed  that  the 
semiconductor  side  still  has  the  above  properties,  Eqs.  (2.41)  and 
(2.57)  would  still  predict  accurate  values  of  Jg.  However,  the  total 
contact  Jg  would  be  reduced  because  the  oxide  barrier  has  transparency 
less  than  100%. 

For  a  sample  with  a  given  Cg  and  Jg,  the  ratio  of  the  measured 
Jg  to  that  predicted  by  Eq.  (2.41)  or  (2.57)  for  (j>g  =  0.9  eV  with  the 
same  Cg  gives  an  indication  of  the  transparency  of  the  interfacial 
oxide  layer.  If  the  WKB  approximation  is  used  to  calculate  the 
transparency  of  a  rectangular  barrier  as  in  Figure  5.2,  an  estimate  of 
oxide  thickness  can  be  obtained. 

For  example,  using  data  from  Sample  6-lB,  the  curve  fit 
yielded: 
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Cs  =  1.68  X  10i5/cm3 

(?B  =  1.22  eV 

E«  =  8.43  X  10-2  eV 

Js  =  1.544  A/cm2 

under  Fe  at  T  =  23“  C.  Using  the  above  Cg  aind  (jig  =  0.9  eV,  tne 
theoretical  Jg  as  calculated  using  Eq.  (2.41)  is 

Jg  =  64.31  A/cm2 

Therefore,  the  ratio  of  the  objserved  Jg  ,  with  the  assumed 
interfacial  oxide,  to  the  value  of  Jg  for  no  oxide  and  ij)B  =  0.9  eV  is: 

Ratio  =  =  2.4  X  10-2  (5^2) 

This  should  equal  the  transparency  of  the  oxide  layer  at  the 
i nterface. 

From  a  WKB  approximation,  the  transparency  of  a  rectangular 
barrier  of  height  V(x)  is  given  by  (63): 

yO 

T  =  exp  (  f  [  V(x)-E  ]  dx  }  (5.3) 

h  xi 

where  xj  and  X2  af'e  the  points  where  V(x}  =  E.  For  this  case  E  =  np  = 
Fermi  level  is  te=aken  as  the  reference  energy.  For  a  constant  oxide 
barrier  height  measured  with  respect  to  the  Fermi  level, 

V(x)  -  E  =  q<(»ox  »  integral  becomes: 


X2 

I 

xl 


[  V(x)-E  dx 


(5.4) 


where  tgx  =  X2-^l  width  of  the  oxide  layer.  Substituting  Eq. 

(5.4)  into  Eq.  (5.3)  and  equating  this  to  Eq.  (5.2)  gives  the  result: 

♦ix^  ’^ox  “  1.379  X  10-Mev)^/2  m  (5.5) 

Therefore,  if  an  approximation  of  »  2  eV  is  used: 

^ox  -  9,75  A 

or  the  oxide  thickness  is  only  on  the  order  of  10  A  or  about  two  mono- 
layers.  For  other  estimates  of  <{iox  (5.5)  the  value  of  tgx  will 
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change  but  the  order  of  magnitude  of  tox  is  still  only  a  few  mono- 
1 ayers. 


It  is  very  likely  that  an  oxide  of  this  thickness  could  form  on 
the  surface  during  contact  processing  before  the  metal  evaporation. 
Therefore,  the  observance  of  410  values  slightly  greater  than  0.9  is 
probably  due  to  the  presence  of  a  thin  interfacial  layer  that 
originated  from  an  oxide. 

The  application  of  the  simple  Schottky  barrier  model,  with  the 
increased  40,  to  the  prediction  of  R^-  vs.  T  is  thus  reasonable.  It 
must  be  remembered  that  the  values  of  Rc  are  measured  near  zero  bias 
and  therefore  the  electrons  are  tunneling  through  a  nearly  unperturbed 
barrier  at  energies  about  0.9  eV  from  the  top  of  the  barrier.  This 
tunneling  is  affected  mostly  by  the  thickness  of  the  barrier  and  not  by 
its  shape.  Therefore,  an  ideal  Schottky  barrier  with  increased  40  or  a 
two  piece  barrier  as  shown  in  Figure  5.2,  with  the  same  effective 
tunneling  width,  should  give  very  similar  results.  As  seen  in  the 
plots  of  Rc  vs.  T,  the  use  of  an  increased  410  in  the  standard  Schottky 
barrier  model,  although  not  absolutely  physically  accurate,  can  still 
give  quite  good  predictions  of  contact  performance  since  the  increased 
40  does,  to  first  order,  account  for  the  additional  tunneling  barrier 
thickness. 


5.3  General  applicability  of  the  Schottky 
barrier  tunneling  model 

From  the  fair  curve  fits  obtained  to  experimental  I-V  data,  it 
is  apparent  that  the  use  of  the  Schottky  barrier  model  and  the 
resulting  I-V  equations  for  FE  and  TFE  type  tunneling  can  characterize 
contact  behavior  sufficiently  accurately  for  contact  modeling.  The  use 
of  the  model  is  further  supported  by  the  fact  that  the  parameters 
obtained  from  the  I-V  curve  fit,  Cj  and  4)0  ,  can  adequately  predict 
contact  resistivity  with  respect  to  another  independent  variable, 
tanperature,  as  seen  in  the  plots  for  R^  vs.  T. 

The  variables  €5  and  410  ,  as  calculated  from  the  use  of  the 
model,  can  be  considered  to  be  accurate  since  they  could  be  verified  by 
other  measurements.  The  values  of  Cs  were  verified  by  a  calculation  of 
Cs  based  upon  resistivity  measurements  on  the  diffused  layer.  However, 
as  explained  before,  it  was  necessary  that  some  additional 
consideration  be  given  to  the  estimation  of  carrier  mobility  in 
developing  the  p  -Cs  curve.  The  values  of  4)3  were  verified  from  the 
data  of  other  researchers  wherein  4)0  «•  0.9  eV  was  given  as  an  expected 
value.  Values  of  4)0  obtained  from  the  curve  fit  data  in  this  research 
were  somewhat  greater  than  0.9  eV  but,  as  was  explained  in  section  5.2, 
the  discrepancy  arose  from  the  presence  of  thin  interfacial  oxides  and 
the  manner  in  which  4)3  was  obtained  from  the  curve  fitting  process. 
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This  model  then  implies  that  to  obtain  non-alloyed,  Au-GaAs 
contacts  with  "ohmic"  properties,  i.e.,  in  the  10"^  s2-cm^  range,  FE 
type  tunneling  must  be  present.  This  then  puts  very  specific 
requirements  on  the  values  of  ijig  and  Cs  which  would  be  necessary.  If 
an  actual  (|)g  of  0.9  eV  could  be  realized,  meaning  necessary  care  in 
processing  to  eliminate  any  interfacial  oxides  or  insulating  layers, 
then  values  of  Cs  in  the  lO^^/cm^  range  are  required.  This  has  further 
implication  for  the  type  of  processing  steps  which  would  be  required  to 
obtain  this  level  of  Cs. 

5.4  Contact  resistivity  of  Au-Sn  doped  GaAs  contacts 

For  the  non-annealed  Au-Sn  doped  GaAs  contacts  studied  in  this 
research,  the  minimum  value  of  Rc  obtained  was  in  the  10”  n-cm^  range. 
When  the  contacts  were  analyzed  using  the  Schottky  barrier  model,  the 
main  limiting  factor  preventing  the  realization  of  lower  R^'s  was  the 
apparent  limitation  in  Cs.  Incorporating  Sn  into  the  GaAs  surface  by 
diffusion  produces  a  maximum  Cs  of  only  on  the  order  of  1  ~  2  x  10^^/ 
cm2,  -piis  means  that  even  if  the  surface  oxides  which  gave  an  apparent 
(fig  “  1.2  eV  had  been  removed,  allowing  a  larger  Js  and  the  more  normal 
Og  “  0.9  eV  to  be  realized,  the  value  of  R,-  would  still  be  only  in  the 
10-3  _  lo-**  ii^cm^  range.  This  is  nearly  two  orders  of  magnitude  larger 
than  the  values  of  R^  reportedly  being  obtained  for  various  alloying 
methods. 

Solid  solubility  data  for  Sn  in  GaAs  are  not  readily  available 
but  other  researchers  (17,  44)  give  values  of  maximum  Cs  for  Sn 
diffused  layers  in  the  high  lO'^^/cm^  low  10^^/cm^  range.  If  this 
represents  the  maximum  Cs  due  to  a  solid  solubility  limit,  then  R^ 
values  lower  than  the  10"2  a^cm^  observed  here  are  not  possible  using 
Sn  diffusion.  Other  possible  dopant  materials  will  also  be  under  solid 
solubility  restrictions  and  if  they  result  in  only  ■»  10^ ^/cm^  values 
for  Cs  then  their  use  will  also  not  result  in  lower  values  for  R^. 

The  above  considerations  mean  that  if  a  non-alloyed  contact  is 
desired,  then  a  method  to  dope  the  GaAs  surface  to  Cs  »  lO^O/cm^  is 
necessary  to  get  R^  -  10-^  n-cm^.  It  appears  that  a  predeposition 

type  of  diffusion  of  dopant  is  unsuitable  since  solid  solubility 
effects  are  likely  to  limit  Cs  to  less  than  the  102°  range. 

Some  very  recent  research  (30)  has  produced  low  R^  (10-®  n-cm^ 
range)  contacts  using  an  initial  diffusion  of  Sn  but  including  a 
subsequent  laser  treatment  of  the  diffused  surface  before  application 
of  the  contact  metal.  Most  likely,  the  laser  treatment  results  in  a 
localized  melting-regrowth  condition  on  the  surface  and  under  such 
conditions  dopant  segregation  during  regrowth  occurs.  For  Sn  in  a  GaAs 
melt,  the  Sn  collects  in  the  melt  during  regrowth  and  since  the  surface 
would  be  last  to  solidify,  the  Sn  would  be  collected  at  the  surface  and 
result  in  the  large  Cs  needed  for  low  values  of  Rq.  Another 
possibility  is  that  the  laser  treatment  allows  tin  atoms  to  go  from 
interstitial  to  substitutional  sites. 
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One  problem  with  creating  a  surface  with  Cs  ■»  10^°/cm^  is  that 
the  dopant  concentration  is  becoming  an  appreciable  percentage  of  the 
semiconductor  material.  This  can  result  in  large  changes  in  band 
structure  and  semiconductor  properties  in  general.  The  modeling  of  the 
MS  contact  would  become  extremely  difficult  in  such  cases  and  little 
could  be  quantitatively  explained  about  contact  performance  other  than 
that  it  exhibited  a  low  value  of  R^. 

Alloying  a  contact  metal  into  a  GaAs  surface  to  create  a  low  R,; 
contact  is  also  presumed  to  work  according  to  some  type  of  dopant 
redistribution  and  regrowth.  The  problem  is  that  trying  to  fit  the  FE 
or  TFE  model  to  such  contacts  is  difficult  when  the  fabricated  R^  is 
very  low.  For  such  cases  the  expected  exponential  variation  of  I  and  V 
is  seen  only  for  I  values  larger  than  those  appropriate  for  useful 
experimentation  and  therefore  effective  values  of  $8  and  Cg  cannot  be 
determined.  Also,  without  suitable  knowledge  of  iJib  Cj,  tunneling, 
if  it  occurs,  cannot  be  characterized. 

In  summary,  the  verification  of  the  FE  and  TFE  Schottky  model 
implies  that: 

1.  For  non-alloyed  MS  contacts  with  large  Cg,  TFE  or 
FE  tunneling  is  responsible  for  the  low  R^  values 
and  "ohmic"  behavior  and  contact  performance  can 
be  predicted  using  the  appropriate  tunneling  equa¬ 
tions. 

2.  To  obtain  extremely  low  R^  values  (10"^  n-cm^ 
range)  for  non-alloyed  contacts,  Cg  must  be  very 
large,  in  the  lO^'^/cm^  range,  and  such  a  value  is 
not  likely  to  be  attained  by  a  diffusion  type  of 
process.  Therefore,  another  type  of  doping  pro¬ 
cess,  perhaps  such  as  ion  implantation  or 
subsequent  heat  treatments  of  the  surface,  along 
with  their  inherent  complications,  must  be 

uti 1 i zed. 

3.  The  processing  involved  in  making  non-alloyed 
contacts  is  critical  in  regard  to  the  surface 
preparation  of  the  semiconductor  just  before  the 
application  of  the  metal.  Interfacial  oxides  can 
adequately  be  modeled  by  an  increase  in  apparent 
(pg  and  if  present  the  increased  (jjg  results  in 
increased  values  of  Rc  and  poor  contact  proper¬ 
ties. 
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SECTION  VI 


APPLICATION  OF  THE  SCHOTTKY  MODEL 
TO  ANNEALED  CONTACTS 


6.1  Au-Sn  diffused  layer  annealed  contacts 

Having  shown  in  the  previous  section  that  the  Schottky  barrier 
model  and  the  associated  FE  and  TFE  tunneling  equations  can  adequately 
characterize  non-annealed  (non-al loyed)  contacts,  the  model  was  then 
applied  to  attempt  to  analyze  annealed  structures.  Selected  samples  of 
the  Au-Sn  diffused  GaAs  contact  structures  were  given  sequential 
annealings  in  flowing  N2  at  temperatures  to  425“  C.  The  annealing 
step  was  performed  for  5  minutes  at  each  temperature. 

After  each  anneal,  I-V  curves  were  taken  in  the  same  manner  as 
previously  explained.  Figures  6.1,  6.2,  and  6.3  show  the  results  of 
the  I-V  tests.  From  the  figures  it  Is  apparent  that  the  curves  have  an 
overall  exponential  shape  but  bend  over  due  to  the  effect  of  Rs  of  the 
conducting  layer,  which  doesn't  change  with  annealing.  The  effect  of 
the  successively  higher  temperature  anneals  is  to  push  the  I-V  curves 
upward. 


When  the  theoretical  equations  for  TFE  or  FE  tunneling  were 
applied  and  curve  fitted  to  the  experimental  data,  values  of  (})B  and  Cs 
were  determined  for  each  of  the  contact  I-V  curves.  Figure  6.4  gives 
the  variation  with  annealing  temperature  of  4>b.  as  obtained  from  the 
curve  fits.  Figure  6.5  gives  the  variation  with  annealing  temperature 
of  Cs . 


One  apparent  fact  is  that  for  T  up  to  275°  C  or  300°  C  both  (j)B 
and  Cs  are  seen  to  decrease  with  successive  anneals.  The  values  of  $3 
move  nearer  to  the  value  of  0.9  «  1.0  eV  normally  associated  with  a 
Au-GaAs  Schottky  barrier.  The  values  of  Cs  also  move  closer  to  the 
value  of  Cs  determined  from  p  measurements. 

This  decrease  can  easily  be  related  to  a  reduction  in  the 
thickness  of  the  interfacial  oxide.  As  mentioned  before,  the  presence 
of  an  interfacial  oxide  on  the  order  of  10  to  20  A  was  suspected  as 
resulting  in  apparent  experimental  values  of  (Jib  >  0.9  eV.  The  thinning 
of  this  oxide  layer  with  successive  anneals,  due  to  Au  migration  or 
other  effects,  results  in  an  increase  in  Js  since  this  additional 
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Figure  6.2  I-V  response  of  Au-Sn  diffused  sample  #  6-4C 
after  successive  anneals.  T  =  23°C. 
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Figure  6.3  I-V  response  of  Au-Sn  diffused  sample  #  7A-1C 
after  successive  anneals.  T  =  23®C. 


tunneling  barrier  is  decreased.  An  increase  in  Jg  is  reflected  in  the 
analysis  as  a  decrease  in  ((ig. 


The  decrease  in  apparent  value  of  C5  is  also  due  to  a  thinning 
of  the  oxide  layer.  With  a  thinner  oxide  layer,  less  of  the  applied 
contact  bias  is  dropped  across  the  oxide.  Therefore,  more  of  the 
applied  bias  appears  across  the  space  charge  region  and  gives  rise  to 
an  increase  in  current.  The  slope  of  the  exponential  region,  ,  is 
proportional  to  .  With  the  effects  of  an  increased  IR  drop  in  the 

conducting  layer^s  taken  into  account,  it  is  seen  that  as  T  increases, 
the  value  of  E'  decreases  which  implies  a  decrease  in  the  apparent  Cg. 

Past  the  300"  C  point,  (})g  and  especially  Cg  exhibit  erratic 
behavior.  The  value  of  4)g  increases  again  and  Cg  becomes  very  large. 
This  is  clearly  opposite  to  the  trend  up  to  this  point.  It  is  felt 
that  this  behavior  is  due  to  the  non-uniform  formation  of  contact 
pinholes. 

A  contact  pinhole  is  considered  to  be  a  small  area  in  which  the 
interfacial  oxide  has  been  essentially  reduced  to  zero.  At  this  point, 
a  near  ideal  MS  contact  is  made.  The  <1)B  for  such  points  should  be  near 
0.9  eV.  The  result  is  that  this  point  is  a  much  better  conducting 
point  for  current  flow  and  the  current  applied  to  the  contact  is  there¬ 
fore  concentrated  at  this  pinhole  point. 

The  increase  in  conductivity  at  these  points  results  in  a  large 
increase  in  Jg.  However,  the  localization  of  the  current  into  a  single 
entry  point  into  the  conducting  surface  layer  radically  alters  the 
model  of  the  contact  structure.  The  effect  is  to  insert  a  large 
additional  resistance,  due  to  Rg,  into  the  path  of  current  flowing  from 
the  pinhole  to  the  x  =  0  contact  edge.  The  increase  in  current  path 
resistance  flattens  out  the  exponential  portion  of  the  I-V  curve  due  to 
the  additional  IRg  drop.  Since  the  curve  fitting  procedure  is  based  on 
a  normally  distributed  Rg  and  uniform  contact  conduction,  the  flattened 
curve  is  interpreted  as  an  increase  in  E'  and  therefore  an  increase  in 
Cs. 


The  increase  in  413  is  related  to  the  increase  in  Cg  since  ^g  is 
obtained  by  a  calculation  involving  the  experimental  Jg  and  an 
exp(-  (jig/Eo)  term.  If  Jg  increases  then  4)3  should  decrease,  but  since 
an  effectively  higher  Cg  was  used  in  the  calculation,  an  effective 
increase  in  413  is  obtained  as  the  result. 

The  point  is  that  if  the  Schottky  model  is  used  in  determining 
413  and  Cg  when  contact  conduction  is  non-uniform,  misleading  informa¬ 
tion  results.  However,  the  presence  of  some  type  of  interfacial  oxide 
was  confirmed  and  the  implied  trend  is  that  with  ideal  processing  this 
interfacial  oxide  could  be  reduced  to  zero  uniformly  across  the  con¬ 
tact.  The  contact  would  then  exhibit  properties  related  to  conditions 
near  413  »  0.9  eV  with  Cg  near  the  value  determined  by  p  measurements. 
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One  important  consideration  for  the  above  is  that  the  annealing 
must  not  proceed  so  far  as  to  drastically  change  the  essentially 
non-alloyed,  abrupt,  nature  of  the  MS  interface.  It  was  found  that  for 
T  >  423°  C,  significant  melting  and  balling  up  of  the  Au  occurred. 
Presumably  this  would  also  affect  the  nature  of  the  interface  other 
than  by  merely  thinning  the  interfacial  oxide  to  zero.  The  end  result 
of  such  high  temperature  treatment  is  usually  more  highly  conducting 
contacts  with  higher  Jg  and  lower  R^-  but  the  mechanism  responsible  for 
the  better  properties  would  not  be  a  simple  reduction  in  the  (fig  of  a 
standard  tunneling  Schottky  barrier.  Furthermore,  at  temperatures 
above  400°  C  substantial  atomic  migration  and  departure  from 
stoichiometry  occurs  (26). 

The  effect  of  increasing  the  annealing  period  at  each 
temperature  was  also  examined  for  one  sample,  #5-4D.  In  this  case  the 
contacts  were  annealed  at  250°  C  for  increasing  periods  and  then  at 
300°  C  for  increasing  periods.  Figure  6.6  shows  the  I-V  response  of 
the  contact  after  the  various  annealing  temperature  and  times.  The 
maximum  temperature  was  limited  to  300°  C  to  avoid  the  previously 
mentioned  problem  with  possible  pinholes. 

Curve  fitting  of  the  theoretical  TFE  equations  was  done  and  the 
results  are  summarized  below,  and  also  presented  in  Figures  6.7  and 

6.8. 

1.  At  a  given  temperature,  subsequent  anneals  result  in  an 
upward  shift  of  the  I-V  curve,  in  the  same  manner  as  the 
higher  temperature  anneals  used  in  the  previous 
experiment.  Eventually  the  contact  seemed  to  stabilize  to 
a  final  Cg,  ifg  combination  for  a  given  temperature. 

2.  At  300°c,  Cg  and  <>b  appear  to  increase  from  the  250°  C 
data  due  to  the  probable  contact  nonuniformity  problems 
mentioned  earlier.  However,  further  heat  treatment  at  the 
same  temperature  seems  to  stabilize  the  change  in  Cg  and 
ijig  but  both  values  at  300°  C  are  higher  than  at  250°  C. 

The  implication  of  the  above  is  that  the  reduction  in  the 
interfacial  oxide  can  be  partially  accomplished  by  longer  annealing 
times  at  lower  temperatures.  However,  it  appears  that  the  effective 
elimination  of  the  oxide  depends  on  physical  changes  related  to  the 
magnitude  of  the  annealing  temperature  and  not  just  to  the  annealing 
time. 


6.2  Au-Ge-GaAs  annealed  contacts 


The  previous  contact  model  and  associated  FE  and  TFE  equations 
were  also  used  to  attempt  to  analyze  the  behavior  of  Au-Ge-GaAs 
contacts  after  successive  annealing-alloying  steps.  An  Au-Ge 
combination  is  the  most  widely  used  material  for  the  formation  of  a  low 
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Figure  6.6  I-V  response  of  Au-Sn  diffused  sample  #  5-4D 

after  successive  anneals  at  the  same  temperature 
but  for  increasing  times.  T  =  Z2°0. 
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Rc  contact  to  n-type  GaAs  using  the  alloying  method.  In  such  a  method 
a  contact  structure  is  delineated  on  the  GaAs  using  an  Au-Ge  eutectic 
alloy  or  a  Au-Ge-GaAs  layer  sandwich  and  heated  above  the  eutectic 
point  in  an  annealing  furnace  or  with  a  laser.  The  Ge  is  believed  to 
go  into  the  surface  forming  a  highly  doped  layer  which  results  in 
tuneling  and  a  low  R(;. 

For  this  experiment,  n-type  bulk  GaAs  substrates  with  the 
following  specifications  were  used. 

Si  doped 

p  ■  2.6  X  10-^  il-cm 

Cg  “  9  X  10^ ^/cm^ 

thickness  15  mils 

To  make  the  Au-Ge  contact  a  method  similar  to  Anderson,  et.al. 
(21)  was  used.  On  a  cleaned  wafer,  a  few  hundred  angstrom  layer  of  Ge 
was  evaporated  and  over  this  a  layer  of  Au  was  evaporated.  This 
differs  from  Anderson's  method  where  the  evaporated  Ge  layer  was  first 
formed  into  a  crystalline  epi  layer  by  substrate  heating  prior  to  the 
Au  evaporation.  After  the  contact  structure  was  delineated  using  a 
photoresist  lift-off  technique,  individual  strucutres  were  scribed 
apart  and  given  successive  anneals  and  tested. 

The  resulting  I-V  curves  from  selected  samples  are  shown  in 
Figures  6.9,  6.10,  and  6,11.  The  curve  fitting  of  the  theoretical 
FE  and  TFE  equations  to  this  data  is  much  simpler  than  previously  since 
the  IRs  drop  in  the  conducting  substrate  is  negligible  since  Rj  is  very 
small  being  -  0.607  iVsquare.  Therefore,  the  exponential  portion  can  be 
directly  fitted  to  a  straight  line  of  slope  1/E',  and  current  squeezing 
effects  in  the  contact  are  not  important.  For  the  As  deposited 
unannealed  curves,  the  I-V  response  bends  over  slightly  due  to  the 
effect  of  the  GeLayer  acting  as  a  slight  additional  series  resistance. 
For  these,  a  fit  was  made  to  the  lower  current  portion  of  the  curve 
where  IR  drops  are  negligible. 

Similar  to  the  annealed  contacts  discussed  in  Section  6.1,  these 
contacts  have  I-V  curves  which  shift  generally  upward  with  increasing 
annealing-alloying.  However,  in  many  cases  the  I-V  curve  shape  is  not 
strictly  exponenial  implying  that  alloying  is  modifying  the  interface 
characteristics  so  that  perhaps  a  simple  Schottky  barrier  cannot 
accurately  characterize  its  performance.  Curve  fitting  to  these  curves 
was  done  on  the  lower  portion  of  the  I-V  curve  above  the  knee.  In  all 
cases,  since  the  curves  shift  upward,  higher  values  of  Js  are  implied 
and  since  is  roughly  proportional  to  I/J5  increased  alloying  gives 
lower  values  for  Re* 

The  I-V  data  was  curve  fit  and  analyzed  to  give  values  of  (fg  and 
Cj  which  are  plotted  in  Figures  6.12  and  6.13.  The  (^g  and  Cj  data  for 
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igure  6.12  Variation  of  0g  for  Au-6e  samples  after 
successive  anneals. 
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Figure  6.13  Variation  of  for  Au-Ge  samples  after 
successive  anneals. 


tne  onjindl  unannealed  contacts  are  very  similar  for  dll  samples  and 
yield  a  Cs  value  whicii  matches  closely  with  the  original  specified  Cj. 
The  values  of  are  also  in  agreement  out  are  slightly  lower 

tnan  the  0.9  eV  standard  value. 

The  variation  of  C5  with  annealing  indicates  an  increase  in  0$ 
with  increasing  temperature.  This  is  to  be  expected  if,  as  usually 
assumed,  the  Ge  is  moving  into  the  GaAs  and  acting  like  an  n-type 
dopant  giving  rise  to  a  highly  doped  layer  and  increased  tunneling.  A 
maximum  temperature  of  =  425°  C  was  for  the  one  case,  used  since 
balling  up  of  tne  Au  was  observed  for  higher  temperatures.  The 
eventual  value  of  0$  reached  at  a  temperature  of  =  300°  C  was 
=•  5  ~  6  X  10* ^/cm-^  and  was  nearly  the  same  for  all  samples. 

One  problem  with  the  data  is  the  extraordinarily  high  values  of 
(Ptj  which  resulted  from  the  curve  fit.  They  are  clearly  out  of  line 
with  a  physically  valid  Schottky  model.  However,  it  is  felt  that  these 
high  values  are  due  to  the  presence  of  localized  conduction  pinholes 
and  arise  from  the  manner  in  which  the  data  is  interpreted  when  pinhole 
conduction  is  neglected. 

Other  researchers  (19)  nave  shown  that  alloyed  Au-Ge  contacts 
have  a  problem  with  localized  formation  of  various  Au-6e-GaAs  phases 
with  their  differing  conduction  properties.  A  conduction  pinhole  in 
this  case  may  be  regarded  as  a  localized  area  where  rapid  incorporation 
of  Ge  has  occured  giving  rise  to  a  highly  doped  point  in  the  GaAs 
surface.  This  point  then  allows  easier  tunneling  than  the  rest  of  the 
surface  and  such  points  would  effectively  act  as  a  current-hogging 
points.  Any  current  through  the  contact  would  essentially  be  confined 
to  these  easily  conducting  pinholes  and  the  measured  I-V  characteristic 
would  then  indicate  only  the  properties  of  the  pinholes  and  not  of  the 
entire  Au  contact  area. 

The  problem  concerning  the  data  analysis  arises  when  a  value  of 
Js  is  calculated  for  a  given  I-V  experimental  curve.  The  actual 
quantity  measured  is  a  value  of  saturation  current  1$  and  a  value  of 
saturation  current  density,  Js  ,  is  calculated  by  dividing  Is  by  the 
measured  contact  area.  However,  if  the  current  is  flowing  mostly 
through  isolated  pinholes  then  the  I-V  response  reflects  the  values  of 
Cs  and  (t>B  applicable  to  only  those  pinholes.  For  any  single  pinhole 
supporting  conduction  of  the  test  current  I,  the  value  of  Js  is  quite 
large  implying  a  much  smaller  4)8.  However,  if  the  measured  value  of  Is 
is  divided  by  the  total  contact  area,  the  main  portion  of  which  is  not 
supporting  conduction,  the  apparent  value  of  Js  is  reduced  which 
results  in  a  larger  apparent  value  of  ^ib- 

Pinhole  conduction  does  not  affect  the  calculated  value  of  Csas 
much  as  it  did  in  the  previous  experiment  on  annealed  Au  contacts  to  Sn 
diffused  layers.  The  reason  is  that  for  the  Au-Ge  contacts  the 
substrate  is  much  more  highly  conductive  and  the  IRs  drops  due  to  the 
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current  flowing  from  an  Internal  pinhole  are  negligible.  That  is, 
these  IRs  drops  do  not  cause  the  I-V  response  curve  to  be  bent  and 
therefore  the  slope  is  only  slightly  affected  by  the  doping  level. 

The  observed  high  values  of  <{ib  for  the  Au-Ge  case  are  probably 
not  due  to  the  effects  of  an  interfacial  oxide  as  was  assumed  to  be  the 
case  in  the  diffused  Sn  layer  case.  One  reason  is  that  the  unannealed 
contact  data  yielded  a  iJ>b  lower  than  expected.  The  presence  of  an 
oxide  should  have  made  this  value  greater.  Also,  the  observed  values 
do  not  progress  through  an  initial  period  of  decrease  in  ^b  ^nd  Cj  as 
was  seen  before.  Only  increases  are  observed.  The  evaporation  of  the 
Ge  appears  to  effectively  remove  any  effects  of  an  interfacial  oxide 
1  ayer. 


The  Schottky  barrier  model  can  be  applied  as  a  general  guide  to 
the  analysis  of  Au-Ge  alloyed  contacts,  but  some  strong  reservations 
must  be  made  in  the  interpretation  of  the  data,  as  regards  0$  and  0b» 
because  of  the  non  exponential  shape  of  the  I-V  response  and  the 
probably  formation  of  isolated  tunneling  points. 
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SECTION  VII 


OPTICAL  MEASUREMENTS 


The  argon  ion  laser,  dye  laser,  cavity  dumper,  and  active  mode 
locking  unit  arrived  in  October  of  1979.  The  assimilation  of  these 
components  into  a  synchronously  pumped  cavity  dumped  dye  laser  began  in 
late  November  and  was  completed  in  July  of  1980.  The  installation  of 
the  accousto-opti c  modelocker  causes  the  argon  ion  laser  to  have  a 
pulsed  output.  The  acousto-optic  modelocker  is  driven  by  a  stable  RF 
synthesizer.  This  RF  signal  produces  an  acoustic  standing  wave  in  the 
acousto-optic  crystal  when  the  RF  matches  the  resonance  frequency  of 
the  crystal.  Since  the  acoustic  wave  velocity  and  crystal  dimensions, 
which  determine  the  resonance  frequency,  are  temperature  dependent  a 
stable  temperature  controlled  oven  is  installed  around  the  crystal.  By 
modulating  the  cavity  losses  the  mode-locker  causes  approximately  1000 
modes  to  have  a  fixed  phase  relation.  The  interference  of  these  modes 
produces  the  pulsed  output.  The  pulses  are  ~  160  psec  in  length  and 
the  interpulse  spacing  is  12  nsec.  The  maximum  average  power  is  IW  but 
is  usually  run  at  600  mW  for  stability  and  pulse  length  considerations. 
Thus,  the  peak  power  is  45  W  when  the  wavelength  is  514.4  nm.  The  dye 
currently  being  used  in  the  dye  laser  is  Rhodamine  6G.  This  dye  will 
produce  output  in  the  range  from  550  to  650  nm  with  the  peak  power  near 
600  nm.  The  laser  system  is  termed  as  being  synchronously  pumped,  this 
is  achieved  by  lengthening  the  cavity  of  the  dye  laser  until  it  matches 
the  length  of  the  argon  ion  cavity.  Thus  the  round  trip  time  of  pulses 
travelling  in  both  cavities  is  the  same.  The  pulse  length  of  a 
model ocked  laser  is  determined  by  the  number  of  longitudinal  modes 
which  can  be  locked  together.  By  synchronously  pumping,  the  gain 
bandwidth  of  the  dye  laser  is  larger  at  the  time  tht  the  dye  laser 
pulse  hits  the  dye  jet.  Thus  the  dye  laser  pulse  can  be  and  is  much 
shorter  i.e.  pulse  lengths  on  the  order  of  1  psec  are  possible.  Peak 
power  in  these  pulses  extractable  through  a  10%  mirror  are  on  the  order 
of  5U  W  at  80  Mhz  repition  rate.  However,  also  installed  in  the  system 
is  a  cavity  dumper.  The  main  mechanism  of  the  cavity  dumper  is  an 
acousto-optic  crystal  which  like  the  modelocker  is  driven  by  an  RF 
signal.  However,  this  RF  signal  is  pulsed  producing  a  travelling  wave 
in  the  crystal.  This  RF  pulse  though  longer  than  a  round  trip  time  of 
the  cavity  is  synchronized  with  the  light  pulse  such  that  only  one 
pulse  is  dumped  out  per  RF  pulse.  Since  with  the  cavity  dumper  the  dye 
laser  cavity  is  defined  by  only  high  reflectance  mirrors  higher 
intracavity  power  is  possible  and  a  much  greater  fraction  of 
intracavity  pulse  power  may  be  switched  out  giving  a  factor  of  30 
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increase  in  output  pulse  peak  power,  i.e.  the  peak  power  is  1.5  kW  with 
a  repition  rate  up  to  4  Mhz.  The  average  power  with  the  cavity  dunper 
is  specified  to  be  25-3U  nW.  These  peak  powers  are,  with  optical 
focusing,  sufficient  to  do  optical  excitation  studies  in  solids. 


In  the  proposed  experiment  small  changes  in  the  reflectivity  of 
the  material  will  be  measured  by  a  technique  based  on  the  principles  of 
el  1 i psometry .  A  typical  block  diagram  of  an  ellipsometer  is  shown 

in  Fig.  1.  The  plane  polarized  light  from  the  laser  is  converted  to 
circularly  polarized  by  the  first  quarter  wave  plate,  this  ensures  that 
the  intensity  does  not  vary  as  the  polarizer  P  is  rotated.  The  plane 
polarized  light  emerging  from  P  then  passes  through  the  second 
retardation  plate  C  which  is  oriented  such  that  its  fast  axis  is  at  40° 
to  the  horizontal  plane,  which  is  also  the  plane  of  incidence.  The 
emerging  light  is  generally,  elliptically  polarized.  On  reflection 
from  the  sample  surface,  the  two  components  of  the  electric  field  of 
the  incident  light  Ej^  and  Eu  undergo  different  changes  in  their 
respective  amplitudes  and  phases.  The  reflected  light  passes  through 
the  analyser  A  and  is  detected  by  the  photomultiplier  tube  (PMT).  In 
the  experiments  P  and  A  are  rotated  to  such  positions  that  the 
elliptically  polarized  light  from  C  becomes  linearly  polarized  after 
reflection  from  the  sample  surface  and  can  then  be  extinguished  by  the 
analyzer  A.  In  this  configuration  the  angles  that  P  and  A  make  with 
the  plane  of  incidence  can  be  related  to  the  (complex)  refractive  index 
of  the  sample.  However,  in  the  proposed  experiment,  the  quantity  of 
interest  is  the  change  in  the  refractive  index  of  the  material  caused 
by  the  excitation  of  the  electron-hole  plasma.  It  can  be  shown 
that  for  small  changes  6n  in  the  refractive  index,  the  intensity  of 
light  transmitted  through  the  ellipsometer  arrangement  is  proportional 
to  (6n/n)^.  Thus  if  P  and  A  are  adjusted  for  extinction  in  the  absence 
of  the  pump  beam,  then  the  change  in  the  signal  reaching  the  PMT  when 
the  pui.i^,  beam  is  turned  on  can  be  used  to  calculate  the  change  in  the 
refractive  index  and  thus  the  change  in  the  electron-hole  plasma 
density.  This  plasma  density  can  be  measured  as  a  function  of  the 
delay  between  the  pump  pulse  and  the  probe  pulse  yielding  information 
about  the  relaxation  time  constants  of  the  plasma. 


The  sensitivity  of  this  procedure  as  reported  in  literature 
is  6n/n  <  10"^  and  an  initial  change  of  5%  gave  sufficient  information 
to  obtain  characteristic  time  constants  in  the  reported  experiment. 
Taking  into  account  the  coupling  of  the  pump  beam  energy  to  GaAs  and 
its  absorption  constant  at  the  laser  wavelength  (600  nm)  the  initial 
plasma  density  excited  in  the  material  can  be  estimated.  Using  the 
simple  Urude  model  to  relate  the  plasma  density  to  the  change  in  the 
refractive  index,  it  was  estimated  that  an  energy  density  of  5  mJ/cm^ 
would  be  required  to  produce  the  initial  5%  change  in  the  refractive 
index.  The  dye  laser  system  is  capable  of  providing  energy  of  30  nJ 
per  pulse.  Thus  in  order  to  obtain  the  required  energy  density  it  is 
necessary  to  focus  the  pump  beam  to  a  spot  with  dimensions  of  a  few 
microns.  It  is  also  desirable  that  the  probe  beam  should  probe  the 
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SAMPLE  C  P 


Figure  7.1  Block  diagram  of  the  ellipsometry  arrangement 
for  measuring  small  changes  in  the  reflectivity 
of  the  sample  as  a  function  of  delay  between  the 
pump  and  the  probe  beams. 


region  of  the  sample  where  the  pump  beam  excites  the  plasma.  This 
makes  it  necessary  to  spatially  overlap  the  two  spots.  Also  in  order 
to  calculate  the  time  constants  of  the  material  it  is  necessary  to  know 
the  zero-delay  position  of  the  delay  line. 

The  spatial  and  temperal  overlap  was  obtained  by  maximizing  the 
second  harmonic  signal,  generated  by  a  combined  effect  of  the  electric 
fields  of  the  two  beams  in  reflection  from  a  (110)  surface  of  GaAs. 

In  order  to  evaluate  the  sensitivity  of  the  PMT  (RCA  4837)  in 
the  uv  range  and  to  gain  familiarity  with  the  process  of  generating 
second  harmonic  light,  a  one  beam  experiment  was  initially  performed. 

A  block  diagram  of  the  experimental  arrangement  is  shown  in  Fig.  2. 

The  vertically  polarized  beam  was  reflected  from  the  sample  and  the 
reflected  primary  beam  was  filtered  out  using  two  uv-pass  glass 
filters.  The  transmitted  light  was  detected  by  the  PMT.  The  dye  laser 
beam  was  chopped  with  a  mechanical  chopper  which  also  provided  the 
reference  for  the  lock-in  amplifier. 

To  confirm  that  the  signal  thus  observed  is  due  to  the  300  nm 
second  harmonic  and  not  due  to  the  residual  transmittance  of  the  600  nm 
primary,  the  dependence  of  the  signal  on  the  angle  'V  between 
the  (10.0)  axis  in  the  plane  of  the  sample  and  the  plane  of  incidence 
was  studied,  w)  The  SH  signal  has  two  components,  one  polarized  in 
the  plane  of  incidence  and  the  other  polarized  perpendicular  to  the 
plane  of  incidence.  The  reflected  SH  intensity  can  be  calculated 
according  to 

I{2v)  U  (2v)  +  III  (2v)  (7.1) 

where  Ii  (2v)  is  proportional  to  (Pi’^^^)^  and  In  (2v)  to  (Pl|NI-S)2^  pNLS 
are  the  non-linear  polarizations  induced  by  the  electric  field  Ei  of 
the  primary  beam.  The  dependence  of  P*'**-^  on  T  is  given  by: 

Pj^NLS  =  3/2  X^L  e!  sin2  4,  ^os  'll  (7.2) 

and 

p,NLS  =  1/2  xNL  E?  sin  (1-3  cos^  y)  (7.3) 


Here  is  the  non  zero  term  in  the  second  order  non-linear 
susceptibility  for  GaAs. 

Figure  3  shows  the  experimental  data  of  I(2v)  as  a  function  of 
plotted  along  with  the  theoretical  curve  obtained  from  Equations  2 
and  3.  The  theoretical  curve  has  been  seal  led  to  account  for  the 
proportionality  factor.  Agreement  between  theory  and  experiment  is 
excellent  indicating  that  the  observed  signal  is  indeed  the  second 
harmonic. 
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Figure  7.2  Block  diagram  of  the  experimental  arrangement 

for  second  harmonic  generation  for  a  GaAs  crystal 
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ANGLE  FROM  (100)  AXIS  IN  DEGREES 


Figure  7.3  The  intensity  of  the  second  harmonic  signal  generated 
from  GaAs  as  a  function  of  the  angle  ¥  between  the  (100)  axis  and  the 
normal  to  the  plane  of  incidence  and  reflection. 


After  observing  the  single  beam  second  harmonic,  a  two  beam 
experiment  was  performed.  The  angle  63  at  which  the  SH  signal  appears 
is  given  by  ^ 


sin  63 


sin  01  +  sin  02 


2 


where  01  and  02  are  the  angles  of  incidence  of  the  two  beams. 

In  the  initial  experiment  a  beam  splitter  giving  approximately  equal 
intensities  for  the  two  beams  was  used.  6i  and  62  were  chosen  to  be 
of  opposite  sign  so  that  the  SH  signal  appeared  near  the  film  normal. 

A  block  diagram  of  the  arrangement  is  shown  in  Fig.  4. 

The  observed  signal  was  easily  confirmed  as  being  the  SH  by 
qualitatively  observing  its  dependence  on  'V  and  comparing  it  to  Fig.  3. 
To  further  evaluate  the  noise  levels  in  the  system  and  to  check  the 
alignment  of  the  delay  line,  this  SH  signal  was  measured  as  a  function 
of  the  delay  introduced  in  one  of  the  beams.  The  resulting  curve 
duplicated  the  auto  correlation  trace  obtained  from  the  commercial  unit 
(Spectra  Physics  409)  as  shown  in  Fig.  5.  The  actual  geometry  to  be 
used  in  the  ellipsometry  was  then  set  up  with  the  pump  beam  incident 
normally  and  the  probe  beam  with  an  angle  of  incidence  of  70®.  Due  to 
this  shallow  incidence,  positioning  the  spots  on  the  sample  became  more 
difficult  requiring  harder  focusing  and  smaller  lens  mounts.  A 
variable  beam  splitter  was  used  in  place  of  the  fixed  one  in  Fig.  4. 
After  finding  the  SH  signal  with  the  two  beams  of  approximately  equal 
intensity,  the  ratio  can  be  changed  while  tracking  the  SH  signal.  At 
the  sample  the  probe  beam  can  be  made  to  be  a  factor  of  twenty  weaker 
than  the  pump  beam. 

To  numerically  subtract  the  background  signal  caused  by  the 
probe  beam  scattered  into  the  PMT  and  by  the  residual  transmittance 
through  the  analyzer,  a  shutter  controlled  with  the  stepper  motor  is 
placed  in  the  pump  beam.  This  way,  for  each  psoition  of  the  delay  line 
the  output  of  the  PMT  can  be  measured  with  or  without  the  pump  beam 
striking  the  sample. 

In  any  laboratory  there  are  nwny  pieces  of  apparatus  necessary 
to  perform  the  experiments  successfully.  In  addition  to  the  standard 
needs  of  lens,  mirrors,  mounts,  and  positioners  several  other 
mechanical  and  electrical  devices  have  been  acquired,  built,  and 
investigated  for  use. 

Crucial  to  the  experiments  performed  is  the  operation  of  an 
optical  delay  line.  The  difficulties  are  mechanical  stability, 
reproducibility,  and  trueness  of  travel.  The  probe  beam  in  the 
experiments  passes  through  this  delay  line  and  then  another  meter  to 
the  sample.  At  the  sample  the  beams  are  focused  to  20  micron  diameter 
spots.  Thus  it  is  necessary  to  have  the  delay  line  translate  and  be 
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DEL/^y  TIME  IN  psec 


Figure  7.5  Top:  Second  harmonic  intensity  as  a  function  of  time  delay 
between  the  two  pulses  in  Figure  7.4. 

Bottom:  Auto  correlation  trace  produced  by  a  commercial 
system  which  employs  second  harmonic  generation  by  trans¬ 
mission  (20  psec/div.  horizontal). 
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aligned  very  true  in  order  to  keep  the  two  spots  overlapped.  Also 
tnis  delay  line  is  coupled  to  a  stepper  motor  and  provides  the 
independent  variable  in  the  experiment.  This  has  been  accomplished  i)y 
rebuilding  a  commercial  translation  stage  and  prism  table  with  much 
stiffer  springs,  developing  an  alignment  algorithm,  using  flexible  and 
sliding  coupling  mechanisms,  and  building  a  positioning  device  to  mare 
tne  micrometer  and  stepper  motor  axes  colinear.  A  significant 
advantage  also  is  that  the  lens  tends  to  reduce  the  magnitude  of  spot 
movement. 

A  device  which  was  needed  was  designed  and  built  here.  In  order 
to  get  small  spot  sizes  it  is  necessary  to  use  short  focal  length 
lenses.  This  causes  a  space  restriction  near  to  the  sample.  Mounts  to 
hold  small  lenses  near  the  sample  and  yet  have  the  necessary  degrees  of 
freedom  were  not  commercially  available.  A  commercial  system  was 
reduced  in  size  and  more  degrees  of  freedom  were  added  to  result  in  a 
positioning  device  1  inch  square  and  providing  5  degrees  of  freedom. 

In  order  to  facilitate  a  study  of  different  samples  and  the  use 
of  separate  materials  for  delay  line  calibration  and  actual 
measurements,  a  sample  mount  that  allows  the  easy  exchange  of  samples 
while  retaining  optical  alignment  by  always  bringing  the  front  optical 
surface  the  sample  into  contact  with  a  reference  plane  was  designed  and 
built.  The  reference  plane  can  then  be  adjusted  to  bring  its  normal 
parallel  to  a  rotational  axis  through  the  center  of  the  sample.  This 
rotational  axis  can  then  be  swept,  keeping  a  fixed  orientation  to  the 
sample,  to  give  optical  alignment. 

As  expected  we  have  found  the  signal  to  be  small.  Therefore 
the  use  of  some  type  of  signal  recovery  was  demanded.  Two  methods  have 
been  employed.  One  being  the  use  of  a  lock-in  amplifier.  A  chopper 
wheel  which  also  generates  a  timing  signal  was  built  to  modulate  the 
signal  at  a  known  frequency.  The  lock-in  then  only  amplified  signals 
which  had  the  same  specific  frequency  and  phase.  However,  the  lock-in 
detects  mainly  DC  signals  and  so  it  only  detected  the  average  of  one 
pulse  over  the  duty  cycle.  This  reduces  the  signal  level  by  a  factor 
of  lOU.  Depending  on  the  signal  strength  only  marginal  to  acceptable 
results  were  obtained. 

The  second  method  used  was  a  Boxcar  integrator.  The  aperture 
over  which  the  signal  is  integrated  was  held  fixed  in  time  with  respect 
to  the  pulsed  output  and  had  a  duration  of  10  ns.  With  a  pulse  width 
of  10  ns  and  needing  only  an  integral  of  the  pulse  this  method  should 
be  much  better  than  a  lock-in  since  it  is  not  affected  by  the  low  duty 
cycle.  However  the  sensitivity  of  the  Boxcar  is  much  worse  than  the 
lock-in  and  in  fact  was  not  sufficient.  A  preamp  for  the  Boxcar  is 
being  obtained  which  should  correct  this. 

Before  the  completion  of  the  stepper  motor  assembly,  a  few 
runs  for  data  acquisition  were  made  by  manually  controlling  the  delay 
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line  and  shutter.  Although  the  data  showed  the  expected  trends,  the 
scatter  was  too  large  for  quantitative  analysis.  The  actual  data 
taxing  has  so  far  been  beset  with  problems,  some  of  them  minor, 
requiring  modifications  and  rebuilding  of  machined  parts.  But  other 
problens  involving  breakdowns  in  the  laser  system  or  malfunctions  in 
the  detector  electronics  required  shipment  of  parts  to  the  manufacturer 
causing  delays  of  several  weeks  at  a  time. 

Besides  the  el  1 ipsonetric  technique  described  above  in  detail, 
otner  ’lethods  for  measuring  the  life  times  are  also  being  considered, 
^reparations  for  one  type  of  experiment  based  on  the  Haynes -Shockl ey 
techqniue  are  well  under  way. 

The  experiment  would  measure  the  transit  time  t  of  the  minority 
carriers  from  the  point  of  injection  to  the  collector.  The  excess 
minority  carrier  density  moves  with  a  velocity  v  between  the  emitter 
and  collector  which  are  a  distance  d  apart.  Hence  vt  =  d.  A 
homogeneous  field  E  is  assumed  between  the  emitter  and  the  collector. 

If  u  is  the  minority  carrier  drift  mobility  then  v  =  uEoru 

The  traditional  Haynes -Shockl ey  experiment  generates  minority  carriers 
by  applying  a  short  current  pulse  to  an  emitter  contact  whereas  in  the 
proposed  experiment  minority  carriers  are  generated  using  a  picosecond 
optical  pulse.  The  injection  point  and  hence  the  distance  to  the 
collector  can  be  varried  by  steering  the  beam  across  the  sample. 

The  sample  consists  of  a  rectangular  block  of  GaAs,  length  0.5 
cm;  width  -  25  pm  and  thickness  ~  5  pn.  Assuming  a  resistivity  of 
(.01)  )1  cm  a  resistance  of  4  is  obtained.  A  pulse  generator 
providing  200  volt  pulses  was  constructed.  This  will  provide  an 
electric  field  of  400  V  cm"^.  The  rationale  behind  preparing  samples 
ofas  small  a  width  and  thickness  possible  is  to  obtain  a  high 
resistance.  This,  along  with  a  pulsed  field  minimizes  the  heating 
effects  As  soon  as  techni strip-Au ,  an  etchant  for  gold  overlays 
existing  on  present  samples,  or  a  new  batch  of  samples  arrives  the 
experiment  will  be  performed. 

The  other  method  being  considered  involves  the  use  of  a  CW 
source  for  the  probe  beam.^  The  pulsed  pump  beam  and  the  CW  probe 
beam,  which  can  be  obtained  from  a  diode  laser,  hit  the  same  area  of 
the  sample.  The  reflectivity  of  the  material  would  change  as  a 
function  of  time  after  the  arrival  of  the  pump-pulse.  The  intensity  of 
the  reflected  probe  beam  follows  this  change  in  the  reflectivity  (Fig. 
6).  This  beam  is  mixed  with  a  portion  of  the  dye  laser  output  in  a 
nonlinear  crystal.  The  sum  frequency  light  generated  in  the  crystal 
can  be  measured  as  a  function  of  the  delay  between  beam  1  and  the  pump 
beam  to  obtain  information  about  the  decay  of  the  reflectivity  of  the 
semiconductor.  This  procedure  would  have  the  advantage  of  having  the 
delay  line  on  a  beam  that  does  not  have  to  be  focused  as  hard  as  in 
el lipsometry.  Thus  any  misalignment  of  the  delay  line  would  not  have 
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ds  pronounced  an  effect  on  the  data  as  observed  in  the  llipsometry 
arrangement.  Also,  if  the  CW  laser  is  chosen  to  have  energy  less  than 
the  band  gap  of  the  material,  this  experiment  can  easily  be  modified 
for  transmission  measurements.  Towards  performing  this  experiment, 
diode  laser  sources  for  the  CW  beam  are  being  evaluated  and 
calculations  to  find  the  most  suitable  non-linear  material  and  its 
orientation  for  the  cross-correlation  are  being  undertaken. 
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SECTION  VIII 


CONCLUSIONS 


8.1  Research  on  non-anneal ed  contacts 

The  aim  of  this  research  was  to  investigate  the  possible  method 
for  fabricating  non-all oyed,  low  resistance  (ohmic)  contacts  to  GaAs 
and  to  determine  whether  the  performance  of  such  contacts  could  be 
accurately  characterized  and  predicted  using  a  Schottky  barrier  model 
with  the  inclusion  of  electron  tunneling. 

The  initial  step  in  the  contact  fabrication  was  the  formation, 
by  high  temperature  diffusion,  of  a  heavily  doped  Sn  layer  in  the 
surface  of  an  p-type  GaAs  substrate.  Testing  of  the  characteristics  of 
this  layer  before  the  application  of  contact  metal  revealed  the  follow- 
i  ng 

1.  The  Sn  diffused  layer  can  be  accurately  character¬ 
ized  with  an  erfc  profile. 

2.  Use  of  this  erfc  profile  shape  in  the  theoretical 
prediction  of  average  layer  resistivity,  p  ,  with 
respect  to  surface  dopant  concentration,  Cj,  will 
give  results  which  can  be  experimentally  verified. 

However,  some  additional  consideration  is  needed  in 
assuming  a  significantly  lower  value  for  electron 
mobility  in  the  diffused  layer  than  would  be 
expected  for  bulk  material. 

3.  The  assumption  of  full  ionization  of  Sn  donors  is 
valid,  even  though  the  material  is  degenerate, 
apparently  because  of  the  effects  of  band  tailing 
and  impurity  band  formation. 

4.  At  the  diffusion  temperature  used,  950®  C.,  a 
maximum  Sn  surface  concentration  of  Cs  ”  2  x 
lO^^/cm^  was  observed,  implying  a  solid  solubility 
limit  for  Sn  in  GaAs  near  this  level.  This 
presents  a  direct  limitation  on  the  lowest  value  of 
contact  resistivity,  R^,  obtainable  for  a  contact 
made  to  this  doped  surface. 
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After  the  Au  contact  metal  was  applied  and  the  contact  structure 
delineated,  the  contact  was  tested  for  its  value  of  specific  contact 
resistivity,  R^,  using  the  transfer  length  method.  Samples  were  tested 
for  over  a  temperature  range  of  T  =  -15°  C  to  T  =  25°  C  to  determine 

the  Rc  response  to  temperature.  The  current  voltage  ( I-V) 
characteristics  of  the  contactwere  measured  and  used  as  base  data  for  a 
curve  fit  to  theoretical  I-V  equations  as  developed  from  a  Schottky 
barrier  model  for  the  contact.  The  theoretical  equations  gave  the  I-V 
response  for  the  contact  when  field  emission  (FE)  or  thermionic  field 
(TFE)  type  tunneling  was  present.  These  equations  also  formed  the 
basis  for  the  theoretical  equations  predicting  R(;  for  the  contact. 

The  above  investigat f on  revealed  the  following; 

1.  The  theoretical  equations  predicting  contact  TFE 

or  FE  I-V  response  can  be  accurately  fitted  to  the 
experimental  I-V  data  for  the  contacts  and  this 
procedure  yields  values  for  the  controlling 
parameters,  Cs,  and  barrier  height,  (pg.  The 
values  of  6$  so  obtained  were  verified  from  the 
previous  p  -Cs  measurements.  The  values  of  (jig 
were  larger  than  expected  and  varied  from  the 
usually  reported  constant  value  of  pg  =«  0.9  eV. 

The  higher  values  of  <i)g  are  attributed  to 
monolayers  of  interface  oxide. 

2.  Use  of  the  Cs  and  (jjg  values,  obtained  from  the  I-V 
curve  fit,  in  the  calculation  of  a  theoretical  R^ 
vs.  T  response  yielded  results  which  fairly  well 
matched  experimental  data.  This  means  that  the 
Schottky  barrier  model  can  adequately  characterize 
non-alloyed  contact  tunneling  behavior  and  can  be 
used  quantatively  in  the  study  of  such  contacts  to 
GaAs. 


8.2  Research  on  annealed  contacts 


The  above  procedure  was  also  applied  to  the  analysis  of  selected 
Au-Sn  doped  GaAs  contacts  and  Au-Ge-6aAs  contacts  which  were  given 
successive  anneals  from  250°  C  to  »■  425°  C  and  revealed  the  following: 

1.  Annealing  tests  on  the  Au-Sn  doped  GaAs  contacts 
confirmed  the  existence  of  a  thin  (»  lOA)  oxide 
layer  between  the  Au  and  GaAs.  This  layer  acts  to 
increase  the  effective  (>3  ,  and  therefore  raise  Rq 
but  does  not  substnatial ly  effect  the  analytical 
application  of  the  simple  Schottky  barrier  model. 
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2.  Annealing  the  Au-Sn  eloped  contacts  up  to  =>  300°  C 
gradually  removes  the  interfaclal  insulating  layer 
or  perhaps  makes  it  a  conductor. 


3.  Annealing  above  300°  C  results  in  the  formation  of 
small  iolated  current  conduction  areas  (conduction 
pinholes)  which  prevents  further  quantitative  use^ 
of  the  Schottky  model  for  accurate  determination  of 
Cg  and  (Jig  . 


4.  Results  from  the  annealed  Au-Ge-GaAs  contacts  also 
were  affected  by  the  apparent  formation  of  contact 
pinholes  but  it  seems  that  successive  annealings 
result  in  increasing  Cj,  due  to  the  doping  action 
of  the  Ge,  up  to  0$  -  10^^  /cm^  .  Results  on  (pg 
were  not  conclusive. 


5.  For  both  cases,  when  annealing  temperatures  above 
425  ~  450°  C  were  used,  Au  balling  up  and  contact 
surface  disruption  were  seen. 

6.  The  Schottky  tunneling  model  has  limited  usefulness 
in  the  analysis  or  prediction  of  contact  perform¬ 
ance  of  annealed  or  alloyed  contacts  unless  the 
processing  is  uniform  over  the  interface  area  of 
the  contact.  However,  the  trend  implied  from  the 
above  investigation  is  that  for  alloyed  contacts 
tunneling  is  the  major  phenomenon  involved  in 
causing  ohmic  performance.  Whether  or  not  the 
contact  behaves  as  predicted  by  the  simple  Schottky 
model  cannot  be  quantitatively  verified  with  the 
experimental  procedure  used  in  this  research. 


8.3  Implications  for  improvement  in  contact  performance 

The  verification,  by  this  research,  of  the  applicability  of  the 
Schottky  tunneling  model  in  the  analysis  of  nonanneal ed  contacts  to 
GaAs  allows  for  its  straightforward  use  in  trying  to  realize  contact 
improvement.  In  order  to  obtain  non-alloyed  contacts  with  lower  Rq, 
the  model  and  theoretical  equations  point  to  two  possibilities,  lowered 
4ig  or  increased  Cj. 

In  view  of  the  evidence  favoring  pinning  of  (j>g  ■»  0.9  eV, 
lowering  (^g  appears  unlikely.  However,  this  research  has  shown  that 
the  presence  of  interfacial  oxides,  on  the  order  of  only  10  A  ,  can 
effectively  increase  (pg  and  therefore  increase  R^.  The  implication  is 
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that  contact  processing  is  therefore  extremely  critical  in  the  case  of 
non-alloyed  contacts  to  avoid  the  formation  and  subsequent  detrimental 
effects  of  thin  interfacial  oxide  layers. 

If  a  <j>B  value  of  =■  0.9  eV  is  the  lowest  value  obtainable  then 
only  an  increase  in  0$  will  give  lower  values  of  for  non-all^ed 
contacts.  As  mentioned  before,  for  a  desired  range  of  R^  “  10"O  tJ-cm^, 
with  iJb  0.9  eV,  the  above  implies  that  Cg  on  the  order  of  lO^^/cm^  is 
needed.  Obtaining  this  large  of  a  value  by  normal  diffusion  or 
implantation  techniques  with  tin  seems  highly  improbable  because  of  the 
solid  solubility  limits.  Perhaps  a  combination  of  short  period  dopant 
diffusions  (predepositions)  or  implantations  and  subsequent  short 
period  surface  melting  and  regrowth  would  concentrate  the  dopant  into  a 
thin  surface  layer  with  a  Cj  above  the  equilibrium  solubility  limit. 

It  may  also  be  possible  to  find  other  donors  or  combinations  of 
donors  with  higher  solubility  limits.  Since  the  literature  on 
solubility  limits  in  GaAs  is  somewhat  limited  this  would  be  a  good  area 
for  further  research. 

Another  very  important,  but  experimentally  and  theoretically 
difficult  topic  for  further  research  is  the  nature  of  the  band 
structure  in  the  transition  region  between  the  metal  and  semiconductor. 
Even  for  very  moderate  annealing  temperatures  it  is  known  that  there  is 
not  a  sharp  transition  from  metal  to  gallium  arsenide,  but  rather  a 
complex  interpenetration  for  a  distance  of  at  least  a  few  attice 
constants.  For  such  a  transition  the  potential  barrier  would  probably 
not  be  sharply  peaked  and  under  some  conditions  might  be  reduced  to  a 
low  value  resulting  in  good  quality  contacts. 
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